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1.- INTRODUCTION 


Wind weighting, the adjustment of rocket launcher or impact parameters 
to accommodate prevailing winds, has for years been more of an art 
than a science. The main reasons are the non-steady characteristics 
of the wind profile and the complicated nature of the rocket response 
during the wind perturbation period. 

Conventionally, the problem of wind weighting on a rocket has been 
phrased in two ways: (1) for a given launch elevation and azimuth, 

evaluate the impact range and bearing* (2) for a given impact range 
and bearing, predict the required launch elevation and azimuth. 1 The 
former is straightforward. Given the physical and aerodynamic 
characteristics of the rocket and a wind profile, the equations of 
motion may be numerically integrated to achieve impact. The latter 
involves an iterative procedure, and therefore, several trajectories 
may have to be integrated. Iterations may or may not be convergent, 
and in many cases the procedure can become unacceptably time con- 
suming . 

In order to avoid integrating the equations of motions several times, 
the usual practice is to simulate the total wind effect on the com- 
plete trajectory rather than evaluating the continuous wind response 
along the path. This is achieved by defining a "ballistic wind 
velocity" and a "unit-wind effect", such that their product yields 
the effective displacement of the impact point due to wind. 

A computational procedure, SENS-5D, has been developed and pro- 
grammed in Fortran-Y to calculate the wind-weighting functions needed 
to evaluate ballistic wind' and unit-wind effects. Coriolis deflections 
and the first derivative of impact range with launch elevation are 
also calculated. If desired, ,a detailed time history of the coordinate 
and aerodynamic variables computed during integration of the equations 
of motion and a summary of salient variables at the final burnout, 
apogee and impact points (payload and spent stages) are available. 

These data may be used, with a known field, to solve the wind-weighting 
problem in either of its two forms. 
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■ 2. ASSUMPTIONS 


For simple and accurate simulation, this program utilizes the 5 
degree-of “freedom rigid-body dynamics model from first stage ignition 
to the last stage burnout. After the last burnout, the .program 
shifts to the 3 degree-of -freedom point mass dynamics . 

To save time the program assumes linear aerodynamics, i.e., the 
aerodynamics coefficients are functions of Mach number only. 

For stability concern, . both the aerodynamic damping and the jet 
damping terms are used for angular motions. 

Only axially symmetric rocket is considered. No attempt has been 
made to include the effect of thrust or fin misalignment and center- 
of-gravity deviation from the symmetric axis. 

Earth model is a rotating sphere. A 1962 standard atmosphere has 
'been chosen for atmospheric model. 
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3. COORD INS!TES- AND TRANSFORMATIONS^ 


3»1 COORDINATES SYSTEMS 

The coordinate systems used in this trajectory simulation are -defined 
below. 


Inertial Coordinate System (X^ Z) t 
Origin — at earth center . 

X - on the earth equatorial plane, pointing to the zero 
longitude at launch . 

Y - on the earth equatorial plane, pointing to the 90° 
longitude at launch, 

Z - perpendicular to the equatorial plane, pointing to the 
North Pole. 

Earth Coordinate System (X^, Z^,) • 


Origin — at earth center, 

Xg - on the earth equatorial plane, always pointing to the 
Greenwich longitude. 

Y^ - on the earth equatorial plane, always pointing to the 
90° longitude. 

Z„ - perpendicular to the equatorial plane, pointing to the 
North Pole. 

Instantaneous .Topo.c.entri.c Coordinate Systejn, .(x^, y^, z,j) ; 


Origin — at the projection point of the moving rocket on 
earth surface. 

x„ - on the local horizon plane tangent to the instantaneous 
projection point of the rocket, directed along the local 
geocentric north. 

y„ - on the local horizon plane tangent to the instantaneous 
projection point of the rocket, directed along the local 
geocentric east. 

Zrn - perpendicular to the instantaneous local tangent plane, 
directed along the geocentric radius vector and pointing 
toward the earth center. 
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Body Axis Coordinate System (x^, y^, ^g) : 


Origin — at the center of gravity of the rocket. 

Xg “ along the rocket principle (longitudinal) axis, positive 
forward, 

Yg “ normal to the Xg - Zg symmetric plane, completing the 
right-hand system. 

Zg - in the principle plane of symmetry of the rocket, 

perpendicular to the Xg axis and positive downward. 


3 . 2 TRANSFORMATION 

. Transformation from Body Coordinate System to Instantaneous 
Topocentric Coordinate System 


’^t' 


’^b' 

Yt 

[^] B->T 

^B 

_^T 


/b_ 



where 

B^-T 



translation matrix 


(3-1) 


(3-2) 


Cos9Cos\^ Sin4»Sin0Cos^-Cos^Sin!j; Cos(J>SineCosii) +Sin(j)Sinij; 
Cos0Sini|^ Sin(j)Sin0Sini|j+Cos(fiCosTj^ Cos(j)Sin0Sini|) -Sini^Cosi|; 
-Sin0 Sint()Cos0 Cos(J)Cos0 

= rotation matrix (3-3) 

0 - Euler angle for pitch 
- Euler angle for yaw 
<j) - Euler angle for roll 
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However, for the present 5 -degree -of -freedom case, the roll angle is 
always equal to zero. Thus, the rotation matrix can be simplified to 


[h] b.. 


Cos0Cosi|^-Sini|; Sin6Cosi|J 
Cos0Sinijj Cosi|J Sin0Sinip 
-Sin0 0 Cos0 


(3-4) 


. Transformation from Instantaneous Topocentric Coordinate System 
to Inertial Coordinate System 



X 






Y 



Yrj, 

^ ["] T.I 


Z 



_^T_ 

• 

where 


r* 





,-SinyCosX -SinX -CosyCosX 


R 

T-+I 

-SinySinX CosX -CosySinX 

' 


Cosy 


0 -Siny 



rotation matrix 
CosyCosX 
Rg CosySinX 
Rg Siny 


translation matrix 


(3-5)'' 




(3-6) 


(3-7) 


y - latitude of the Instantaneous Topocentric Coordinates 
origin, 

X - longitude of the Instantaneous Topocentric Coordinates 
origin with respect to the rotating Greenwich longitude 

X - longitude of the Instantaneous Topocentric Coordinates 
origin- with respect to the inertial zero longitude at 
launch , 

X = X + Wgt 

o)g - angular velocity of the earth rotation. 
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, Transformation from Body Coordinate System to Inertial Coordinate 
System 



where 



^11 ^12 ^13 

^21 ^22 ^23 
^31 ^32 ^33 


= rotation matrix 


- -SinyCosXCos 8 Cosi|> -SinXCos0Sinti» +CosyCosXsin 6 
Ri 2 = SinuCosXsinij;~SinXCost|J 

= -Sin]iCosXsin9Cosiji “SinXsin0Sint|» -CosyCosXCos 0 
^21 ” '-SinySinXCosOCoSTjj +CosXCos0Simj^ +CosySinXsin0 
R 22 = SinySinXsini|^+CosXCosT|^ 

R 23 = -SinySinXsinOCos^' +CosXSin 0 Sint|; -CosySinXCosO 

R^j^ = CosiaCos 6 Cosi|) +SinySin9 

R ^2 = -CospSint|; 

r ^2 - CosijSin0Cost|) -SinviCos0 


(3-8) 


(3-9) • 


NOTE: Here again the Euler angle of roll (■(}>) has been 

eliminated from the transformation components. 


[’^J B-I 


R CosyCosA 
R CosySinX 
R Siny 


= translation matrix 


(3-10) 


R = (Rg + h) = (X^ + y^' + 


(3-11) 
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3.3 BODY AXES DIRECTION COSINE VARIATIONS 


The angles of the direction cosine are sometimes undefined, e.g., 
the heading azimuth becomes indeterminate if the flight elevation 
reaches 90°, Thus, the elements of the rotation matrix (e.g, 3-9) 
are calculated from the orientation angles only during problem starts . 
Afterward, the elements of the rotation matrix are determined by 
integrating the time rate of change of the direction cosine elements . 

The time rate of change of the rotation matrix elements can be cal- 
culated from pitch (P) , yaw (Y) , and spin (S) as follows: 


• • * 





%1 ^12 ^13 

* 


R^l R ^2 %3 


0 -Y P 

^21 ^22 ^23 


^21 ^22 ^23 


Y 0 -S 

^31 ^32 ^33 


^31 ^32 ^33 


-P S 0 


Since three of the nine rotation matrix elements depend on the re- 
maining six, it is sufficient to integrate only six rotation matrix 
§leme]^ts a^d calculate the remaining three from the relation 


^13 


^21 ^ 32 . " ^31 ^22 

^23 

= 

^31 ^12 “ ^11 ^32 

^33 


^11 ^22 "■ ^21 ^12 


(3-13) 



4.- ROCKET MOTION- RELATIVE TO THE AIR 


4.1 WIND -profile 

The wind profile data is usually gathered and made available to the 
range by on-site or nearby weather station. The composite wind data 
are smoothed and edited into the format of . local north wind component 
(Wjj) and local east wind component (Wr) vs. altitude (h) . Due to its 
random nature, the vertical drift- or gust of the wind is totally 
neglected in the present simulation. 

Since the wind data are in local topocentric coordinates, transfor- 
mation must be performed in order to re-orient its components into the 
body axes coordinates . 

Assume [uw, vw/ are the wind velocity components in the body axes 

coordinates , then 




1 


T->B 

Wg + Ug R Cosy 

_V 


0 


where 



t^-b 


rotation matrix from topocentric coordinates 
to body axes coordinates 


(4-1) 


4.2 Combination of the wind to the rocket motion 

The rocket velocity relative to the true air is equal to the rocket 
absolute velocity relative to earth plus the wind velocity, i.e. 


'^ra 


1 

+ 

d 

1 

v 


V + Vr-r 

RA 


W 

_^RA 


+ 

1 
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where (u, v, w) are the velocity component in body axes coordinates. 
Finally, the total velocity relative to the true air becomes 


V 


RA 


- [u. 


RA 


+ V. 


RA 


^ 2 . 1/2 

+ «EA I ■ 


( 4 - 3 ) 


and 


Dynamic Pressure = q = 


Mach No. 


M 


V, 


RA 


a 



2 


( 4 - 4 ) 

( 4 - 5 ) 
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5. AERODYNAMICS 


In SENS-5D forces are assumed to vary linearly with the aerodynamic 
angles arising from vehicle motions relative to the air mass. 

These angles are defined asj 

Angle of attack = a = Tan 

Angle of sideslip e 3 = T 

The resulting forces are* 

Drag (force along X-axis) = D = C^qs 


'^RA'^^RA 


Fw^raI 


(5-1) 


Normal force due to a (along Z-axis) e 

Normal force due to 3 (along Y-axis) E 

Jet damping force (along Z-axis) = F. 

3 z 

Jet damping force (along Y-axis) e f. 


^3 ^ ^ 3 *^® 

E mL^^P 


(5-2) 


-Because an axially symmetrical vehicle is assumed, we have 
°Na “ ‘^NS " '^0,6 “ “ ” “• 

The moments arising from these forces are; 

Pitching moment (about Y axis) = M = N L _ 

^ y -.a sm 

Yawing moment (about Z axis) EM = N„Ii 
^ z 3 sm 

Pitch damping moment (about Y axis) E (Pd/2Vj^) qsd (5-3) 

Yaw damping moment (about Z axis) E = C (Yd/2V„j. ) qsd 
Jet damping moment (about Y axis) e m^^. = 


Jet damping moment (about Z 


axis) -= M._ 
J ^ 
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In the above equations • 


Again due to syininetry, C^p = 

q = dynamic pressure 
s = reference area 
d = reference diameter 
m = time rate of change of mass 

L = static margin (5-4) 

sm ^ 

distance from nozzle exit plane to c.g. 

Cp = drag or axial force coefficient 

C„ = normal force coefficient 
N 

Cmp= pitch damping force coefficient 

P = pitch rate (angular velocity) 

Y = yaw rate (angular velocity) 
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6. EQUATIONS OF MOTIONS 


For the purpose of achieving both time saving aiid sufficient accuracy, 
the present trajectory simulation employs 5-D rigid-body dynamics 
from the first ignition to the last burnout. After the last ‘burnout, 
3-D point-mass dynamics is adopted. 


6.1 5-D RIGID BODY DYNAMICS 

6.1.1 Equations of Linear Motions 

The equations of linear motion are written in the inertial system. 
Let - ' 

r = + Y^ + 


be the distance of the vehicle from the center of the earth. The 
linear motions are then given by ■ ' ' 



where 




E 

m 


(6-1) 


GM 



B-vl 


product of universal gravitational constant and mass of 
the earth , 

rotation matrix from body to inertial system (see 3-9) , 


m 

E 


total mass of the vehicle, 

external force on the vehicle in body system. 


Th - D 




N 

a 



( 6 - 2 ) 


with Th and D as thrust and drag forces respectively, and other 
terms are defined in Equation 5-2. 


6.1.2 Equations of Angular Motions 

The equations of angular motions are written in body system. It is 
assumed that rocket is axially symmetric and that roll is neglected. 
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Let ^(P,Y) be the angular velocity vector about Y and Z axes 
(i.e., pitch, and yaw) respectively. Then the angular motions are 
given by 

_ S ■ 

at2- I ' 


where 

I - moment of inertia about Y or Z axis, 

-y 

M = external moment (torque) on the vehicle, 


with the terms above defined in Equation 5-3, 




z mz j z 


6.2 3-D POINT MASS DYNAMICS 


The equations of linear motions are written in inertial system. 

No angular motions ..are considered. The equations are given as in 
(6-1) , but external forces are now computed in topocentric system, 

i .e. 




R 


£ 


T->-I 


m 


(6-4) 


where 

T-j-i “ rotation matrix from topocentric to inertial system, 
■■ and 


E 


(Th-D) Cos e Cos a 
(Th-D) Cos e Sin a 
-(Th-D) Sin e 


(6-5) 


with e and a as flight path angles (elevation and azimuth respectively.) 
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7. NUMERICAL' INTEGRATION 


3 4 

7.1 PREDICTOR-CORRECTOR METHOD ' 

The integration, technique is Adams -Bashforth ' s fourth-order,, modified 
predictor-corrector method for the solution of general initial-value 
problems . 


Given k first-order differential equations 



fi (x,y^) ; i = 1, k 


(7-1) 


and values of y, and therefore f(x,y), at four equally- spaced inter- 
vals (h) of X, the result for y at the fifth point of x is obtained 
by the following scheme: 

Predictor: '55fj^-59£j^_^+37fj^_2-9f^_3l , (7-2) 


.. -(0) (0) 251 , (0) , 

Modifier: - 27TT <^.1 ' ’ 


Corrector: =• y„+^ (x^.^3^,?^»[)+19£^-5f^_3+f^_2] , 


(7-3) 

(7-4) 


Final Value: = y^=> + (y^"[ - y^=>) (7-5 

AjDOve the subscript i is dropped because (7-2) to (7-5) hold for 
any i . 

Expressions (7-2) and (7-4) have local truncation errors. Errors 
are proportional to h 5 times the fifth derivative of y. Assuming 
that this fifth derivative does not vary appreciably in the interval, 
they are given by the second terms in (7-3) and (7-5) for the pre- 
dictor and for the corrector respectively. The second term in (7-5) 
is, therefore, a measure of the truncation error in yn+1/ and there- 
fore control of accuracy and adjus'tment of step size h can be done 
by generating the following test value- 


6 


1 k 

^ i=l '"i 



,(c) 

n+1 


(7-6) 


where Wj^ are error weights specified in the input. 
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If 5 is greater than the upper error tolerance, step size h is 
halved. If 6 is less than the lower error tolerance/ step size is 
doubled. 


7.2 RUNGE-KUTTA METHOD 

Starting the predictor -cor rector method requires the functional and 
derivative values at four equidistant points, yo .and fCxQ^yo) are 
given initial conditions. - For computation of Yi, f(xi,yi), 

Y2, t Y3 adjustment of the step size h 

to accuracy requirements, a special fourth-order Runge-Kutta pro- 
cedure suggested by Ralstori*^ is used. 

For initial control of accuracy and adjustment of step size h in 
starting the integration, a result for y^^^= y(xo+h) is computed 

( 2 ) 

using step size h, and then another result = y(x 
puted using step size h/2 twice. If the relative error, defined by 


,h.h. . 

3 + 2 + 2 -) xs com- 




(7-7) 


-5 (2) 

is less than 10 , then solution y 2 is assumed correct. If 

5>10“5^ step size is halved, and the procedure starts again at point 
xq . In case step size is less than a given minimum, the last 
solution y^ 2 )^oul <3 assumed correct. 


See Reference 3, Page 200, Expression (5.6-49) 
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8-. WIND-WEIGHTING FUNCTIONS 


1 ' 


8.1 BALLISTIC WIND FUNCTION 

Ballistic wind, velocity W is a hypothetical wind which is constant 
in direction and magnitude from the ground level to- a defined upper 
limit of the effective atiposphere SENS-5D, = 100000 ft. 

The integrated effect of W on the rocket impact is equal to the 
integrated effected of the actual wind. 

Ballistic wind is defined by 

= s' d.W. ; (8-1) 

i=l ^ ^ 


,N 

S 

i=l 


d. 


1 


1, z = z„, 

' max N ' 


( 8 - 2 ) 


where W^ is the wind velocity at altitude levels i = 1,2,3, N, 

and di are the weights defined by 



(8-3) 


such that li is the impact vector with winds between ground and 
the ith altitude levels only, d^'s are sometimes called wind- 
weighting fvinction, so called delta-f(z) curves. 

Another wind-weighting function, which is of interest is f(z) curve 
defined by 


f. = 
1 



0 



(8-4) 


8.2 UNIT-WIND EFFECTS 

The unit-wind effect 6(e) for a given launch elevation e, is the 
magnitude of the impact displacement vector due to a unit ballistic 
wind to the height Zj^ax* Assuming that the rocket response to wind 
is linear, the impact displacement vector due to the winds is 
6(e)W, The definition of 6(e) is extended to provide for head, tail 
and cros s winds . 
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9. INPUT DESCRIPTION 


As with any aerodynamics program, the most difficult problem for 
the user is understanding and preparing the voluminous input required. 
We have attempted to simplify the preparation of the input data. 

The input data are split into two sets. Set 1 is a series of four 
program control lists, which use the FORTRAN ‘NAMELIST* facility. 

Set 2 encompasses the aerodynamic tables . The data is read with 
FORMAT (V) , a list directed read operation which may be peculiar to 
Honeywell Series 6000 FORTRAN. 


SET 1 


All the entries in Set 1 are entered in ‘NAMELIST? foannat, on file 
code 05. 

Set 1 consists of four lists, all of which must be present. The 
lists are* DLIST^ . BLIST^ FLIST and ULIST, DLIST may be empty if one 
desires to keep all of its default values, ■ The other lists may be 
empty only if they are to be bypassed. All the lists have default 
values defining some or all the input variables. They may be over- 
ridden by re-defining them. 

Note that the units used in this program are feet, pounds, seconds, 
and degrees . 

Definition of the Variables in DLIST : 


The format of DLIST is • 


^LIST 

in 

Colimns 2 

‘ through 7, and then 





GDLATL 


Geodetic 

la'^itude of launcher 


. 



LONGL 

= 

Geodetic 

longitude of launcher 





T2ERO 

= 

Initial 

time 





AZERO 

= 

Initial 

altitude 





VZERO 

= 

Initial 

velocity 





ITAPEl 

— 

Number of the file code used for the 
This is a dummy entry. The program 
read from file code 05, 

input 

forces 

of Set 1. 
Set 1 to be 

ITAPE2 


Number of the file code used for 
Set 2. 

the 

input of 

data 
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ITAPE3 = Number of the file code for storing thrust-mach tables 
after cross -interpolation, to get one time table and one 
mach table. 

ITAPE4 = Nximber of the file code for storing data to be used for 
spent stage calculations . ■ • • 

DTLA = Step size to be used from launch through apogee 

DTAI = step size to be used .from apogee through impact 

DMIN = Minimum step size allowed 

DMAX = Maximum step size allowed 

EPTINY = Minimum error allowed during integrations 

EPBIG = Maximum error allowed during integration 

ERW = Vector of error weights used in the predictor-corrector 

VOLIM = Minim\im initial velocity' 

VLIM = If current' velocity < VLIM, apogee is assumed. 

HLIM = If current altitude (after apogee) < HLIM, impact is 

assumed . 

ORDER = Array of leng-th 14 storing -the orders of polynomials used 
in interpolations to create one time table and one mach 
table . 

IPUNCH = Punch-out option. 

= 0 No punch-out. 
f 0 Will be punched-out. 

$ after the last item in -the list. 

Definition of the Variables in BLIST 

BLIST controls to option to calculate burnout flight elevation, 
apogee altitude, spent-stage impact, and payload impact range as a 
function of launch elevation angle and payload weight. 

If "this option is not required, enter BLIST as an empty list with 
^BLIST in Columns 2 through 7 and ^ in Column 10 . 
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The format of a non-empty BLIST is: 


^BLIST 

NPL 

PLM 

AZGDL 

NANG 

ANG 

WIND, 

WNAZ, 

ALOW 

AHIGH 


NLEV, 
ALTW, 
SPEED , 
DIR 


IROT 


IPRINT 


in Columns 2 through 7, and then 
= Number of payloads ■ 

= Array to store NPL payloads 
= Launch azimuth (geodetic) 

= Number of launch elevation angles 
= Array of NANG elevation angles (geodetic) 


= These four variables are used to define a given constant 
wind from one altitude to another, and zero elsewhere, 

WIND = Wind speed (always positive) 

WNAZ = wind azimuth, measured clockwise from North, 

(North-to-South or' East-to-West is positive wind) 

ALOW = Minimum altitude of the wind 

AHIGH= Maximum altitude of the wind 


= These four variables are used for defining a variable ' 
wind, i.e. a wind which varies with respect to altitude 

NLEV = Number of altitude levels 

ALTW = Array of NLEV altitudes 

SPEED = Array of wind speeds (magnitude of wind velocity) 

DIR = Array of wind direction (measured clockwise from 
local North) 

= =0: non-rotating earth 

7 ^ 0: rotating earth 

= Print-out option 

- 0: summary only 

= 1* detail trajectory with summary 
= 2* debug option 
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JSPENT 


= Option for calculating spent stage (s) trajectories 

= 0 : no calculation 
= 01 will be calculated 

0 after the last item in the list. 


Definition' of the Variables in FLIST 

PLIST controls the option to calculate -the wind-weighting factor as 
a fiinction of altitude. 

If this option is not desired, FLIST must still be present as an 
empty list with 0FLIST in Columns 2 -through 7, and 0 in Column 10. 

The format of a non-empty FLIST is; 

0FLIST in Columns 2 -through 7, and -then 

WPL = Payload weight 

AZGDL = Launch azimuth (geodetic) 

ELGDL = Laxanch elevation (geodetic) 

WIND' = Wind speed (always positive) 

WNAZ = Wind azimuth (measured clockwise from North, while 
North-to-South or East-to-West is positive) 

NLEV = Number of wind altitudes 

ALEV = Array of NLEV altitudes 

IROT = 0 ; non-rotating earth 

^ 0 1 rotating earth 

IPRINT - =0: summary only 

— 1; detail trajectory with summary 
=2; debug option 

0 after the last item in the list. 


Definition of the Variables in ULIST 

ULIST controls -hhe option to calculate the coriolis deflections due 
Nor-th and East, range derivatives and unit-wind effects for head, 
tail and cross winds as a function of launch elevation angle. 

As in lists BLIST and FL]|sT, ULIST may be empty if this option is 
not desired. For an empty list, put 0ULIST in columns 2 through 7 
and 0 in Col-umn 10 . 
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The fo 2 nnat of a non-empty ULIST is: 
^ULIST in Columns 2 through 7, and then 


WPL 

= 

Payload weight 

AZGDL 

- 

Launch azimuth (geodetic) 

NANG 

= 

Number of launch elevations 

ANG 


Array of NANG launch elevations (geodetic) 

WIND 

= 

Wind speed (always positive) 

ALOW 

= 

Minimum altitude 

AHIGH 

= 

Maximum altitude 

IROT 

= 

=0; non -rotating earth 
7 ^ 0; rotating earth 

IPRINT 

= 

=0; summary only 

= 1* detail trajectory with summary 
=2* debug option 


0 after the last item in the list. 


Default Values 

Unless otherwise defined, the following values for the variables in 
NAMELIST data will be assumed by the program. 


NAMELIST 

^DLIST 


VARIABLE 



GDLATL 

SS 

37.-848 

LONGL 

= 

-75.4736 

TZERO 

= 

0 

AZERO 

ss 

0 

VZERO 

= 

0 

I TAPE 1 

= 

5 - Thi 

I TAPE 2 


by 

1 

ITAPE3 

= 

2 

I TAPE 4 

= 

3 

DTLA 

= 

0.01 

DTAI 

= 

1.0 

DMIN 

= 

0.0002 

DMAX 


20 .0 

EPTINY 

= 

0.0001 

EPBIG 

= 

0.001 

ERW 

= 

All of 

VOLIM 

= 

0.1 

VLIM 

= 

40.0 
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NAMELIST 


VARIABLE 


^DLIST (continued) HLIM = 

ORDER 

IPUNCH 

^BLIST WIND = 

WNAZ 

alow 

AHIGH 

NLEV 

IROT = 

IPRINT 

JSPENT 

^FLIST WIND 

WNAZ 

NLEV 

ALEV 


IROT 

IPRINT ~ 

^ULIST WIND ■ 

ALOW 

AHIGH 

IROT 

IPRINT 


20.0 

All of them = 1 
0 

0 

0 

0 

100 ,000 
0 
1 
0 
0 

20.0 

0 

12 

0, 50, lOO; 200, .500, 1000, 
2000, 3000, 5000, 10000, 
50000, 100000 
0 
0 

20.0 

0 

100,000 

0 

0 


SET 2 

Vehicle Data 


These data are read from file code ITAPE2 (=1 unless redefined in 
DLIST of SET 1) using "FORMAT (V) " . This is a list directed read which 
may he peculiar to Honeywell Series 6000 Fortran. Data are entered 
on cards using columns one through 72 (1-72) , and are separated by 
commas, spaces, or by beginning a new card (see pitfall 2 below). 
Decimal points should be used in real numbers, deleted in integers. 

All units are feet, pounds, degrees or seconds. 

SET 2 begins with a title card and is divided into a block of general 
data and a series of phases. A new phase is required whenever a change 
in aerodynamics , vehicle geometry or a step change in weight occurs . 
Each phase is preceded by its own title card. Within the general block 
and each phase, data are presented in a specific order as groups or 
tables. Each group or table is also preceded by a title card. The con- 
tent of the title cards is left to the user, but their presence is 
mandatory . 
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in the general data block, the first group contains the total rocket 
weight (less payload) • followed by the time at which point mass 
(rather than 5-D) integration is to begin (normally the start of the 
final phase) , This group is followed by three tables in exactly this 
order , 

1. phase start times (referenced to zero) 

2. time for spent stage separations (these must also be phase 
start times) 

3. discarded weight of spent stages (one for each time in Table 2) 

The first value in each table must be the number of data points 
(an integer) in the table. Each table must have at least one entry, 
zero if no other value is appropriate. If the limits of an indepen- 
dent table are exceeded, the last corresponding dependent values 
are used in subsequent calculations . 

In each phase the first group of data consists of four vehicle di- 
mensions. They are nozzle exit area, length of the vehicle, reference 
area, and reference length. It is followed by fourteen (14) tables in 
exactly this order. 

1. Time (since phase start) for thrust and weight; 

2. Thrust (one value for each time in Table 1. Vacuum reference is 
assumed unless the word "SEA" appears in the title card for this 
table) ; 

3. Expendable weight remaining in this stage (one value for each 
time in Table 1) ; 

4. ' Mach nvimber for drag coefficients; 

5. Drag coefficients (one for each entry in Table 4); 

6. Time (since phase start) for gravimetrics ; 

7. Distance from nose to center of gravity (one for each time in 
Table 6) 

8. Pitch moment of inertia (one for each time in Table 6); 

9 . Mach number for pitch damping coefficient; 

10. Pitch damping coefficient (one value, entered as a positive 
mamber, for each entry in Table 9) ; 

11. Mach number for slope of normal force coefficient; 
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12. Slope of normal force coefficient per unit angle of attack 
(one for each, entry in Table 11 ) • 

13. Mach number for center of pressure location; 

14. Distance from nose to center of pressure (one for each entary 
in Table 13) . 

The rules for table entries in the general block also apply to the 
phase tables. If the limits of an independent table are exceeded, 
the last corresponding dependent values are used in subsequent 
calculations. Data are not saved from phase to- phase. 

Input of SET 2 is greatly facilitated by the use of preprinted foinns 
containing the title cards and spaces to record data. A set of forms 
prepared' for use by the Ground and Flight Safety Section, SQAEB, 

NASA Wallops Flight Center is shown in Appendix A. The user may find 
this input clearer after scanning these forms and the sample input 
for a Nike Cajun shown in Appendix B. 


Pitfalls to Avoid 

1. It is a rule of "FORMAT (v)" that only one delimiter (comma, 
blank, etc) appear between values in tables.. Particularly, lines 
should not end with a comma, and a given value may not be con- 
tinued on another card. 

2. A zero level of thrust is always taken with respect to a vacuum. 
Therefore, during coasting phases, the thrust table (input as a 
single value, zero) must not be titled with the word "SEA", or 
errant forces may be calculated for nozzle exit corrections. 

3. SET 2, each phase, and each data group or table must be immediately 
preceded by a title card and must appear in' the order specified. 

4. In SET 1, all four namelists must be present, and each must be 
closed with a dollar sign ($) . Examples of empty lists appear 
in Appendix A and of used lists in Appendix B. 
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10. OUTPUT DESCRIPTION 


Appendix C shows the output which is produced when the input set of 
Appendix B is used. 

Page 1 is an echo of the Set No. 1 input lists. 

Pages 2 through 5 are an echo of the Set No. 2 tables. 

Page 6 gives a summary for the reference trajectories. ' 

Page 7 gives a summary of the trajectory of the spent stage. 

Page 8 lists the F(Z) and DF{Z) data for the given wind condition. 

Pages 9 and 10 list the unit-wind effects, range derivatives, and 
Coriolis deflections both in English units and in MRS units (which 
are often required for off -range operations) . 

Page 11 is an echo of the punch output of the SENS-5D data which 
may be used in wind weighting. Reader should note IPUNCH = 1 in 
^DLIST of input data (see Appendix B) . 

Pages 12 through 19 list a detailed trajectory which may be. useful 
in comparing the SENS-5D trajectory with another trajectory pro- 
vided by a range user. 
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XI. SUBROUTINE DESCRIPTION 


APCAL: 

APSTEP: 

ATMS PH: 

ATTERP : 

AZRAN: 

BAI; 

BOCAL • 
EFFECT: 

FZCURV: 

IPCAL: 

IPSTEP: 

LINT: 

LIST: 

PCM: 


Calculates current time, altitude, surface azimuth, and 
surface range at the apogee .point. 

Calculates the time that will be taken from current position 
to the apogee point. This is used in adjusting the step 
size until the apogee is approached. 

Computes temperature, pressure, density, and the speed of 
sound using an eighth degree polynominal. approximation 
for the 1962 Standard Atmospheric Tables. 

Calculates the indices within which the current altitude 
falls in the altitude tables (called by ATMSPH) . 

Calculates current. latitude , longitude, surface azimuth, 
and surface range. 

Integrates the trajectory, for a given set of NANG eleva- 
tion angles and a set of NPL payloads , Prints out salient 
features at burnout, apogee, and impact points. 

Calculates current time, altitude, surface range and 
azimuth, and relative velocity at the final burnout. 

Calculates , for a given payload and a given set of NANG 
elevation angles, unit-wind effects for head, tail, and 
cross winds, coriolis deflection, and tower tilt. 

Calculates for a given payload and elevation angle the 
F{Z) curve. Also calculates the ballistic wind factors 
defined by DF(Z) = F(Z) - F(Z-l), 

Calculates current time, altitude, surface azimuth and 
range at the impact. 

Calculates the time that will be taken from the current 
position to the impact point. This is used in adjusting 
the step size until the impact is approached. 

Program for linear interpolation using Aitken's Iteration 
method employing a polynomial of degree up to 10, This 
can be used for extrapolation, but with caution. 

Prints input data. 

Knowing values at a set of four points , it integrates at ■ 
the fifth point (step size DT) , Procedure is modified 
Adams -Bashforth Predictor -Corrector method. If error 
> EPBIG, step size is decreased by a factor of 2; if 
error < EPTINY, step size is increased by a factor of 2. 

It determines if the phase will change after DT, and up- 
dates the data accordingly. (References 3 and 4) 
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RDATA: 


REFAZE . 


RK4: 


ROCKET; 

RTERP ; 
RTIMAT: 

SETIV: 

SOLV3D; 

SOLV5D : 

SPENT; 
SPLASH; 
TPRINT : 
WCONST ; 


WTERP: 


Reads from disc number ITAPE2 the thrust and drag data, 
interpolates for one time table and one mach number table, 
and stores the data on disc number ITAPE3. 

Updates the thrust-drag data, the integration variables, 
and mass at the phase change. Also stores integration 
variables for spent stages. 

Runge-Kutta method of order 4. At launch and whenever the 
step size changes, this subroutine is called. There is 
an option to integrate one step or three steps, (Reference 
3) 

Interpolates aerodynamic variables employing tables created 
in RDATA at current time and Mach No. 

Calculates indices to be used in ROCKET, 

Calculates at the current time the rotation matrix from 
topocentric to inertial axes , Also calculates altitude , 

Initializes the various physical parameters necessary to 
start the integration at. launch. 

Solves for functions P in equation Y* = f in 3D after final 
burnout , 

Solves for functions F in equation Y’ = f in 5D from launch 
to final burnout. 

Calculates trajectory for each spent stage. 

This is similar to BAI, but is called by EFFECT and FZCURV. 

Prints in detail the salient features along the trajectory. 

Constant wind table. Given the input parameters (WIND, 
WNAZ,-ALOW; AHIGH) , it creates table; ALTW(l) = ALOW, 

ALTW(2) = ALOW+1, ALTW(3) = AHIGH? VELN(l) = VELE (1) = 0, 
VELN(2) = WIND x COS(WNAZ+PI) , VELE (2) = WIND X SIN(WNAZ+PI) , 
VELE (3) = VELE (2) , and VELN(3) = VELN(2) . 

Interpolates for wind velocities at current altitude. 
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ACC. 

AD: 

AHIGH : 

ALEV; 

ALOW: 

ALT; 

ALTW: 

ANG: 

APA: 

APAZ : 

APR; 

APREV: 

APT; 


12. PROGRAM VARIABLE DEFINTTIONS 

Rocket thrust and mach niomber tables, A (I,J) . I is the 
element number. J has the following meanings; 

J = 1: Time array for thrust and propellant weight. 

= .2; Thrust table. 

= 3* Propellant weight table. 

= 4; Mach number array for drag coefficient. 

= 5; Axial drag coefficient. 

= 6: Time array for C.G. position (from nose) and 

moment of inertia, 

= 7; Distance from nose to C.G, 

= 8; Pitch moment of inertia, 

— 9; Mach number table for pitch damping coefficient, 
=10; Pitch damping coefficient 

= 11; Mach number table for slope of normal force co- 
efficient, 

=12* Slope of normal force coefficient, 

= 13; Mach niomber table for nose to c.p, distance, 

= 14; Nose to center-of-pressure distance. 

Relative acceleration in inertial coordinate system. 

I 

Axial drag, j 

Maximum altitude up to which wind is used and beyond which 
is zero. 

Altitudes of wind strata for F(Z) curve. 

Minimum altitude below which wind is zero. 

Current altitude. 

Altitude table for winds . 

Array for launch elevation angles. 

Altitude at the apogee point. 

Surface azimuth below the apogee point. 

Surface range below the apogee point . 

Previous altitude for interpolation among the wind altitude 
table . 

Current time at the apogee. 
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ARRAY: Dumiay array to store data or variables. 

AZ: Current azimuth (geocentric). 

AZGD: Current azimuth (geodetic) , 

AZGDL: Launch azimuth (geodetic) , 

B; Array of length 4. In B(J ) , J stands for* 

J = 1; nozzle exit area 

= 2* distance from nose to nozzle exit 
= 3: aerodynamic reference area 

= 4: reference diameter 

BOA: Altitude at the burnout point. 

BOAZ* Surface azimuth at burnout. 

BOEL*. Elevation angle at burnout. 

BOR: Surface range at burnout. 

BOT: Current time at burnout. 

BOVL: Relative velocity (topocentric) at burnout. 

CA: Drag coefficient. 

CF: .Array of conversion factors. Multiplying by CF(J) means: 

J = 1: radian to degree 

= 2: degree to radian 

= 3; nautical miles to feet • 

= 4* fee't to nautical miles 

- 5: inverse of acceleration due to gravity (GO) 

= 6 : feet to meters 

= 7: meter to feet 

= 8; nautical miles to kilometers 

== 9: kilometers to nautical miles 

CMP* Pitch damping coefficient 

CONST: =34,163194 x 10“^ ®K/meter, Constant in atmosphere 

pres sure f ormul a . 

DENS: Current atmospheric density, 

DENSO: Air density at sea level = 2.3769 x 10”^ slug/ft^, 

DIR* Direction table for wind. 

DMAX: Maximum step size during integration procedure. 

DMIN: Minimum step size during integration procedure. 

DPS; Array for storing integration variables S(14). 
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DT: 

DTAI: 

DTLA: 

EL: 

ELGD: 

EPBIG: 

EPTINY; 

ERW: 


ESQ: 

ESQl; 

ESQI: 

F; 


GAMA: 

GAMAL: 

GDLAT: 

GDLATL: 

GM: 

GR: 

GO: 

HALTS: 

HLIM: 

HPI: 

HPREV: 


Current step size for integration. 

Initial step size for integration from apogee to impact, 
initial step size for integration from launch to apogee. 
Current elevation angle (geocentric) . 

Current elevation angle (geodetic) . 

Tolerance used in step size control (if current error 
> EPBIG, step size is decreased) . 

Tolerance used in step size control (if current error 
< EPTINY, step size .is increased) . 

Vector of error weights used in predictor-corrector 
method. These weights, in conjunction with EPBIG and 
EPTINY, are used to control step size during integrations. 

_3 

Square of the earth's eccentricity =6.6934217 x 10 

1 - ESQ 

1/ESQl 

Array to store the functions f in the differential 
equation y' = f. In F(I,K), I denotes: 

I = 1,2,3: inertial velocity components along x,y,z. 

= 4,5,6: inertial acceleration components along x,y,z. 

= 7,8: rate of change of angular momentum about y 

and z. . 

= 9,10,11,12,13,14; time derivative of the body-to- 

inertial rotation matrix. 

Current elevation angle (geocentric) . 

Launch elevation angle (geocentric) . 

Current latitude (geodetic) . 

Latitude at launch (geodetic) . 

Product of the universal gravitational constant and the 
mass of the earth. 

Acceleration due to gravity at the current altitude . 

Acceleration due to gravity at the sea level. 

Geopotential altitude table for the temperature curve 
from 1962 U.S. Standard Atmosphere. 

Altitude from the surface at which the trajectory (impact) 
has to be terminated. 

Half of 7T, 

Previous interpolation value used in atmosphere table. 
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IPAZ: 

IPR: 

IPRINT* 


IPTj 
I PUNCH ; 


IROT: 


ITAPBI: 

ITAPE2: 
I TAPE 3: 


ITAPE4: 

lYY; 

JS; 

JSPENT: 


KPREV: 

KSP: 

KSTEP: 

LAT* 

LATL; 

LCG: 

LONG: 

LONGL; 


Surface azimuth below the impact point. 

Surface range below the impact point. ' 

Option for printout of the trajectory variables: 

- 0 : summary only 

= 1; detail trajectory with summary 
= 2* debug option 

Current time at the impact. 

Option for punchout of range, burnout, flight elevation, 
unit-wind effects, and delta f(z) curves. These are used 
as input to 5D wind weighting. 

= 0 * no punchout 
^ 0 • will be punched out 

Option^ to consider rotating or non-rotating earth: ‘ 

= 0 : non-rotating earth 

^ Oi rotating earth 

Disc number (=5 always) for inputting on cards the 
NAMELIST data. 

Disc number for inputting thrust and mach tables . 

Disc number on which processed vehicle tables are stored. 

The tables are cross-interpolated prior to the trajectory 
simulation phase for faster program execution, (The two 
tables with time as an independent variable are given a 
common time axis; the four tables whose independent variable 
is mach number are treated similarly.) 

Disc number to store data for spent stage calculations. 

Pitch (= yaw) moment of inertia. 

Subscript for separation times. 

Option for calculating trajectory for spent stage (s): 

= 0 : no calculation 

^ 0: will be calculated 

Previous subscript used in atmosphere interpolation sub- 
routine . 

Number of entries in the spent stage table for drag co- 
efficient , 

Number of integration steps required before doubling the 
step size. 

Current latitude (geocentric) . 

Latitude at launch (geocentric) , 

Distance from nose to the c.g. 

Current longitude . 

Longitude at launch. 
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LPREV; 


MACH: 

MASS; 

MOOT: 

MLAST; 

MPL: 

MPREV; 

MPROP : 

MSEP: 

NAME: 

NANG: 

NCON; 


NGOOD: 


NLEV: 

NMj 

NP; 

NPL: 

NPl; 

NPMAX: 

NPST; 

NSEP; 

NSMAX; 

NSP: 

NVAR; 

OMEGA; 

ONE: 

ORDER; 


Array for storing previous subscript used in interpola 
tion of thrust and drag tables . 

Current mach number. 

Current mass of the vehicle. 

Current rate of change of the propellant mass. 

Number of levels in wind tables. 

Mass of the payload. 

Previous subscript used in the wind tables. 

Current mass of the propellant. 

Table of separation masses , 

Name of the vehicle. 

Number of entries in array ANG (launch elevation) . 

Option for initialization for various purposes : 

< 0; Initialize all variables if launch coordinates 
(latitude, longitude) are changed, 

=0: Initialize all variables if launch settings 

(elevation, azimuth) are changed. 

= Of Wind not present. 

>0: Winds present. 

Number of steps integrated at the current time without 
changing step size. 

NT^er of entries in the wind altitude table. 

Array for storing number of entries in array A. 

Current phase number. 

Number of payloads . 

NP + 1 

Maximum number of phases that the program can take. 
Total number of phases . 

Number of weights to be separated. 

Maximum number of separations that the program can take 
Total nxomber of separations. 

Maximum number of entries for array A. 

Earth’s rotation rate (rad/sec) . 

1.0 

Array to store the degrees of polynomial for each A in 
cross -interpolating the time and Mach number tables. 
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PHI; 

PHIL; 

PHT; 

PI: 

PLM; 

PSI; 

PSIL; 

PRES; 

PRBSB; 

PRESO; 

QPI: 

QS: 

R: 

RE; 

RKSTEP ; 


RMASS; 

RQ: 

RTI; 

RXY; 

S ; 


SLOPE : 
SP.; 

SPEED; 

SPM; 


Polar longitude (0° <PHI <360°) , 

Polar longitude at launch (") . 

Array to store times at which a new phase begins » 

Constant it. 

Array to store payload masses . 

Current azimuth (0° <PSI <360°) , 

Azimuth at launch (*') , 

Current atmospheric pressure. 

Ratio of air pressure at a given altitude to. sea level 
pressure (Ref. to 1962 D. S, Standard Tables.) 

Atmospheric pressure at sea level. 

tt/4.0 

Product of dynamic pressure and aerodynamic reference area. 
Current distance from the center of the earth (inertial) , 
Earth's radius at launch site. 

= 1; Integrations performed one step by Runge-Kutta 

method. Used near the phase termination boundary. 

3; Integrations performed three steps by Runge-Kutta 
method. Three steps of Runge-Kutta are used in 
obtaining starting values for predictor-corrector. 

Total mass of the vehicle at the ignition minus payload. 

Equatorial radius of the earth. 

Array to store rotation matrix (from topocentric to 
inertial) , 

Projection of R on the inertial X-Y plane. 

Array to store the integration variables: S(J) 

J = 1,2,3; inertial position coordinates (x,y,z) 

= 4,5,6; inertial velocity component (along x,y,z) 

= 7 ; pitch angular momenta 

= 8; yaw angular momenta 

= 9,10,11,12,13,14; the elements of body-to-inertial 

rotation matrix at the current time , 

Change in atmospheric temperature with altitude. 

Array to store the integration variable (s) for spent stages. 

Speed table for wind. (Wind direction is denoted by DIR, 
and altitude by aLTW) . 

Mass of the spent stage. 
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SOUND: 
SOUNDO ; 
SRAZ: 
SUfiM: 
SY: 

SZ: 

T: 

TAV: 

TBO: 

TEMP: 

TEMPB; 

TEMPO: 

TEPS: 

THETA; 

THETAL: 

THR: . 

TNPl: 

TPI: 

TSEP: 

TSLOPE : 

TSP: 

T3D; 

VA: 

VEL: 

VELE* 

VELN; 

VLIM: 

VPREV; 

VT: 

VOLIM: 
WE : 
WIND: 


Current velocity of sound. 

Velocity of sound at sea level. 

Current surface azimuth (bearing) . 

Rocket mass minus propellant mass of a thrusting stage. 
Angular rate about body y axis . 

Angular rate about body z axis . 

Current time. 

Array to store initial time, altitude and velocity. 

Time at the final burnout. Integration change from 5D 
to 3D here . 

Current temperature in atmosphere. 

Temperature points from the 1962 U.S, Standard Atmosphere 
Tables . 

Temperature (atmospheric) at sea level. 

Minimum time interval allowed for accuracy (not used in 
program) . 

Polar latitude (0° < THETA <180'’) 

Polar latitude at launch (”) , 

Thrust at vacuum. 

Time when the next phase begins , 

2tt 

Mass separation times , 

Slope of temperature with respect to altitude. 

Times at spent stage separations . 

= TBO. 

Components of air-relative velocity in body axes. 

Current relative velocity , 

Table of east components of wind velocity. 

Table of north components of wind velocity. 

If the current vEL < VLIM then apogee is assumed. 

Array for storing previous interpolation value in A, 
Components of velocity in topocentric system. 

Minimum velocity at launch. 

Current east component of the wind velocity. 

Magnitude of wind velocity. 
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WN: 
WNAZ ; 
WPLj 
XO : 
YO- 
ZERO; 
ZO : 


Current north component of the wind velocity. 

Direction of wind velocity measured clockwise, from North. 
Weight of the payload. 

Inertial X component at launch. 

Inertial Y component at launch. 

«-15‘ 

Practical value of zero = 1.0 x 10 , 

inertial Z component at launch. 
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1 p p p 1 p 1 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 1 1 1 1 1 1 1 1 1 1 1 J 1 1 1 

p p 1 p 1 p 1 

i 1 1 1 1 > 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

I.J. J 1 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 < 1 1 1 1 1 1 1 1 1 1 1 

1 1 1 1 1 1 1 

Fi -( Ni Oi Rj Ml Ai If |P|0|R|C|E| j C, 0| E| F| F| S| i i i < i i i i i i i ■ i i i | 

p p. 1 p p p .1 

r 1 1 1 1 1 1 .1 1 1 1 l-.i 1 ! 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . 

P ! 1 

1 1 1 1 1 1 t t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 

P 1 .1 1 p p 1 . 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .1 1 1 1 I 1 1 1 1 

1 1 1. 1 1 Pi 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 

P 1 1 P p p P 

1 1 1 1 1 1 1 1 ! 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 1 1 1 p p p 

1 1 1...1 1 1 L. 1 1 1 ! .1 r 1 1 1 1 1 .1 1 j_ 1 1 1 i_ 1 \ .A i_i 1 1 i_l 1 1 1 1 

1 P P 1 P P P 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


R|0,M, |N|0|S|E|), ,D|IiS|T|A|N|C|E|S| i i i i i i ■ , i i i i i i i i i i i 


1 1 1 1 1 1 1. L 1 1 1 1 1 1 1 1 t 1 1 l_L 1 .1 I 1 1 1 t 1 1 1 1 1 1 1 1 1 1 

BBBIHflfl' 

1 1 1 1 1 1 1 1 1 r 1 1 1 1 1 1 1 [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I ! 1 1 

MHBRHHi 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 f 1 1 

p 1 1 . 1 p p 1 

II I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 1 1 1 F 1 1 ( 1 1 1 1 

1 1 1 1 p i_i_ 

1 [ I...I 1 1 j. 1 1 1 1 . 1 1 1 1 1 .1 1 1 1 .1 1 1 1 1 1 1 1 1. 1 1 1 1 I 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 

imiiH 

1 1 1 1 1 1 1 1 P 1 1 1 1 1 l_i 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1. 1 1 1 1 1. 

I I 1 1 1 !._.l 

1 1 1 .1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 P 1 

1 1 1 1 1 1 1 ' 

t 1 .1. 1 1 1 1 1 1 1 1 1 . 1 1 1 1 t 1 1 1 1 1 1 1 1 1 P 1 1 1 1 1 1 1 1 1 1 1 

1 1 1 1 1 1 1 

1 1 1 I 1 1 I..1. 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 1 p P 1 1 1 1 1 1 1 1 1 1 

1 1 1 1 1 1 1 

-J.J 1 .L_1..L i_L..L J J. 1 1. 1_J 1 p 1 1 1 1 1 1 1 I 1 P _1 1 -P 1 .1 _.l 1 1 1 1 P 

BfllBflHfl 


IDCNTI nCAT50*r 


>733 34 >6} 3IS9 60 e 1 *1 A 4 &> A6 jft * fit 4^>0 71 77 7l 


NASA Ml.gas 









NASA MI*P09 ll<44> 








NASA IVI n<«0) 
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NASA WI.963 t1«OB) 
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NASA Wl*009 
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NASA Wl.fldS 
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FORTRAN 
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NASA Wi*0BS 



fortran 
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APPENDIX C 
SAMPLE OUTPUT 



SENS-5D Ca.CULATIONS BE5IN »= 


[NPUT DATA S ET NO. 1 . . . 

5ULI?r AiE=lO=^0, IPUNCH=1 $ „ 

EBLI3T MPl=1»pL'1»50» AZGDI, = 12D, NANG=11,ANSb 70» 74, 7Ai76,78i§0i 82,84., 
86.88,90, JSPENt«i i ... 

EPLI^T WPL»50, AZGDLe90, ELC0l.=80ji WInD»20»WNAZ»80B, 

NLgislSj At PVbQ. 5 0.100.150 .•200725 0. 80D.i5ub,.5pQOriDO0O,3(iP0O, 

60000 ^ 

EOLjgi WPL*50,AZ(?DL*90> NANGbii,ANG»70, 72, 74,76,78; 85 , 82, 84, 86,38,90, 
WIMOs20»AH1GH=60000 J ' . ’ 


P O 

hg ^ 

10 

§ 

'M 


-C-1 



I/^PJr DATA SET NO, 2 
NlKf (JAJtlN 

ROCKET yEIGHT, HtJRNcOUT TIME 
1^19,0, 22', -0 

N, STJRT TlM3§ FOR PHASER 

4, 50,017, 3,52i,:-l7.,.Di,j2,g 

N, ?T4RT TIMES FOR -SPENT SlAfiSS 
1» ^.?23 Vi 

N, TAiLf Of ^SIGHTSVFOR. spent STAGES - 
li. |68,0 . I'. 

PHASE REGINS.', comments THIS 1 5 THRUSTING PHASE 

NOZiLg PKITarSa^ length from nose to NOZZLEi NEF, ARESt REF, DIAMETER 

i,S. 2?i5, 1.474, 

N, TIME Tabus for thrust and mropellant height 
22 , g.Of a'.Oli .0,04,^ f.gsi'g.g?, O'.ll,. u,|4t l.lAt 4,14, 2,p4,; 2,34 
?.99, 2,64, 2.79, J^^99, 3.Q9, 3,21. 3,29, 3j3*| 3j40.. 3,44, 3t54 

N, THRUSTS at __ "I^SEA UEVgL ^ 

22, 1933;o, 26439';04' 39417 , q* 41628,0, 42397, Oi 4S569lpi 4326J.0 
43531, 0, 44705, 0, 45667 ;b, ^ae^oio. 46723io»"45575,o, “429'74',0 
i7975 ioi 33648,0, 24o34;8,' 16343^0,- iB8S4Jii‘ O',0 

N, PRqPE^LANT HEIGHTS (NOT HAS?J 

22, 73B,o, 737,26, 732, lOi 7|9,98, 72i;29* 7Q7l37tlS5|,07j 481^01 

54§j63i_27^i7ai 203; ?7, 146'.6: I30.42t 9s;727 SJTil. 34:69, i6i60 
f,?2, 4,95i 2,22, 0,68, O.p . _ , 

N, TAijLg Of Mach nos, for axuC drag coeffs, 

11 , g,p, o;?5, i',Qt 1,20, i.flpi 2. o', 2, 4, 2 , a, 3,8,..4,S,' 5,.0 
N, ^XIAL drag ccjeffs, 

U( ?,475. 0.595. ,925, ,874 -,78, In,'. ,661 i Sifl. j I? 455 .. . . 

N, TIME TiBUi for. C,G, (FRO M NOSE) DISTANCE AND PITBH HOHEWT OF INERTIA 
2i 0.1 3,54 

N, C,G, distance from Ngse 
2. 17,383, 16;0 

N, P.ITCM M“MSNTS of inertia , . . . • , . 

2, 1&3P., 1255. . . ‘ . 

N, TaSLL QF,.jJACH. NPS, FOR RIT CW DAMPI N G COEFFS.. . 

12, P,t .9, 1 , 0 , 1,5, 2,(01 215i 3i0i 3,5, 4,Q,j. 4,», 5,g, 20,p 

N, PITCH damping coeffs, . . 

12, 15,89, 19.95,. 19, '72. 21,??, 14;66, 11,35,, 9,?, B.W. 7;.33,_6,98. 

N, TAsig OF MACH NOS. FOR 3Lfl?i-Dr-N0RMAL laBCg CObEfS;;_. 

12..B. B. o;9...1.0, 1 ,5. g,DI 8,5. 3,P. 3.5, 4,0....4l5, S.q, 2d:q 

N, SU0PE OF formal force CQEFE5,._. 

12. 0.27, /,'07. ,325. ,24. . 20, _;i7f....l55t.xll,,.4l5i 4.125, 

N, TABUS OF Mach NQS, for Q.F-* CFROM NOSSl ClSlANetl. . 

• ” ■ 




9,t D,?| l^Qf 3;0i 3 .§j_?.Qi, 1i?f s’,0,.20,.0. 

Mf §*•;. BISTANICES FROW N03B 


IZf jU|Zn»t?;5'33,l9,333il^.ifly.l8.75<17.9 t7« l7.5<17.D93,16. 25«l5 ,a33 

i5,833 _ „ _ ... 

PHAS K NO. 2 BEGIMS . COM MENTS THIS IS COA8T1nS PHASE. 

MOZZLe e*I T AREAi UENlJriTPROM NOSE TO NOZZUEi REF, ARtJ}, REF, DIAMETEf? 
0,y,?,9l0!230,1,54 _ 

M, TI»16 TA0I.B FOR THRUST AND l»BgPEI,l,ANT WEIGHT 

0, , . . 

N, THRUSTS AT VACUUM 


. La.Ax 

N, PRQPpLLANT WEIGHTS (NOT MAS5) 

i< t>! , . ■ ’ 

N, Table of hash nos, for axial drag coeffs, 

7i 2;0» 3,0, 4,0, s,0»'6,0, BjO 

N, AXIAL OHAS COEFFS, ' . . ‘ 

»^2« t^ X«i 1 42 > i36 

N,'TiHE TABLS'F.TR C,G.'(FRpM NOSG) DISTANCG'ANO PITCH MOMENT OF INERTIA 

0, ■ ‘ ■ ■ 

N, G,q, DISTANCE from NOSE 
i, 7,783 

N, PITCH MOMENTS OF (NtRTlA 

It L29*3__ . . . _ 

N, TA8Lp OF HACH nos, FOR P}T5M DAMPING COEFFS, 

9i Z',0( 3.0, 4.0, S(Ot 8,0, 7,0. Bjo^'zO.O 

N, PITCIjl 0AHRIN5 COEFFS, 

9< I7j,4. 157',, 122,1. 104, 6t 88,94, 7B,4Bi 69,74, 69‘,76, 69’, 76 

N, TAgiE OF Mach nos, for slope-of-normal forqI cotFFSj 

’l, i.i§€ Z’lOi.iiPe 4,0, 5,0* 6.1),' 7'.D, 8'.0,'20‘.5'~‘ " 

N, SLQP| of normal force coeffs, 

B|58, ,47, ,33, ,26, t22,..,19, .17, ,16, il6 . 

N, tabls of Mach nos, for c.r-, {from nosej disianchs 

9, It?, Z^O, i|0, 4,0, 5,0* 6|_0, 7,0, B,0* 20,0 

N, C.P, distances from nose 

9i U*g4Z*lal833.1Q,3 j3,.9.9l7,9.5.9.167,a'.875.B;o25,6.625 

PRASE NQ,3 begins, comments THIS IS THMgSIING PHfiSE 

nozzle fXir area* LENOTH from NQSE to nozzle* REF,'AREfli ref; niiMJTER 

a,g2<?,15.y83,0;230,0|&42 
N, TIME TA8L3 FOE THRUST AND PROPELLANT WEIGHT 
21i 0.1 Oi 04, .08, ,15, ,5, ,7, 1,, 1.2, i,4, i,6, 2,^ 2 , n. 2;dt .2 , V8 

. 4*84^.. 3^ti_3^26i.. 3. 4A, 3,3. 5.. 6 .0 

N, thrusts at VACUUM 

21, 2321., 7871,, 767li, 7546., 7671 , a 8i71,, agji,, 87A6, . _ 

9871,1 8971,, 994J,, 9521,, 7371,, 3171,, i57i;j J9_6, 5Z1, 



0,0 

N, PRQPELUANT WpIGHTS (NOT MASS) 

21t 11?|. 118, rti tl7,9:;,» 115i||. 105173, 97,0, 86;j6l 19, 03, .71, 34 

^3,48, 47.54, 31,30* 14,gfl* 6, ,43/ .4,17* 1,46, ,'61* ,24* .j;, 0,,0, 

N, lAgLE QP MACH NOS., FOR A XIAL DRAB COEFFS. 

8, l,j, ^,0, 3,d, 4,0, 5,0* «,0, 8,p, 20,0 
N, fiXlAt, DRAG COEFFS, . 

8, .d',7B, ,62, ,53, ,46* f42* ,38| ,34, ,34 

N, TIME TArtUB FOR C,G, (FROM’ NOSE) DISTaJCE AND PITCH MOMENT OF INERTU 
2i 0,9., S'.'O 

Ni.S.a, BJSlM£!:.F..Rg« NOSE . ; 

2, 7,Z83, 7,15? „ . 

N, R]Tcy M0M3NTS OF INERTIA 
Zt 3|29,?*. 56,3 

N, TASLi OF Mach nos, for pitch camping CQEFFS*. 

9, 1,5, 2 . 0 , 5.0, 4.0, 5,0, $,0, 7.0, 8,g* gC.O 

Nj_aiIC!i QAtlflNQ .QOgFFS. 

9, 174,4, 157,, 122.1, 104, 6; 88,94, 78,48, 69i76| 6?;76j 69,76 

N, Table of mach nos, 'for slope-of-noSmal fo«ce cocffS, 

9, 1,1, 2;g, .3,0, 4.0, 5,0, 6,0, 7,0, 8,0* 29*0 

N, SLOPE OF -normal force COEFPS, 

9, 0',?a, ,47J ,33, ,26 )t ,22* ,19, ,17, ,],6[ ,16 

n; Table qf MA-y ^s, for c.p- , (from nose) distancbs 

9, 1,5'J H7'0, 3,0, 4,0, 5,0, 6,0, 7,0, 8.5, 20,0 

N, C,P, DISTANCES FROM NOSE 

9, U ,0 42, If), 833, 10.333, 9 , 91 7, 9, 5, 9, 167, 8 , 875*8 825, |;625 
PHASE NO, 4 begins, COMMENTS-, THIS IS COASTING A'ND UgST PHASE 
N074UE EXIT area, length FROM NOSE fo NOZZLE^ REFj ARSfej REF, OIAMETER 

aiai0j9i.Di3aDi.C^Q- . 

N, TIME table FO^I THRUST aND PROPELLANT HEIGHT 

I, 0, 

N, Thrusts at vacuum 

1 , 0 , 

N, PRQPfc'wLANT HEIGHTS (NOT MA9S> 

, ... . - 

N, tabl^ of Mach nos, for axial drag coeffs, 

8, 1.9, 2,‘0, 3,0, 4,0, 5,01 6.0, 8,9, 20,9 
N, axial DRA3 coeffs, 

9, 9.?3, '.72, .,59, .51* ,45, ,42, .36, .36 
N, TIME table for C.G, (FROM NOSE} DISTANCE ANO PHCH MSMENT OF INERTIA 

li-Q, 

N, G.G, distance from NQSP 
1, 0, . . 

N, P.ITCH MOHENTS of inertia 



li.Q. 

M, TAflLf OF MACH NOS, FOR PJTCH DAM PING CoSFFS. 

li Qt 

N, PITCI;< OAH^IMG COEFfS, . 

N^'^f aIle of ‘<ach tTris” FoF'sT^FsbF^NORMAjL FpRc£ 'coFtf's',' 
is 0 1 ... , 

N, SuqP|: OF 'I0=?MAL FORCE COEFFS. 

li 0, 

N, TA^LE of Mach nos, for C.P-i (from, N0S.EJ OlSTANCbS 

. ’ II ^ 

N« 5,P, distances from nosb „ 

is Of . 

SPEtiT SIAGfc N0't5 BEGtNS, COMMENTS 

ijMI OF separation, SBP'araTED HEIGHT, REF. AREA 

36AVg, Qip 

N. Tabl e of mac h nos , for axial drag coeffs. 

3i ?,t 2«.0 

N, AXIAL DRAG CIJEFJS, 

2, P.i. 0, 


NASA WAuLops Flight uemter 

WALLOP? island, VIRGINIA 

TRAJeCTDRY SU'INARY AT B09N«0UT, APOGEE AND iHPACT 
VEHIGLg = NIKE CAJUN ’ 

pat wOfiQ i sn'.oo UBS ■■■ 

t,AyNQH AZ = IJiOjoO DEG 

Wien = ZERO 

EARTH = rotating HOOPu 


apogee 


alt Range 

(FT? (NMI 


EL 

THU 

(□EG) 

CSEPJ 

70.00 

72.00 

74.00 
76-jOO 

78.00 

-177;94 

182,59 

106,62 

190,39 

193;37 

00,0 0 

196lo5 

82,00 

64.00 

66.00 

f9fl;l7 
199176 
2f>0 ,35 

0 0; 00 
9O;(J0 , 

201,46 

'2Ci;54 


42R921' 63j'9 

- 559973 ; 59ia 
474109, 54i9 
493497, 4°j4 
519755'. 43 ;? 
525976; __ 36,6 
537700, ■ ■29,4 
546997 j 2?;o 
553402, 14,3 
557300, 6^4 
55R1H5; 1.6 


burn out' impact 

tlHE ‘ ■ AUt ‘RaNSE 9EI FLT/EL FLT/aZ TIHg PaNCE AZ 

(SRC) (Ff> (NHi (FTZSeC) ' (DEG» (DEG> ISEC) (nEG> 


22’, DO 

22 1 o' o’ 

22i0U 
22;00 
22;ou 
22 >00 

55 920', 

569B5; 

57930, 

58705; ■ 

59532', 

60165'. 

4,1 

377 

375 

217 

ll? 

2*1 

-- 577F773T- 
5765ld0 
5776,98 
5789,34 
5799,01 
6000,36 

63, 3o 
6S;96 
60;63 
7i;3i 
74;t)U 
76, 7d 

lir.94 

119,91 

119,87 

119|82 

119,76 

119,68 

359|64 
367; 39 
374 36 
300,39 
385|44 

- 1207 s 

120,0 

iio;2 

99;i 

06;e 

73,4 

121125 

121,33 

121j42 

12l|54 

I2i;69 

22,00 

6 O 684 ; 

61087 

61372; 

61539; 

61509', 


SB16.T5' 

5020,49 

“"797T0 

02;ii 

lt^5 

119,35 

309TJ2 

392:61 

WT2 

121,9 0“ 

22i00 

173 

44 ; 2 

12?i24 

22i0ii 

OiS 

5824i05 

84;83 

87,53 

09|73 

110,94 

394 ; 70 
395179 

20^ 

122195 

22,00 

0,4 

5026;i5 

117,62 

i2;9 

125-,35 

22,00 

OrO 

5026,76 

323,35 

395 '.90 

3.2 

203i?2 







' 


■ 
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NASA W6LL3PS FLIGHT CENTER 
WALLpPS ISLANDi VIRGINIA 

TRAJECTORT SUMMARY FOR SofiNT STAGES, 

VEHiqL| r NIKE CAJLN _ _ 

PAy LOAD 's ■5r.00"LeS fPARENT VEHICLfT ■ ■ '■ ■” “■ 

launch A2 s 120,00 L'EG <;,> 

WIND = ZERO 

EARTH = rotating MODEL 

separati_q,n impact 




LAUNCH 

NO,' 

TIME 

weight 

ALT 

^RE 

AZ 

VEL 

FLT/EL 

FTT/A7" 

Time 

■*TTATn5E~ 

a7" 

EKues) 


(SEC 1 ^._, 

_ ISLUC> 

(FT) 

(NMI 

<pbG) 

(FT/SEcI 

(DEU> 

tUEG) 

<SEC) 

(NH) 

(DFGJ 

70.00 

1 

3.52 

564 05 

5«-<o6 

5762, 

5842; 

59ii; 

0,37 

b|3? 

119;63 
119 S7 
ii9:Si 
11^,4^ ■ 

3308,44 

33oQ;3o 

68,29 

119;62 

199-.64 
202 ; 32 

39.27 

120,59 

72.00 

1 

5:?2 

7D,46 

ii9;58 

35 ‘.93 

120 i 61 

7A,cg 

1 

3,52 

3,52 

564^00 

o;3p 

TSoBJnB' 

72,64 

ii9;s2 

ii9;44 

204’.7g 
206 . 5^ 

32;4o 
287 W 

l2o;64 

76,00 

1 

563, PC 

5974, 

a.zi 

38o7i?6 

74,82 

120.66 

78.00 

1 

3,52 

564 , 00 

6030, 

b.zS 

il9l31 
119 14 

116.90 
118:62 
li? 76 

3357783 

77’, 00 

119;34 

ii9',i9 

2 o 8 ;s? 

24;81' 
20 . 8 q 

120.70 

40.00 

1 

3.52 

564,00 

6076, 

b;i9 

3357.72 

7?;i8 

PlOjZi 

2ii;43 

1Z0.74 

4?. op 

1 

3.52 

563,00 

6115, 

s.'i5 

33o7;62 

Hi57 

iie;98 

16 ’.67 

120 ‘.80 

44,00 

1 

3,52 

564,00 

6145, 

p.u 

33o7:S6 

63,55 

118,62 

ii7;9o 

2iz;40 

12 '.45 

120;9n 

46. op 

1 

3,52 

564.00 

6164. 

dio? 

33o7T48 

55,74 

213 oS 

fl'.l6 

121 ;o9 

48, op 

1 

3.52 

S63.no 

6177, 

OfOP 

Its SB 

3307, 4S 

87,92 


213.^0 

3 ; 03 

121:73 

Vo, op 

1 

3.52 

564,00 

618p , 

O,oti 

355,31 

3307i44 

89,80 

353,52 

?i3;4g 

0'.53 

291'.93 
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JilABA WAuLOPS rulGHT CENTER 
HAU9P§ ISLANP* VIRGINIA 

r(Z)-C\JRve AND DALIISTIC HIND FACTORS 

VEMl CUe s NIKE 

PAy IjOfie 5 50,00 LBS 

UA^/N^H SL = 80,00 LEG 

LAyNgH AZ s 9C,6 o LEG 

WltJD = 20,00 FT/SFC ( 6 , to M/SEC) 30 0,00 2 Eg AZ FROM NORTH 

gAgTy t NON-ROTATINC MOOgU 


WIND alt 


impact 


F(Z) 


- 

<FT) 

tM) 

(NH> 

IKM) 


i 


?. 


140.18 

0., 

8 

i 

i 

i 

6 

7 

5o; 

loo; 

ISO 

2od; 

850, 



30, 

46, 

61. 

76, 

n± 

73,90 
72,35 
71 ,21 
70, 3i 
69,72 
69,17 

136,87 
134 , do 
i3i^e7 
13 q, 2Z 
i89il3 
l2B.il 

0,13472 
0,25345 
0,33946 
0;4082d 
6, 4 5358 
n .49609 

9 

9 

1? 

u 

18 

luoo; 

300 0 ; 
lodee; 
30O0O', 
60000 . 

305, 

914, 

30Ae; 

9144; 

182B8, 

66,07 

64,23 

63,90 

64,47 

62,88 

182,37 

li8i95 

iia,34 

li9,40 

116,46 

0,74ql8 
. 0.86880 
0,91596 
0;87 i 96 
lioodoo 


DF<Z» 

FtZj^FiZsl) 


0,13472 
0,1187^ 
d ,5f600 
di 06674' 
0.6453^' 


WIND ALT 
<FT) <H) 


C.U U A 
2?0, 


0 . 


15, 
3o } 
46« 

76 , 


Q ,24408 

^ f , 

1000 , 

" 305 

0,14863 

30 00 , 

914 

0,02716 

iuodo, 

5 o 4 R 

d ,o 48 oi 

3 ddoo, 

^9144 

0 ,18004 

60000 , 

’10286 


C- % 





NAfA WALLOPS FLIGHT CENTER 
WALLQP? ISLAND, VIRGINIA 

UNlTjWMD E^TECTS, CCRIOLtS DEFLEXION ~AfTP ^ANG'£ DERIVATIVE 

VtHI£Le = HIKE C AJUN 

pAy = 5C.00 UBS 

launch AZ = 90 ,00 'DEG 

WIND * 20,l0, LIV^O *ilC M/SEC}, 

EA^Ty : M3N»R0 T a'T- 1N 6 MQDgL 


F=IJ=S ?YSIEH 

• 

- 







LAUNCH SL 
IDEG) 

NbsWiND-- ’ 

RANGE«NM)‘-j 

impact 

A.Z1DER). 

CORIOLIS 

NORTH(NH) 

DEFLECTION 

'■eaST(nm> 

RANGE DERlV 
<NM/Dp> 

UNIT 

bEAU(NM/FRS> 

. WIND 

tAlL<NM/FpS} 

effects 

cR05SINrt/FP9) 

7l)iCtl 

128.69 

’^>49 

c2i65 

-e;58 

i3;i393i 

0,30899 

•0,40667 

0174411 

72, Op 
74, op 
76, Op 
78. op 
1 ^ 0 , op 
§2,gg 

120,89 
lllISfc 
100. 8i 
. 88.77 

75.69 

61.70 

:-'8?!’43 

*69.35 

B9',75 

89,12 

66,94 

68.68 

92,60 

o2,52 

q2,39 

62,23 

=2|04 

5ii82 

61,06 
61,35 
H;7B 
=?l03 
5? 33 
=3 59 

54,21265 
64,94776 
=5:45032 
56:34817 
s6, 72390 
16007' 

0,4o376 

0j494ii 

0,57050 

d,65901 

6,72354 

6.78922 

■50, 50130 
-0(98969 
od;65790 
=0 ,72564 
50.78550 
-0183298 

oi8014o 

0102962 

01846Q3 

0j86415 

0,80042 

94.09 

96.09 

pe,cp 

»o,og 

46.95 
31j66 
15, 9B 
1.4B 

88.23 
8?, 34 
84,66 
368.00 

:ll8S 
6l ,31 
=iib4 

60,75 

?3}§a 

53199 

E^;pa 

?5.ii 

57,49689 
£7170074 
£7 87882 '’ 
57,96376 

0,83lff6 
0 |86653 
6,89329 
6,90417 

-0.86762 

30j89372 

60,90552 

=0,90612 

0,89212 

D jSPBfn 
0,9 d347 
0 90409 


M=K=S system 









LAUN6H =L 

NOcWlND 

impact 

CORlOUis 

deflection 

RANGE DGrTV 

UNIT 

= hind 

EFFECTS 

(D68) 

RANGECNM) 

AZIDEG) 

NORThI'KM) 

EASTfKMI 

CKM7DEG) 

HEAD(KM/HPS} 

TAILCKM/MPS) 

CR0SS(KM/MPS) 

7g,gp 

238,34 

.8£i49 

54.91 

Si; 0 7 

=5.90660 

1x87721 

_=2j.47lCii _ 

4'.52131 

72.00 

223,89 

89,43 

= 4,52 

si, 96 

57<S94S^ . 

2,45326 

=3,04597 

4 7IJ424 

74, op 

206,61 

89.35 

54,66 

=3149 

=9,16326 

3,00228 

=3,55509 

4,86942. 

76, Cp 

166,73 

89.25 

s4|43 

5^129 

- 40,59399 

3,51549 

=3,99749 

5,01656 

78, op 

164 ,40 

S?!lZ 

w 4 8 X ^ 

i§|77 

511,75680 

4,00422 

=4,40906 

5.14058 

§o,op 

140,18 

88,94 

-3,78 


»12i45266 

4,39634 

=4177282 

5125069 

92iOP 

. 114,26 

88,68 

93i37 

-- 0.79 

513.26044 

4.77107 

=5,05886 

5.34933 

^94,0? 

86,95 

88,23 

=2,92 

57118 

513,88485 

5,05448 

=5,27177 

5,42060 


58,64 

87,34 

=2,43 

si', 53 

=14,26178 

5,26513 

=5,43035 

5146182 

S 0B.OO 

29,59 

84.68 

51,92 

=7,70 

=14159157 

5,42772 

=5,50205 

5148956 



?o,oe 

THiS CgNCLUOES 


a. 75 36B»DD sl^a? ;g;76 3l4»749.?fl 


5,493fl2 


s5, 50568 


THE CALCULAT .t-ON. 



C-IO 


5,49823 



Tf»lPS FOR MTKF CA.tUtJ 

M. LfiUMCH FLPVflTTONS (DFG) FOR IMPACT RANGE, R-0 FLT-FL AND APOGEF ALT 

11 70.0(1 7?. no 74.00 7G.00 7R.00 00,00 02.00 04,00 06,00 

00.00 tiO.OO 

T»PACT ‘JArgfiFO IN NM 

120.1O 120.00 110.18 99.10 86.79 73.4? 89.17 44.21 28.73 

12.94 3.20 

9-0 FLT/eL IN HFG 

6-».30 68.96 68.63 71.31' 74.00 76.70 79.40 82.11 84. H? 

87. RR 89.73 

APOGFF ALTTTUnF TN FT 

429821. 482973. 474109. 493497. 510788. 528876. 837708. 

944997. 883482. 887308. •888125. 

N. AUTTTUnF TA8Ui FOR F(Z) AND DELTA F WJ CUPVFS 

1? 0 . 80 , 100, 150, 200. 280, 

300. moo. 3000. 10000. 30000. 60000. 

F(7) rilRVh 


0. 

(1.1347 0.2838 0.3398 0. 

4082 0. 

4536 

0.4961 

0,7402 

0 .8889 

0.9160 0.8720 1.0000 


- 



nPLTA 2(7) CliRVF 





0. 

n.1347 0.1187 0.0860 0. 

0687 0. 

0454 

0.0425 

(1.2441 

0 .I486 

0.0272 -0.0440 0.1280 





M, LAUNCH 

FLFVATT0M8 (DFG) FOR UNIT-KIND 

EFFECTS. 

RANGE 

DFRIV AND CORIOLI 

11 7(1. on 72,00 74.0(1 76.00 78.00 

80.00 82 

.00 84 

.00 86. 

0(1 

8 8.00 

90.00 





HFBI1 UNIT 

-><'IND FFFECT (NM/FT/SFCI 





n.30« 

n.404 0.494 0.579 0.659 

0.7R4 

0.785 

0.832 

(1.867 

II ,K93 

0,906 





TaII UNIT 

-fciwn pfFFCT (NM/KT/SfC> 





-0.407 

-0.8(11 -0.888 -0.-658 -0,726 

-0.786 

-0.833 

-0.868 

-U.894 

-0.906 

-(1.906 





CROSS IINIT-wIMD FFFFCT (NM/FT/SEC) 





0.744 

0.774 O.BOl 0.826 0.84b 

0.864 

n .880 

0.H92 

0.899 

0.90.-1 

0.908 





DA‘'(,r off 10 (NW/DFG) 





-3.189 

-4.213 -4.948 -5.450 -6,348 

-6.724 

-7.160 

-7.497 

-7.701 

-7,879 

-7.064 





COR TDLJS 

1,'FFI FCTIDN TO NORTH (NM) 





-2.682 

-2.604 -2.819 -2.393 -2.232 

-2.- 04 1 

-1.819 

-1.576 

-1.314 

-1.U38 

-0.792 





rORTOLTs 

11FP| FCl lOH TD fast (NM) 





-(1.877. 

-1.084 -1.347 -1.779 -2.033 

-2.330 

-2.588 

-2.799 

-2.985 

1 

> 

-3.112 








SeN5-5D CAICULATJCNS BEra-y 

INPUT (^AT* SET flO, I 
SDUST AZE’^OsEO $ 

JBUIST NPLsl.PLMsSC, A2'?UL=12n, NaNGsI , ANG«B0 , !PRINT»1 % 
SFLIST I 
$UtlST s 


c- iz 



NASft l<ALLCPS flight ffcNTtP 
wallops ISlauD, VlnMMlA 

THAJfcPTDRY Summary ai bupiiboijt, APocite and Jhpact 

VbtUGLfc = KiKt CAoLN 
Pay load s 50, do LQS 
LAUNCH A7 X 120.00 IjEG 
WIND * 7ER11 

EAUTri X rotating modfl 



OfcTAlUEO F'RIMT-OUT OF TRAJgCTQRV FOR UAyNCH EtEVATJOrj s 80.00- nSG 
UNITS ARE F-P-S-ntGRLE, FXCLBT RAMCECNM) AND ACC<U01. 


TJME. 

HA-NGE 

9C AGING 

AUIUIDE 

rUT/EL 

FU/ 8 Z 

^ELDCIT'C 

ACC 

LiA£bl 

THR 05 T 

-.ofiAJ. 

.J1/PB£S 

£IA_ 

QAUF.A. 

WE,LG„HT 


MINOE 

U <00 

0.00 

lao 

.00 

20 . 

79.66 

119.11 

0 . 

le.Q 

0.00 

29740. 

Q. 

Oi 

0,43 

0. 

4566.9 

0 . 

0. 

il.OO 

0.00 

180 

.00 

20. 

79,64 

119.11 

1. 

18.0 

0.00 

29807, 

0. 

Q> 

0.45 

Q. 

1566.9 

£. 

0. 

0 .00 

0.00 

laii 

«O0 

20. 

79.68 

119.10 

1. 

18,1 

0.00 

29075, 

0. 

Q. 

0,46 

0. 

1566.9 

Q. 

0. 

Ut 00 

o.co 

130 

.00 

20. 

79.62 

119,10 

1, 

18.1 

0,00 

29943, 

0. 

0, 

0,47 

Q. 

,1566.9 

0. 

0. 


4-.40 

180 

.03 

2a. 

79,61 

.119.. liJ. 

1,- 

-18.2- 

0,00- 

-3P010. 

o..— 

-ft 

-0,48 . 

...0*. 

1566.8 

IL- - 

D. 

o.oc 

3.0!) 

131) 

.00 

20. 

79.60 

119.10 

1. 

18.2 

0.00 

30078, 

0,. 

0. 

0,49 

0. 

1566.6 

Q, 

0. 

OlOO 

0.30 

180 

.00 

21. 

79.60 

119.10 

1. 

IB. 3 

0,00 

■30145, 

Q,.. 

. 0. 

0*49 

0. 

1566.0 

Q. 

0. 

U. 00 

0.30 

l8o 

. on 

21. 

79.59 

119,10 

1. 

ia.3 

0.00 

30213, 

0. 

D. 

0*50 

0. 

1566.0 

0. 

0, 

u. 00 

0.00 

218 

. 97 

21. 

79;S9 

119.10 

1. 

IS . 3 

0.00 

30281, 

0. 

0, 

0,5(1 

0. 

1566.7 

D. 

0. 

d -XD 

O.on 

13tl 

.00 

21. 

79 ; 58 

119.10 

2. 

18,7 

0,00 

30BS9, 

0. 

0. 

0.51 

0 . 

1566.5 

0. 

0. 

o-.-oa. 

O-.OO 

lao 

.410 

21. 

- 79;58 

. Ui-10_ 

... ._.o---ia-,8 — u.,00- 

-31024*— 

0*-- 

— 04 — 

- 0.32 - 

. 0- 1566,4 

.-0*. 

0. 

0 . 00 

0.00 

I8n 

. 00 

21. 

79;S8 

119.10 

2, 

IB. 9 

0.00 

31159, 

0. 

0. 

0.52 

0. 

1566.4 

Q. 

0. 

0.00 

0.00 

190 

.00 

20. 

79,58 

119.10 

3, 

19,0 

0.00 

31295, 

0. 

0^ 

0,52 

0 . 

1566.3 

Q. 

0. 

Q.OU 

0.00 

ISO 

. 00 

20. 

79 -.57 

119,10 

3. 

19.1 

0.00 

31430,. 

0. 

0. 

0*52 

0. 

1566. 3 

Q. 

0. 

OfOO 

0.00 

l 80 

.00 

20. 

79.57 

U9.10 

3. 

19.2 

0,00 

31563, 

0, 

0,. 

0,52 

0. 

1566.2 

0. 

0. 

0.01 

0.00 

190 

.00 

20. 

79 ; 57 

119.10 

3, 

19.3 

0,00. 

31700,. 

0, 

0. 

0*52 

0 . 

1566.2 

0. 

0. 


-0.00 

180 

.00 

30. 

79,58 

- 112,10- 

- - .3., 


-0,410- 

31 8 :^ 8 . 

- -O,- 

0-,- 

-0*52- 


-1566-,! 

- .0,.. 

0. 

UtOl 

0.00 

130 

. 00 

20. 

79,58 

119.10 

4. 

19,4 

0,00 

31970, 

0. . 

0. 

0*52 

0. 

1566.1 

0. 

0. 

J.Ol 

o.oo 

130 

.00 

21. 

79;38 

119.10 

5. 

20.2 

0,00 

33197, 

0. 

oi 

0*51 

0. 

1565.6 

0. 

0. 

OfOl 

0.00 

130 

.00 

21, 

79,58 

119,10 

6. 

20,4 

0.01 

33453 , 

0, 

0. 

0*51 

0. 

1565.5 

Q. 

0. 

o«.oi 

0.00 

130 

1 00 

21. 

79,58 

119.10 

6. 

20.6 

0,01 

33723, 

Q. 

0. 

0,31 

0 . 

1565.4 

Q. 

0. 

OtOl 

0.00 

130 

.00 

21. 

79)59 

119,10 

7. 

20.7 

0.01 

3399B, 

0, 

0. 

0,51 

0. 

1565.3 

0. 

0. 

o.-ai- 

0-00 

218 

.ij: 

21. 

29) 59. 

-U9.,in. 

7,. 

20.9- 

.0-Ull- 

■34269, ■■■ 

Q.,.. 

-Oa- 

-.0*50.- 

.0,.. 

15.65.1 


0. 

0.01 

0.00 

307 

,32 

21. 

79 ; 59 

119.10 

7. 

21.1 

0.01 

34539, 

0. 

Oi 

0,50 

0. 

1565. D 

0. 

0. 

J. 01 

0.00 

207 

.82 

21. 

79,59 

U9.10 

0, 

21.3 

0.01 

34809, 

0. 

0. 

0,50 

0 . 

1564.9 

0. 

0. 

0.01 

0.30 

207 

.92 

21. 

79.59 

119.10 

a. 

21.4 

0.01 

350SQ, 

0. 

Oi 

0,30 

0 . 

1564.8 

0, 

0. 

J. 01 

o.oo 

201 

.59 

21. 

79,60 

119,10 

9 , 

21,6 

0.01 

35350, 

0. 

0. 

0,30 

0. 

1564.7 

0. 

0. 

0 » 02 

o.co 

201 

i59 

21. 

79.61 

119.10 

13. 

23.2 

0,01 

37704, 

0, 

0. 

-0.40 

0 . 

1563.0 

0, 

Q. 

J + 02 

0.00 

201 

,59 

21. 

79,62 

119.10- 



23,5 

0 .01 

. 

n . 

n . 

n . 48 

0 1 

18A3.5 

f) » 

0 • 

0.02 

0.00 

201 

,59 

21. 

79.62 

119.10 

13. 

23.9 

0.01 

39865. 

0. 

0. 

■ 0*47 

Q . 

1563.3 

0 

0. 

U « 02 

•0.00 

201 

.59 

21. 

79.62 

119.10 

16. 

24,2 

0.01 

39406. 

0, 

0. 

0,47 

0 . 

1563.1 

Q. 

0. 

0.02 

o.co 

201 

,59 

21. 

79.63 

H9.10 

17. 

24.4 

0,01 

39689, 

0. 

0. 

0,46 

0 . 

1562. B 

a. 

0. 

0 f 0 3 

0.30 

201 

,59 

21. 

79,63 

119.10 

10, 

24,6 

0.02 

39965, 

0. 

0. 

0*46 

0 . 

1562.6 

D. 

0. 

0 . 03 

0.00 

201 

.59 

21. 

^63 

119,10 

19. 

24. B 

0.02 

40242, 

0. 

0, 

0,46 

0 . 

1562.3 

0, 

0. 

0.03 

0.00 

201 


21, 

79,64 

3.19,10- 

- ---20, 

2S-.0- 

0,02- 

40510, 

— Q 8 * , 

0-,- 

- . 0,45 . 


186g. f) 

Q-* 

. Q. 

U f j3 

0.00 

201 

.59 

21. 

79.64 

119.10 

21. 

29.1 

0,02 

40794, . 

. 0,. 

- 1. 

0.43 

0 . 

1561.8 

0. 

0. 

OtOA 

0.00 

180 

.00 

31, 

79,66 

119,10 

32, 

2S.8. 

0.03 

41822* 

.1, 

- .. 1. 

0,43 

1. 

1598,8 

Q. . 

0. 


c- \i\ 



T{Mfe 

flA'IfiE- 

9gW?lNr, 

altitude 

FLT/eL 

FlT/AZ 

i/ELOCITY 

ACC 

HAGM’ 

THRUST 

drag 

d/pres 

ETA 

QAIFA 

Weight 

WINdN 

O.O'l 

0.30 

IB.0.00 

21. 

79,67 

119.10 

33. 

25.9 

0.03 

41046, 

1. 

1. 1 

0.43 

1. 

1558.5 

0. 

o.as 

0.00 

ISO. 03 

21. 

79.67 

119.10 

34. 

25,9 

0.03 

41070* 

1. 

1. 

0.43 

1. 

1558.2 

0. 

0^05 

3.6ft. 

i8o.no 

.21*,.. 

. 79*67 

119.10 

35, 

.25.9 

-0*03L. 

,A1824* - 

— -1*-. 

- --1.- 

0*42 

1. 

1550.0 

0, 

0>0S 

0.30 

l8a,0Q 

21. 

79.67 

119.10 

36. 

25.9 

0.03 

41918. 

2. 

2. 

0.42 

1. 

1557.7 

0. 

O.uS 

0.30 

180,00 

21. 

'79,67 

119.10 

37. 

25.9 

0.03 

41942, 

2. 

2. 

0.42 

1. 

1557.4 

0. 

3.05 

0.30 

ISO, 00 

21. 

79,67 

119.10 

30. 

26.0 

0.03 

41966, 

2. 

2. 

0.42 

1. 

1557.1 

0, 

0.35 

0.00 

180.00 

21. 

79,67 

U9.10 

39. 

26.0 

0.04 

41990* 

2. 

2. 

0.42 

1. 

1556.9 

0. 

0.05 

0.00 

180*00 

21. 

79,67 

119.10 

40. 

26.0 

0,04 

42014. 

Z. 

2. 

0.42 

1. 

1556.6 

0. 

0.46 

. . ft^ftft. 180,00- 

22,— 

—7.9 ',68- .149*10 

- . 5ft* 

24*7- 

_ft*ft4_- 

.42231*- 

3*— 

3*- 

ft* 41- 

1. 

1554.1 

0. 

0.07 

0.00 

180,00 

22. 

79,68 

119.10 

52. 

26.2 

0.05 

42279* 

3, 

3. 

0.41 

1. 

1553.6 

0. 

0.07 

0.-00 

180,00 

22.. 

79,68 

119.10 

54. 

26,3 

0.05 

42327. 

3. 

3. 

0.41 

1. 

1553.0 

0. 

0.37 

b.co 

iao-,00 

22. 

79 ’.69 

119.10 

56. 

26.3 

0.05 

42375. 

4 . 

4 , 

0*41 

2 . 

1552.5 

0. 

0.07 , 

0.00 

165-,22 

22. 

^69 

119.10 

50. 

26.3 

0.Q5 

42403* 

4, 

4. 

0*41 

2. 

1551.9 

0. 

Q.38 

0.00 

1651,22 

23. . 

^69 

119.10 

60. 

26.4 

0.05 

42411. 

4, 

^ t 

0.40 

2 . 

15.51.4 

0 . 

-UAjia 

0*00 

16 S , 2 ? 


-...is-Us-- 

119.10 . 

* 62. 

? 6.4 

_J1*Q6 

42419. 

5. 

5* 

0*40 

2 , 

1550.8 

0 . . 

Q.08 

3.00 

16S;22 

23. 

79 -,69 

119.10 

65, 

26.4 

0.06 

42429. 

5. 

5, 

0.40 

2 . 

1550.3 

0 . 

a. OB 

0.00 

162.43 

23. 

79^69 

119.10 

67. 

26.-4 

0.06 

42434 ..., 

.. .5, 

5 . 

0.40 

2 . 

1549.7 

0 . 

Q.U 

0.00 

142., 65 

25. 

79;70 

U9.10 

06. 

26.5 

0.00 

42500.- 

. 9, 

?, 

0.40 

3. 

1544.7 

0 . 

0.11 

0.00 

.147,65 

25. 

79;?0 

119,10 

90. 

26.6 

0.00 

42524. 

9, 

10. 

0.39 

4. 

1543.6 

0 . 

0.12 

0.00 

152.10 

26. 

79;?0 

U9.10 

94. 

26.6 

0.08 

42540, 

10. 

11. 

0.39 

4. 

1542.5 

0 . 

—ft. 12 

0 .00 

lAl*-63 . 

... 

- 29 Zft- 

U4*lfl- 

91* .26*6. . 

.-0*09 

425S6.*,.- 

— 11*. - 

. .12* 

0*39 

5. 

1541.4 

0. 

Q 

0.00 

141,63 

27. 

79 ; 70 

119.10 

103. 

26.6 

0.09 . 

42572, 

.12. 

13. 

0.39 

5. 

1540.3 

0. 

0.12 

0.30 

.145,94 

27. 

79.70 

119. 10 

107. 

26 »7 

0.10 

42509, 

- 13. 

14, 

0.39 

5. 

1539.2 

0. 

0.14 

0.00 

137,96 

2S. 

79.70 

U9.10 

112. 

26,7 

D.lG 

42596, 

JL4, 

15. 

0.39 

6 . 

1530.0 

0. 

0.14 

0.30 

137,96 

28. 

79i?0 

119.10 

116. 

26.7 

0*10 

42601. 

16. 

16. 

0.39 

6 . 

1536.9 

0. 

0,15 

0.00 

141.63 

29. 

79,70 

119.10 

120 , 

26,7 

0.11 

42606, 

-17. 

17. 

0*39 

7. 

1535.8 

0. 

ft* IS. 

0.00- 

- 130,.9.1 

35,*.. 

--71*71 

.U9*i0_, 

159.. 

8.0— P.Q t.1.4. . 

-42650. . 

29, 

. . 30. 

Q*3ft . 

.. -11. 

1525.7 

D. 

g.2o 

0.00 

132,37 

37, 

79 ; 71 

119.10 

169. 

27.0 

0.15 

42660. 

32. 

33. 

0.30 

13. 

1523.5 

0. 

0.21 

0.30 

126.74 

38. 

79.71 

U9.10 

176. 

27.0, 

• 0.16 

42670. 

36. 

37, 

0*38 

14. 

1521.2 

0. 

0.22 

3.00 

126,29 

40. 

79.71 

119.40 . 

185. 

.27,1 

0,17 

42600. 

39. 

4l. 

0.30 

15.' 

1519. Q 

,0. 

0,23 

0 . 00 

126,08 

42. 

79,71 

119.10 

194. 

27.1 

0*17 

42690. 

43. 

45. 

0.38 

17. 

1516.7 

0. 

y *24 

0.00 

125,60 

44. 

79;?2 

119*10 

203. 

27,2- 

0,18 

42700* 

-47.. 

49, 

0*37 

16. 

1514.5 

0. 

-.0*25. 

0.00 


■ 46. 

—29472 

U9-*10-. 

211*- 

-27..2- 

-0*19 - 

42710*-- 

51*-. 

S3. 

0*37 _ 

20 . 

1512.2 

0 . 

0.26 

0.30 

125,22 

48 , 

79;?2 

U9.10 

220 . 

27,3 . 

0.20 

42720. 

53* 

57. 

0.37 

21. 

1510.0 

0 . 

C.27 

0,00 

123,62 

• 50. 


119,10 

229. 

27,3. 

0*21 

42730. 

■ 60 * . 

62 , 

0,37- 

23. 

1507.7 

0 . 

0.36 

o.co 

121,51 

■24, 

79,72 

119,10 

30l. 

27.7 

0*20 

42021* 

107. 

113. 

0.35 

40. 

1487.6 

0 , 

0.36 

0.00 

120,09 

Bo, 

79,73 

119.10. 

326. 

27.8 

0.29 

42041. 

120*- 

126. 

0.35 

44 . 

1463.1 

0 . 

0.43 

0.00 

119,81 

87. 

79,73 

119.10 

344. 

27.9 

0,31 

42661. 

133. 

141, 

0.34 

48. 

1478,6 

0 . 
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WINqE 


oooooooooooooooooooaoooo^ooooooooooooo 



TIME 

ftAilGE 

BEARINO 

altitude 

FIT/EL- 

TL-T/AZ- 

velocity 

ACC. HACU- 

TJIfiUST 

DRAG 

0/PRES 

-ETA .OALrA 

weight 

WINON 

0+42 

0.00 

120,02 

94. 

’ ' 79 ; 73 
-79113 
70;73 

119.10 

362. 

25.0. 0,32 

42 V 82 + 

147.“ 

156. 

0+34 

52. 

1474,1 

0. 

U .44 

O.OQ 

120.24 

loi. 

'119.10 

. 380. 

28.1 0,34 

42902+ 

161. . 

171. 

0.33 

57. 

1469.6 

0. 

0 . 4 6_ 

. 0.6-fl 0 

0,00 

. 119+5,5— 

- - - 

119.10 

388. ?R.g n.36 42123+ 

171+ 1M+— 

n,3g 

..61+ . 

1465.1 . 

-...0, 

C,48 

180,14 

X17. 

79;75 

149,10 

417. 

88.3- 0+37- 

4*443. 

193+ . 

. 2C6-, 

(U32 

65. 

1460.6 

0. 

0+50 

0.00 

119,90 

183+ 

-79,73 

119,10 . 

.435, 

,88,4 . 0+39 

4*964+ 

20 9+- 

-224+ . 

0.31 

69, 

1456.1 

0. 

0f52 

0,00 

119,73 

134. 

79,73 

U9.10 

453, 

28,5 - a+4i. 

42954+ 

226+. 

243+ 

0.30 

72. 

1451.7 

0. 

(j+70. 

0.01 

119,49 

A29. 

79,71 

-119.11 

621. 

29,3 .0,56- 

4.31Z1+ 

412+.. 

455. 

-5+16 . 

83. 

1411.3 

0. 

0+74 . 

0.01 

119+43 

254. 

79;xi- 

119.11 . 

658, 

89,5 . 0+39. 

43213+ 

461... 

511,. 

0.15 . 

76. 

1402+3 

0. 

e .78 

n + rn+. 

lli+AO- 

119,39 

1R1 , 

79, 7fi 


69* , 

2«,7 0,69 

.43255+—. 

■513. 

577. 

0 ■I'' 

. 05 + 

1/193.3 

n. 

0.82 

.0.01 

S09. 

79 68 
...79+67 

us! 11- . 

735. 

29.9 d+66- 

43297+. 

567+ 

j636. 

0.08 

51. 

1354.4 

0. 

0.86 

0.01 

119^54 

639.- 

119. U- 

. .773. 

30+1-0+69 

63339+ 

.623+ 

J. CUL 4 +. 

.0+05 

32. 

1375+ 3 

. 0 , 

0+94. 

0.01 

119.,24 

*05. 

79,64 

ll9.il. 

851.- 

04-J!6+ 

43423+ .. 

261+ . 

551+ 

0+02 . 

14. 

1352+1. 

Q. 


0.9i . 
0..98 

-i.ao.- 

2*17 

J*C5 

3.52 

4 *C7 

5*U3 

6*;.S- 

7* u7 

ii*C3 

9*i;3 

10 > Ci3 

11 .03 

■ 4»-r*7. 
13,15 
14-^03 

15.03 

16.03 
17. OC 
•ta.oi 
19.09 
2(1. 02 
21.11 
22 . 01 
23.52 


o.tri 

OvOl 

-(1,01 

0.07 

0.14 

0.19 

0.24 

0.34 


119 10/ „ 
119.39 
11S-.-41-.. 
119*18 
139^l7 

119.14 

115.14 
119,16 


—0,40 -_4rl.a.^lB-. . 
0.52 115,19 


420* 

437. 

2*02. 
4543, . . 
4876. . 

ta32*.- 

1*164, ... 
-1 »1 6. — 


0.6o 

g.69 

0.77 

0.05 

-0.-94 

1.01 

1.07 

1.15 

1.22 

1.29 

1^37 

1.51 

1.68 

1.92 

2.11 

2.44 


119.21 

119.22 

119.23 

119.24 
144i-*26 
119|27 

119.29 

319.30 

119.31 

119.32 
419^,34 
119,35 
119,39 
119.42 
119.44 
119,47 


3*519. 

19*38. 

2*553. 

7S?66. 

2**89. 

29*55,. 

3*947, - 

32*89, 

34*13. 

34*77. 

38*31. 

44549. 

48953. 

5SU5. 

66*65. 

68043 . 


.79.63 

Z9.UJ. 


29;£3 
J9^48 
. 79., 13. 
. 79 ' 03 . 
- I b.,- 9 - 2 


il9,U 871,. 30.5 d,.7A 

119,12. . 89i; 30. 6- 04.aa. 

.U9.1B . -3209, 2*.,4. .2.t-42- 

U9..19. 3308, U-6. J.,.03 

4.19,21.. 3198, ,5.9 1^95 

.U9...23, . 3028. a. 1- 2^.82. 


m 6*.i. . . -S35.1 .A9 Qi- 

43.445+ . ..AOS*. .931+ 

4H4R6, 074 , 07?, . 

■464831.-5510.. . 5239. . 
33 054+ . .9501.— 10693, 
.638+ .9.445+ lQiL44 . . 
. -0+. -1319. ..2606, 

0-t. .1109,.- 7554+. 


O+.03 

0.04 

-O+OS— 

.0.-00.. 

-o+ioa 

- 0 + 00 - 

0+05. 

-OvOO. 


26. 

37. 

-42,. 

12 . 

12 . 

а. 
7. 

б, 


1352. 6 
1345.1- 
174.3.5 
1070.8 
869.7 

031.0 

267.0 
267.0 


0 . 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


7a;j9. 
78.117. 
20,54 
78,. 40 
7 8 2*. 
7 8 ,05 


-Z7;i0 
- 7X1 75 
77;S7 
77,35 
77., 19 
-7X,oi 


ll2,-43 2231^- 


7.6 59 
76,88 
7617.5 
X6.,70- 
.-76; 60 


149 r46., 

-119., 15. 
.119,51 
119+53 
119.56, 
1 19,» . 


S«2 -2419 


0,.. .3*1, 


- 1294 + 

*019. 


- 0 + 00 - 

OrOO 


3.. 

3. 


.26Z.-O 
267,. 0 


- 0 . 

0 . 


119., 41- 
U9L,65 
119.66 
112,65. 

119,21 


HIN05 

0 , 

0 . 

0. 

0 . 

0. 

0 .. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 

0. 


’■1*1 9* 

yfl»1 , 

4.4 

7 , 70 

0 1 

9^0- 

—6425+— 

.-0+00- 

5, . 

26.7+4 

267.0 

0. 

0. 

119.29 . 

2734; 

3+9 

• f ^■V"- 
2+60 . 

- -0+- 

755+- 

. J319 

0+00 

5. 

0. 

0. 

119.32 

2619. 

3.5 

2+52 

-0 + 

. 674. 

-449.0+ 

0+00 

5. 

267.0 

D. 

0. 

119,34 

2513, 

,3+2 ,2+44 

. . 0 + 

-579,- 

3794. 

0+00 

4. 

267.0 

0. 

0. 

U9.37 . 

2416,- 

2.9- 2+37 

.. -0 + 

499+ 

3228+ 

-0+00 

4, 

267.0 

0. 

0. 

119,40 

23*9. 

2.6 .2,30 

o«- 

432 + 

27j12, 

0,00 

4, 

267.0 

0. 

0. 


0. 

0 . 


2107 . 

2063 , 

1983 . 

1929 . 

7 R 4 A 

2.1 
. 1.9 
1+8 

4+7 

,1 

2 + 14 - 
2+10 
2 + 05 ■ 
1+99 
2,91 

_ . . .0+. 
oi. 
0 + 

-• - 0 +- 
010 ® r 

-. 257 .+-. 

. 252 + 

221 . 

194 , 

778 , 

4252 , 
1552 , 
1348 . 
1167 , . 

J 774 . 

0 + 05 
O + OD 
0.00 
0+00 
0^ 9 n 

3 . 

2 .’ 

2 , 

2 , 

2 + 

267 . 0 , 
267.0 
267.0 
267.0 
- 214+1 — 

0 . 

0 . 

0 . 

0 . 

_ 4 l+ . . 

0 . 

0 . 

0 . 

0 , 

0 . 

4143 . 

43 ;s- 

4+28 

- 8 . 9 il+. 

645 +. 

. 4054 ;. 

. 0+66 

3 . 

191.9 

0 . 

Q, 

5654 . 

48+0 

. 5,54 

-. 5031 + 

54.4 + 

6116 . 

0+00 

3 . 

1 S 2.9 

0 . 

0 . 

5822 , 

0 + 8 - 

- 6 + 04 - 

. .458 + . 

. 422 + 

4829 . 

.. 0+00 

2 . 

148,1 

0 . 

0 . 

5808 .. 

0 + 3 - 

- 6 +OU. 

. 435 +— 

. 330 +. - 

- 3376 . 

-0 + 00 

1 . 

140.0 

Q. 

0 . 

5664 , 

2+a_ 


- --- 0 +- 

213 + . 

2386 + 

— 0 . 

0 . 

148.0 

0 . 

0 . 


- l l 



TIME 

RA2GE 

BEARING 

autttude 

FLT/El 

flt/az 

velocity 

ACC 

MACH 

thrust 

DRAG 

D/PRES 

ETA 

OALFA 

WEIGHT 

windn 

winde 


2.67 

119.49 

7A532. 

76.53 

119.74 

55B6. 

2.1 

5.74 

0. 

173. 

1759a 

0. 

0. 

148.0 

0. 

0. 

25i0< 

2,77 

119,50 

76932. 

76,50 

119.75 

5554. 

2.0 

5,69 

0. 

151. 

1533. 

0. 

0. 

140.0 

0. 

0, 

E£*ii4_ 

.1,11 

U¥,52 . 

. 81127. 

7.6,33 

U9.Z9 

9-462- , 

1-.^ 

-.3*57 

.. a*_ 

..jZiA.. 

. ..212a 

.„0* 

0. 

. lAfl.O 

Da 

0. 

27.26 

3.24 

119,53 

88884. 

76,34 

119.00 

5430. 

1.5 

5,52 

0 . 

83. 

826, 

0. 

0 . 

148.0 

0. 

0. 

30^00 

3.98 

119,58 

1C7179. 

76.09 

U9.08 

9279. 

1.2 

5.29 

0i 

33. 

330, 

Oi 

0. 

148.0 

0. 

0. 

32.06 

. 4.24 

119,60 

118907. 

75,99 

119.91 

3232. 

1.1- 

. 5.10 

Q, 

24, 

237 , 

0. 

0. 

148.0 

0. 

0. 

33.04 

. 4.44 

119,61 

118564. 

75.92 

119.93 

5198. 

1.1 

5.10 

Q. 

19, 

196, 

0- 

0. 

148.0 

0. • 

0. 

35,92 

5.04 

-ii9,65 

132940. 

75,69 

120. QC 

5100. 

1.0 

4.89 

Q> 

10. 

23, 

0. 

0. 

140.0 

0. 

0. 

n ■ 6 ^ 

5.17 

-13,9.. 66 

136896. 

.. iS-Ul. 

1 ?n , pi 

-- 5n6Q , 

-4.*0- 

4,84 

(1 - 

8.1 . 

.. --aii- 

-—II... . 

. . n. 

J-49.-D 

0. 

..JIa. 

37,20 

5.30 

U.9,66 

135X38, 

75,59 

120.03 

5059, 

1,0 

4.80 

0. 

7. 

70. 

0, 

0 . 

148.0 

0, 

0. 

3B.48 

5.56 

U9.68 

149483, 

75,49 

120. Q6 

5018, 

1,0 

4.71 

0. 

6. 

53, 

Oa 

0 . 

148.0 

0. 

0. 

45,52 

7.01 

115,76 

178880. 

74,88 

120.21 

4798; 

1.0 

.4,47 

0* 

1. 

13, 

0- 

0 . 

148.0 

0. 

0. 

46.16 

7,14 

119,76 

186038. 

74183 

120.23 

4778, 

1,0 

4,46 

Q. 

1, 

-12. 

0- 

0 . 

148.0 

D. 

0. 

47.44 

7.40 

119.78 

187914. 

74,71 

120.25 

4738, 

1,0 

4.46 

a* 

1, 

9. 

- 0* 

0 . 

.149.0 

0, 

Q. 


_ 2..5J- 

.11SU22.. 

1 , 

74-61. 

■122., 22_.. 

- -4I1E— 

-I -.0 ■ 

.4, .45. 

IU_ 

li 

a.- 

lU- - 

.. .0- 

— lii.a 

Oa 

-Da 

53.84 

e.7o 

119,85 

216322. 

74.09 

120.39 

4943 ; 

1.0 

4.50 

04 

Q. 

3 , 

0. 

0. 

148.0 

0, 

0. 

55.12 

8,96 

119.86 

221089, 

73 ',97 

120.42 

4904. 

1.0 

4.53 

Da 

Q. 

2, 

0. 

0 . 

140.0 

0. 

0. 

56.40 

¥.22 

119,87 

227804. 

73,63 

120.44 

4466, 

1,0 

4,55 

Oa 

.. 0, 

2. 

0- 

0. 

148,0 

0. 

0. 

57.66 

9.48 

119,89 

232069. 

73.70 

120.47 

4427, 

1.0 

4,50 

Qa 

0. 

1. 

0- 

Q. 

148.0 

0. 

0. 

53.96 

9.73 

119,90 

23BE82. 

73', 56 

120.90 

4388. 

1.0 

4.61 

Oa 

0. 

X 1 

0* 

0. 

1.48 , 0 

0. 

0. 


-11,05 , 

. latuox 

..216999. 

. Z2;22 

12n-Zi 

, -.4a41-,. 


. 4.-57.. 

- 2-a 

. -JL_ 

Q-* 

Q.^- 

a. 

149LlD 

Q. 

0. 

73,04 

12.57 

120,04 

296652. 

71.89 

120.79 

3967. 

1,0 

4.4.9 

Oa 

- .0, 

0. 

0. 

0 . 

146,0 

0, 

0. 

75.6c 

13.08 

120.06 

Sp4902. 

71,54 

120.04 

3891. 

1.0 

4.49 

Qa 

0. 

0. 

-0. 

0 - 

148.0 

0. 

0. 

70,16 

13.59 

120,09 

313B49. 

7lUs 

120.69 

5815. 

1.0 

4,49 

04 - 

0, 

Oa 

0- 

0 . 

148.0 

0. 

0. 

00.72 

14 .10 

120.11 

322492. 

7o;bi 

120.94 

3740. 

1.0 

4.49 

Oa. 

0.- 

. 0. 

.0- 

0 . 

148.0 

0, 

0. 

¥6.06 

17,14 

120,25 

373892. 

66,22 

121.22 

3292. 

1.0 

4.49 

Oa 

0. 

0, 

Q- 

0 . 

140.0 

0. 

0, 

.?a,.64_. 

.-17.65- 

-120.28 . 

■380017. 

■ 67.,.22- 

_i 2 i-az_. 

.. -32l4-_ 

< .n 4-49 

- - n , 

n , 

n . 

n, 

0 . . 

1 4ft . n 

0. 

. .. Da 

161.20 

10.15 

120,30 

388940. 

67 19 

121.31 

3145. 

1.0 

4,49 

Oa 

Q. 

Qt 

0. 

0 . 

148.0 

0. 

0. 

163,76 

10.66 

120.32 

399661, 

66^64 

121,36 

3072. 

1.0 

4,49 

Da 

0. 

Oa 

0. 

0 . 

148 . 0 

0 . 

0 , 

106.32 

19,16 

120,34 

4(,270O . 

66,06 

121.40 

3000, 

1.0 

4,49 

0, 

0, 

0 

0 . 

o' . 

148 . 0 

0 . 

0 , 

129,36 

23.67 

120.54 

457797. 

59,26 

121.70 

2367. 

1.0 

4,49 

0, 

0. 

-Q, 

Oa 

0 . 

146.0 

0. 

0. 

134.46 

24.66 

12C!58 

467015. 

57;26 

121,05 

2233. 

1.0 

4,49 

Ui 

0. 

0, 

0. 

0 . 

148,0 

0. 

0 , 

132,6c 

25,66 

120 -.62 

472331. 

54-56 

121-93. 

. -2102— 

-IJl- 

A-A9. 

... ._iU 

— Jll . 

(U- 

0- . . 

. . 0.- 

0 

0 . 

0, 

144.72 

26.65 

120,66 

485446. 

52', 40 

122. QQ 

1976. 

l.Q 

4,49 

0, 

0. 

0 1 

Oi 

0 . 

148 . 0 

0 . 

0 , 

149,04 

27.64 

120,70 

493361, 

49,46 

122.07 

1854. 

1.0 

4,49 

0, 

0, 

0, 

0- 

0 . 

146 . 0 

0 • 

0 , 

200.05 

37.33 

121,04 

52S430. 

■ 5 , 66 

122.64 

1207. 

0.9 

4,49 

0,. 

0. 

0, 

0. 

0 . 

146. 0 

0 . 

0 . 

2p9.05 

39.07 

121.09 

523028. 

<-18,10 

122.71 

1264- 

0,9 

4,49 

Qa 

0. 

-0, 

0. 

0 . 

148 , 0 

0 . 

0 , 

227,05 

42.54 

121 il9 

511132. 

>38.00 

122.05 

1526. 

0.9 

4,49 

0, 

0, 

0, 

0. 

0. 

146.0 

0. 

0. 

■ 


r « - 

— 

- ■ 

• 

— 

— 

— 





— 

— ... 

- 
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17 



-TJME 

' 4 < 

ha.mge 

BEARUJG--4WT-1JUDE 

FLT/EI 

fLT/AZ 

VEUDCITY 

ACC 

HACH 

THRUST 

ORAG 

D/RRES 

eta 

QALFA 

WEIGHT 

windn 

iiNnE_' 


49.53 

121.37., 

4Sg?ia. 

-59,25 

123, Q2 

2373. 

l.Q 

4.4? 

0, 

Q, 

o', 

0. 

0, 

148.0 

0, 

oi".'. 

365., 13% 

70. Q3 

121,67 



-76,26 

128,89 

926?. 

0.5 

5.30 

0, 

78, 

771, 

0, 

0. 

148.0 

0. 

04 

3.67^05- 

. 2IU44-. 

12J.,67__ 

..msia.,. . 

»76,4n 

.122 ,.86. 

... .5307. . 

^Q.1.2.. 5 » 4.2 

_Q. 

-_12fi.,_. 

-.-laZii.. 

- 60 -. 

. .0. 

148.0 

0._ 

„ 0o._ 

368,05 

70.64 

121.67 

70759. 

-76,46 

122.87 

5307. 

0.3 

5.45 

0. 

16^, 

1632, 

0. 

0. 

146.0 

0 . 

Q. 

369i05 

70.84 

121.67 

66602. 

f76,5S 

122,87 

9297. 

0.5 

5,46 

0, 

210, 

2Q?lr 

0. 

0. 

148.0 

0 . 

0. 

37C-55 

71.15 

121,67 

06695. 

b76,65 

122.86 

5262. 

1.1 

5.44 

0, 

301, 

2992, 

Q.. 

0. 

148.0 

0. 

0, 

371,05 

71,25 

121.67 

56940. 

bZ6;69 

122.86 

5242. 

1.3 

5.42 

0, 

339. 

3356, 

0. 

0. 

148.0 

0, 

0. 

372,05 . 

71.44 

121,68 

56261. 

b76,76 

122.85 

5191. 

1.9 

5.36 

0.. 

424. 

4195. 

0. 

0. 

148.0 

0, 

0, 

iZliJlS.. 

__21,64. 

.1214 69 . 

ilMiiU.. 

o76, 83 

_ 122,84 

Jll?-... 

■ 2.6 ■ 


Hi. 

.__527.^ 

51 67 . 

__Jla 

. n » 

146,0 

0. 

0. . 

374,05 

71.62 

121.68 

40899. 

v76,?1 

122.84 

5023. 

3.4 

5.19 

0, 

647, 

6327, 

0. 

0. 

148.0 

0. 

0, 

375,05 

72.01 

121,68 

39464. 

876,96 

122,83 

4699. 

4.3 

5,06 

Ot 

762, 

7987, 

0. 

0. 

l4fl.O 

0, 

0. 

3Z6i,05 

72.19 

121,68 

39?63. 

i77,Q6 

122.82 

4745. 

9.2 

4.65 

0, 

-920, 

0713, 

0. 

0. 

148.0 

0 . 

0. 

377.05 

72.36 

121,68 

26324. 

-77,14 

122.82 

4563. 

6.1 

4.57 

Ot 

1044, 

9S3ii 

0. 

0. 

148.0 

, 0 • 

0. 

576,05 

72,52 

121,68 

24674. 

*77,22 

122.81 

4355.. 

6.S 

4,28 

0) 

1152, 

.10157, 

0. 

0< 

148.0 

0 . 

0. 


....22^62 

121^68 

. -26234- 

.-»274.31_, 

_i2z,ao_ 

4126- - 

7.4. 


n. Ig36. tn523. 

--.Qa 

-0- .. 

146. Q 

0, 

0. 

360 <05 

72.6? 

1?1.68 

16628. 

-77,40 

122,80 

3881, 

7.6 

3,70 

Q< 

1304. 

10612. 

0. 

0. 

148.0 

0 . 

0. 

361,05 

72.95 

12ll68 

13164. 

.77,50 

122.79 

3629. 

8.1, 

3.40 

0. 

1336. 

10430. 

0. 

0. 

108. 0 

0 . 

0. 

362.05 

73.06 

121.66 

9752. 

«77;60 

122.76 

3365. 

Q.l. 

3.12. 

0, 

1337. 

10015i 

0. 

0. 

148.0 

0 . 

0. 

363,01 

73.19 

121.68 

6686. 

i77-,7l 

122.78 

3115. 

B.O 

2.66 

Q< 

1322. 

9442. 

0. 

0. 

148.0 

0 . 

6. 

364,01 

73.29 

121,66 

8765. 

-77 • B3 

122.77 

2863. 

7.7 

2,60 

0. 

1267. 

87114 

0. 

0. 

148.0 

0. 

0. 

i8S-,0t. 

, -2i,39. 

121,62. 

- 4U45-.. 

..b77v97_ 

.. 122 ^ 21 .. 

2421-:- 

..TUi _2-36_ 

lU 

■_1225. 

. 7906a. 

_a. 

. 0 . 

148.0^ 

0. . 

0. 
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EL 

time 

. A1 I 

range - 

. TIME.. . 

-ALT 

RANGE 

VEL . 

FLT/EL. . 

FLT/AZ 

(DEG) 

(SEC) 

)?. 

196'. 05 

trT> 

(NM) 

(SECl 

(FTl 

(NHJ 

(FT/SEcl 

(DEG.) 

(DEG) 

80.00 

525i76» 

36.6 

22.Qfl 

60165. 

2.1 

g.e.09r36 

76,70 

119,63 


p A_ C, I. ... 

a 

time' range AZ .._1 . 
(SEc) <NM) ■ (DEG> 


385. 4A 73.4 121.6?. 


THIS concludes the CALCULATION, 






APPENDIX D 


CONTROL CARDS 


Deck set up with object program on tape. 


1 

8 

16 

0 

IDENT 

153202, SENS-5D 

0 

EXECUTE 

DUMP 

0 

LIMITS 

100, 30K,,5K 

0 

TAPE. 

R*, XID,, 

0 

FILE 

02, X2r, ~ 5 i : 

0 

FILE 

03, X3R, 5L 

0 

DATA 

05 

0 

- INCODE 

IBMF 

Input 

data with 

NAMELIST Format; 

0 

DATA 

01 

0 

INCODE 

IBMF 


Input data with FORMAT (V) Format 


^ ENDJOB 


* **eof 



Deck set-up with object program on cards. 

1 

■ 8 

16 


IDENT 

153202, SENS-5D 


OPTION 

FORTRAN 

Object 

Deck 


0 

EXECUTE 

DUMP 

0 • 

LIMITS 

100, 30K,, 5K 

0 

PILE 

02, X2R, 5L 

0 

FILE 

03, X3R, 5L 

0 

DATA 

05 

0 

INCODE 

IBMF 


Input Data with NAMELIST Eonnat 


$ DATA 01 . 

^ INCODE IBMF 

Input Data with FORMAT (V) Format 

^ ENDJOB 

***EOP 



APPENDIX E 
SOURCE PROGRAM 



label main 


O O 
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1 

2 

3 

4 

5 

6 
7 
e 

9 

1C 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 
29 
3n 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 


CHAIN 

C*«<to SSNB**50 t #•••*• **•••###*•»#»**####**•?*•#***•#* *#••*#»*»••* 

C AUtRoBS sa s; P. STNGH AND R, GOOKi 
C date — sJUNE 23, 197?; 

C 

C * THIS PROGRAM CALCULATES THE 'FOLLOWING •- 
C (1) SUHmaB* at BURSiOuT, APOGEE, AND IHPaCT IN '0A!», 

C (2) F<7) and DF{7> curves JN 'FZCURV'i 

C (3) UNIT-WIND BFFECTjCORIOLIS DEPLN AND RANGE DERIV IN lEFFECT*, 

C 

C » the following units ARB ADOPTED — 

C (1) INPUT/OUTPUT *- FoOTaPOUND-sEcDND-DEGRfcE, 

C EXCEPTtOEl -o surface RAfJfiE iN NaUTICAL MILES IMmI, 

c {2> CALCULATIoSIAL STEPS -- FOOT-PoUnD-SECOND-RADI AN, 

C 

C » V,IND PROFILE ——IT _ . 

C (1) NoRTHpTo-SOUTH DR SAST-TObHeST DIRECTION IS POSITIVE, 

C (2) NO WIND MEANS NO WIND FIELD AT ALL. 

C (3> CONSTANT MInD HEAnS ONE Wi?gD FIELD FRQM GROUND TO SOME 

c ALTITUDE. Parameters used are (hind, hnAz, alow, ahighi, 

C denoting wind speed* NOSTh aZiMUTHt minimum altitude, 

C AND maximum AlTITUPJ RESPECT tVSl^,. 

c (4) variable wind HgANS WINQ VELOCiTJf VARIES WITH ALTITUDE, 

C parameters USED ARE inlev.alTw^speed'.'dIr) , denoting 

C NUHSEB OF altitudes, ALTITUDES* SPEED, AND NORTg 

C . AZIMUTH MEASURED CLOCKWISE, 

C 

C * THE following are THE LIMITS FOR ENTRIES •• 

C (1) 20 PAY PAYLOADS, 

C (2) 61 elevation ANGLES, 

C (3> 30 THBUSTING aNd COASTING PHASES, 

C (4) 20 SPENT STASES, 

C (5) 99 ENTRIES IN THRUST AND MACH Ta0LEs“ 

C (6) 9? WIND altitudes: 

C 

C * THE FOLLOWING TAPES OR FILES ARE DESIGNATED 
C (1) ITAPEJ ULWaVS |5) IS USED FOR READING DATA IN NAMELISTS, 

C (2) ITAPE2 IS used fOR READING VEHICLE THRUsT-HACH TABLES, 

C (3) ITAPE3 IS USED flNTERNALI WRITING AND READING THRUST TABLES, 

C <4> ITAPE4 IS USED UNTEBNAUI Fop SpENt STAGES, 

C 

e*«*o »#«•»•»*•«••*»«»«•••••«««•»•••••»«**»•*«•«««••••*«••«««««•••«•»*•«« 


MAIN 


PAGE 
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label main 


Page 


z 


AZ 

44 

45 

46 

47 

48 

49 

90 

91 
9Z 

93 

94 

95 
S6 
97 
SS 
39 
6 g 
61 
62 

63 

64 
85 
66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
7B 
79 
ap 
61 
§2 
63 
84 


C**#* 

integer RKSfBP 
integer ordbb 

REAL UATL»LQNGt,LATiLbNG 

REAL MASS, MPPdP.HPLfHDOt ♦MACH ,TcSVCcff7TY¥ 

REAL IPTflPAMPR 
REAL MSEP 

double PRECIIIQN DPsIDPC 
CHARACTER NAME»88,i:£Pf*6l,RlGfclT»i9 

C ■ ' 

DlMENiiON segED(ioO)jDlRU60> ' 

C 

COhMoR/STOROO/ ZERb,0NE,PIiOPl»HP.iiTPi;cP(9) 

COHMOR/STOROI/ 66,Ga;0HEGA;GH7^E«'R0»ESQ,BSnl«ESQl 
COHMON/STORP}/ pRESp;DENSOTTeMPO,90UHDO,PRESP<.9)»TEMP8{9(, 

H»LTBt9).TSLOPE<9> .BONSt,PR6S? DENS,TEHP«SOUND 
COmMON/SToR'cS/ X0,roIZ6;R,RXt7ALt»3R;yLlfl 

C0HMON/STOR047 GDLilTLVUATLiLONGL.THESALil^HiL^QDLATiLAT, LONG, THETA, 
, PHUSAHAL*PSIL,GAMA*PSi;ELQbL,AZQ0L7ELQD,AZGD,EL, 

az.sr*: 

COMMON/STORp5/ VT ( 3)IVA < 3) iVELf ACC^ VgLIMfVL'l.H 

COMKiON/gTORgS/ Rtli?) , ' 

C0MM0N/STOHp77 T;DTjTNPi;DTUA>DTAIfTao^RKSTE^,EPBIG,BPflNYi 
, DHlR,QHAXiKSTBP,NQO0D 

COMMON/SToROa/ Stl4i;P(l4,5);tlR8(l4)lDPC)14)tERW<14);sP(6,20i 
COMMDN/STOR09/ MASSrMPRoP p MPL i SUI h 632 | lybof i THR, HAC hJcA , AD , oS, 

SLOP5;CHP*LO«;i,CP,|tYJsytSZ , ‘ 

COHMON/STORlp/ HE. HN.'HLa9T,HINP. WNAZilAL 0W, AHIQH, NLEVi^ALEVUOP i , 

• ALTMttb0)iVILltl60tp:v^LN(i0e) 

COHMON/STORU/ KPRBV,'HPREV>«PBEV,AfiREv;LPREV<14),yPRBVll4(, 

I 08DIMH4) , ■ . 

CQHH0N/ST0R13/ nPmAXInSmAX-, NPL,NANOjtNVAR*NP,NPl» JSPERTi jSi JSC, NCON 
COmMON/SToR13p' APTJAPA2,APA,ARRiBoT,0OA,Bofi,'4oVLiBOE|:,BQAznRT. 

'cOHMON/STORia/ lTtht;iTAPE2,lTAPE^~,TTAPl4Pri’’RW,lTr^ 
COMMON/STOR9T/ NAmE| ARRA Y { 8 00 ) 

6 

COMMON NPST,NS6P;N5li;NH(14J,KSP(20},"MASS,f3B,TAX{3),ARG<61),B(4), 

, PLH(20)7pHTiS2l,TSER{S2>;MsEP(32>;TSPl2D),SPHt80l,A(9t,14> 

C 

equivalence UORAYtli,SpeED<lh;uTRAY(iOi)',"fiR(l5r , 

C 
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LABEL MAlN PAGE , 3 


65 

66 
87 
68 
69 

90 

91 

92 

93 

94 

95 

96 

97 
96 
99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 


113 

114 

115 

116 
117 
116 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

C***» VALUES oF C 0 NSTA 8 TS, 

data zero, one/ l.oEirlS. 1.0E*00/ _ 

DATA PI,QPI,'WPt,Tp!/ 3,14159265* 7 » 85398 l 64 E' 01 » l, 57 o 79633 i 
, 6 , 28316531 / 

data ESQrESQl#es&I/ 6;693421 ooE-03i 9,933o6579E"01i 1,00673853/ 
data GOiOHEGAiOMiBEiRQ/ 32fl74» 7.292llE*o5* 1 '. 4076576E+16 1 
; 2 ,O 09 O 578 OE*O 7 , 2 , 09257410 E* 07 / 

data CF/ 57 |? 957795 iI i', 74532925 E- 02 * 6 , 0761 t 550 E *03 J 
, 1 . 64578635 E 504 , 3 , 1060 ^? 70 E- 02 i 0 t 3 “ 46 » 3 ; 2 'oS 3990 * 

■ , 1.6521 5.S9956804E-01/ 

data NPMx'»NSMAx;NVAR/ 30» 20* 99/ 
data dummy/ 99999,0/ 


C 

C**«* U, S, 1962 standard ATMOSPHERIC T*BlE, 
data C0NST/3A,16$19ABS3/ 

data PRESo/2U6;2/;DENSO/2,3769Ei3/;TEHP(l/28e,l5/,SOUNDO/lll6;45/ 
data PRESB/ l,O.2;233dlE»l|5,4O328B-2,0,96663E-3»l,o9465Eo3i 
; 5, 822896-4, l,797l8E-4Ji,0a4iE-»5il,6223E»6/ 

data temps/ 388, 15*'216, 65, 216, 65, 828, 65, 270,66*270. 65, 252, 65, 

. 180 , 65 , 180 , 66 / 

data halts/ 0,0,il.oE«3.20;oE*3;32ioEi3,47tgE*3,52;o^',3,6l',OE*3, 

; 79 , 0 E* 3 , 88. 7436 * 3 / 

data TSLOPE/ -6a6E-S»C>0*liOBn3*2,8E*'3«O,Oi*3iOE**3,-4,0E«3,0,,0i/ 

C 

C«»«« DEFAULT values OF SOME INPUT DATA, 

data QDlAtL/57,8486/; L0NQL/t75.4736/ 
data tzero/65/t azBro/o;/, vzeXo/oi/ 

DATA ITAPE1,!TAPE2;iTaPE3,ITABG4/ 5, 1, 2, 3/ 

data DTLA/0,01/, DTA1/i;oO/, DMIN/O.OoBZ/, OPAX/ZOtO/ 

DATA EPTtNV/0,0001/, EPBlG/O.OOl/ 
data VOLIM/O;i/, VLIM/40,0/, MLIH/20,0/ 

data ERM/ 1691.0/ 

data IPUMCH / 0/ 
data order/ 14*1/ 
data jspent/0/ 

data WlND/ 26 ;o/, WtJAZ/0,/, ALOM/O.V, AHIGH/IOOOOD ,0/ 

data nlev/12/, ale9/o; , 60,?ipo,' * 266, :soo'.,io8o,,to0o:, 3000., 6000'., 

, 10000,r50()00',)r8’st00050i/ , 

e 

c 
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LAREL HAlN 


127 

126 

129 

130 

131 


HI 

133 


13< 

135 


139 

137 

136 

139 

KO 

HI 

HZ 

H3 

194 

145 


146 

147 

148 
' 149 

130 

13 1 
192 

• 193 
1?4 
135 
196 
157 


138 

199. 

160 

141 

142 

143 

164 

165 

166 
167 
166 


C 

B#*** NAMEUlSyS 

NAHELIST/’DLIST/ SDLATt:,LONGL,TfBRO(:A?EftOf VTEftO,- 1 TAREl , 1TAPE2 * 
i -I TAPE 3^ lTAPE4,nTLAtD TA!,D HtN.DH AXj EPT IN Y , BPB 1 6 , 

I ER|5j9oL1M,VLIM»HL1H,OrHW,IP0NCH 

NAMEUST/BLIST/ fJPL.PUM, AZSDLiNAN6.ANG,WiNh, ' NAZ*ALOW, AHIGH. 

; NLfev;ALTW,3PeEb;BlRrlR0T;iPRlNt.JSPEtJT 

NAMEUIST/FLIST/ WPL;4ZQDL»ELGDL;tll5lbiHNAt,NLEv,ALEV,lROT, IPRlNT 
NAMELIST/ULliiTR WPLi AZQdL >NAMG> A^G^W t ND, ALOH»AH 13U , iRof , IPRI I^T 
c _ 

c ' 

C***« INITTAUZE THE PBOfiRAMi 
PRINT 1000 
SAVEW » WIND 
SAvEA <■ UNAZ 

LEVEL » NLEV , . _ „ 

SAvEL ' ALOM 
SAVEH B AHI5H 
RBAD(ITAP.EliDLtST) 

I s 0 

READ(ITAPEliBl.lST» 

IFMPBlNT.EQiZS 1*3. , ... _ 

REaB( ITAPEIiFLISTI 
lF(lPRJNT;Ea,2> !«3 
REaD< ITAPE ltULlSn 
IF(lPfilNT.EQ?2} Ib8 

IPRlNT ■ I 
PRINT lOQi 
CALL LlSTdTAPElI 
PRINT 10B2 
CALL L1ST(ITAPE2> 

GDLATL » QDLATL*BF{e» 

LONGL » LONI3L*CF«2l 

TAVill 8 TZEBO 

TAV<2) « AZSRO 
TAVI3) * VZERO 

call rdata 

NCON • si 
CALL 5ETIV 

DBCOpE tNAHf, 10031 LEFT,R1BHT 

IPdPONCM.NBJOl PUNCH lOOOj LEFT 

6 


PAGE 
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169 


SUMMARY AT BURS-flUf, APOGEE, AND IMPACT. 

170 


5ROT » 1 

171 


NLEV » 0 ■ 

172 


MLAST » 0 

173 


wind * 0,0 

174 


UNA2 ■ SAVEA 

175 


alow b savel 

176 


AHIGM a SAVEM 

177 


IPRINT » 0 

176 


A2GDL » gUMHY 

179 


READ(ITAPE1i6l1ST) 

leo 

' 

IF( AeS{AZGDU»DUMMy} ,LT’. l,0E-05> Go TO 200 

ifii 


iFtuiND) iiouaoliio 

182 

110 

MLAST » 3 

183 


WNAZ » WNAZ»CFt 2 > 

184 


CALL HCONST 

185 


GO TO l^O 

186 

120 

IFCNLEVJ 200il40*l30 

187 

130 

MLAST « NuEV 

lee 


DO 131 1*1, MLAST 

109 


WAZ = DIR(1>»CP<2> 

190 


WAZ = AHOD(HAZ*PI,TPn 

191 


VELE(I) » SPEED! I )*S!NIWAZJ 

192 


VElNC) * SPEED! I >*ces<WAZ) 

193 

131 

continue 

194 

140 

DO 150 JfliNPL 

195 


PLHUI » PLM(Ji»eFf9J 

196 

150' 

continue 

197 


AZGOL = AZGDL*eF!2> 

198 


DO 160 JsliNANG 

199 


ANGiJi * ANGtJU0F(2l 

200 

160 

continue 

201 


CALL BAIURPT) 

202 


IEND » NSP 

203 


IFUSPENT .EQ, 0) 93 TO 200 

204 


IF<1'PRINT,EQ.-2J ipftINt*! 

205 


call SPENTMBOT) 

206 


IEND * 0 

207 

C 


208 


F!ZJ and DFIZ) curves; 

209 

300 

continue 

210 


IROT * 0 • 


LABEL MAIN 


page 



label main 


PAGE 
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211 

212 

213 

214 

215 

216 

217 

218 

219 

220 
221 

.222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
23a 

239 

240 

241 

242 

243 

244 
2'45 

246 

247 

248 

249 

250 

251 

252 


NLEV * level 
mind ■ SAVEM 
WNA 2 » SAVEA 

IPRlNT « 0 ....... 

A2GDL »■ DuMMS 
REABUTAPElifLISTl 

ir{ABS<A 2 GDUi*D 0 HMY} iLTl 1 , 0 Bb 05 I GO TO 30 O 
MPL » WpC»CPt5J 
A 2 GDL » AZGDL*CFt 2 i 

ELGDL « ELGDL*GFI2i • . .. 

WNAZ « HNA'ZfCFIZ) 

HLA6T • 3 

. CALL F 7 C()RVnR 0 T> 
lEND ■ NSP 

TONCR TILTi UNlr-MlvD EFFECTS, AND CQRlOLlS DEFLBCTlON? 
300 ‘ CONTINUE , ■ 

IROT ■ 0 
wind » s*v'ew. 

ALOW ■ SAVEL 
AMISH ■ SaVEN 

IPRINT *„0 _ . . . . 

A 2 GDU ■ DuHM* ■ 

RBADUTAPEliULlST) . 

if(Abs(a 2 gdl»dumHy} ,l;t. i,oeb 05 » So to Ago 
MPL = Mpt,»cr{5) 

AIGDL ■ AZGDL»CFI2) 

DO 310 J* 1 ,NANS 
■ANG<J> » ANC{jUBFt2) 

310 continue 
MLA 8 T » .3 
CALL EFFBCTIIRBTI 
I END ■ NSP 

C . . 

C»*«« 

400 IFIIEND lEOt 0 ) SO TO 420 
D'O-AIO JfliN*P 
READCTARE 31 B( 3 ) 

REAOUTADES) MAx;{MIi 4 ),I«l,MAX) 

REA8UTARE3) HAjfiixlli'SjilPltHAX!! .. 

410 CONTINUE 

420 RBAD(nARE 3 |ENO«)l 60 ieRRi 45 D} J 
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LABEL MAIN 


PAGE 


233 

450 

PRINT 1020 

234 


GO TO 500 

235 

460 

PRINT 1030 

236 

300 

CONTINUE 

297 

C 


296 

lOOO 

F0RMATCim,48H SENSwSD CALCULATIONS BEGIf^ — 

299 

1001 

FORMAT!///, 21H INPUT BATA SET NO. 1,//J 

240 

1002 

F0RMAT<1H1,21H INPUT DATA SET No, 2,//J 

241 

1003 

FORHATIA61, A195 

242 

1020 

FORMAT!///, ?5H eND-0F;ni.g NOT F.0U)iD ON.JtAPBS.) 

263 

1030 

F0RHATI19HSENS-5D TABLES FDR ,A6i) 

264 

1050 

F0RHaT(///i33H THIS concludes THE CALCULATION,) 

265 


STOP 

266 


END 
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1 

2 

3 

4 

5 

6 
7 
e 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
2B 


SUBB0)JT|NE APCAL 

6 CAUCUtATJON AT THE APOGEE POINT', 

C 

integer RKSTEP 

REAL LATU,LQNGL,LAT,LONG*lPT*fPAZ,IP« 

DOueuE PRECISION DRSiDPC 

COmMON/sToRQI/ G0,eR,0MEGAfGK4RE,R(},ESQ|9SQUE3Ql 
COMMOS/STOR03/ XO , YO IZO , R , RXY , ALT , Sr tHLl H 

COHMON7STOR047 GDLATL,‘LATUi'LDNGLi tMeTaL,PH!L* GOLAT, LAT5 LONGi THETA» 
; PHIT8AMAL;PSIL»GAHA*P9r;ELQDL*AzGDL,ELG0|AzGD,EL, 

. AZ.SRAZ 

COHMON/STOROT/ T;DTiTNPl-DTLA|DTAI»TB0,RKSTEP,EPBI6,EPTINY, 

, DHlNiDMAX*KSTEP»N6oOD 

COHMON/STOR087 S J 14 » • F { l4 * 5 ) »PPS< 14 ) »DPCl 14 » iERW ( 14 ) VsP (6# 20 } 
C0HM0N7ST0R1S/ APT. APAZ , APAi ApR;BOT,aoA,SOR , S0VL*B0EL ,B0'A?, IPT, 

. iPAZiIpR 

0 

RATIO » RE/8 
XS * RArl0*3Cl) 

Y9 ■ RATt0f3(2> 

ZS ■ RATlOftSO) 

call AZRANINSiYSJzSi.TJ 

APT » T 

APA b alt 

APR * SR 

APAZ e SRAZ 

RETURN 

ENO 



ORIGINAL PAGE IS 
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93302 01 11 « 29;75 16;120 


label azran page 


A3 ARI b UAT 

44 Aft2 B LATL 

45 AR3 B UONG » LONSL 

46 OG P DCOS(ASl) 

47 ■ DE » DG#DSI8UB3) . „ , ^ 

46 DF ■ DC0S<ARa)«DSll4(ARl) - DG#DS IN (AB 2 i *DCOS UR3 » 

49 SRaZ * DATASStDEiDP) 

50 1F(BRA2.1T,0'.> sBa2«braz*tbi 

51 RiniBN 

32 END 



Y 


43302 01 11«29«79 lOtlZS 


label 0ai page 


subroutine 0AlU0Ot> ^ .H- 

-CiLGUL-ATtONS FOR BURNiOUTi APOGEE AND IMPACT i 

intEgep BKSTEP _ , . 

R6AL LATL,LflllGLiLAY*LONG;MAS9*HPRoP;HPLrffDOT«M*BfliLC5,iiep;iY?tipT. 

, tPAJi ieP»MSEP 

DOUBLE precision dRb.DPC 
CHARACTER NAHE»60 

External soLV3biSoiv$o 

dimension SAJEHljlSl ; ^ ,, 

COhMON/SToROC/ ZEROiONEiPIiOPI »HP1 *TP^i7o1^(9) 

COMHON/SToRei/' G 0 .GAi 0 HEGA|GHiRS,R(}tGSO|BsniiESQt 


COHMON/SToRO?/ Xe,Vo;ED,RiR)iY,ALT,SRlHUR 

cOhMOS/STDROA/ GDLATt';LATLfLON®L,TMEfAL,PHfL»GOLAT,LAT*LONG, T heta, 
; PHIi'8AMAL;PSIL,GAHA,'P9l,tLQ6L»AZQDL,IL4D|AzGDi'EL, 

’cOHMON/STORp?/ VT<l)VVA<3>rVE\.ACC."vtfUH*vn'^’ 

COHHON/SToROT/ t;dT,TnP1;DTLAiDTA1ItBo,RRSTEP,EPBIG,EPTINY, 

, DMlNiDMAXjKSTEP^NQOOb 

. COMHON/SToROS/ s<l4»;F(l4,5)»DPS(iA)lDPClH4)»ERH<14)V8P<6»20 J 
COMMON/sToRO?/ H*sS,MpRoP,MPLi3UBM{32)7MDof,THft,HACH7tA,AD,QS, 

i . . SL0P87CMP,LCQiLCP';iIii.SYLS2_, 

COmHON/SToRIO/ WE,HN;MLASTjWlND*WNAr;AL0ff7AHT0H,NLEV,‘Ai:Eva0bJ, 

• ALTrtllOO>;vELEClOb)ivELN(*Obl 

COmMON/SToRI?/ NPMAX;(JSMAXyNPL»NANQ,NVAR»NP,NPl, jSPENTiJS, JSC iNCON 
CDhHON/SToRIS/ APtrAPAZ,iPA,APR;BoT,BOA,BOR,0OVL,BOEL,BOAZ,lRT, 

, iPAtilPR. 

COMMON/5TOR14./ I TAPfil.'l TAPEg«lTAP&3. lTAP.e4:iP.R INT, IPUNCH, INQ.EX „ 
COmMON/SYoR??/ NAH?i'#RRAY(800J _ 

common NPST|NSEP;NSR*NH(14|,KSP<20l»RMASS,T3°,TAX<3)7AflO<61IrB(4i» 
, PLH<tO),PHTt52),T8EP<32UMgBp<32)»TSP(20).SPMtCO),A<9t,14) 

EQUIVALENCE tAaRAYtlI,'SAVEtl»l» » 

SAVED ■ 0ME6A . . . . _ . . 

IFdROT.EQ.Ol OMEGArriDEvlO 
DO 700 iPaliNPL 
MPl » PLHnPI 
PRINT 1000 

PRINT 1010, NAME;hPL'‘ODi*ZGDL*CF<i| 

IFIMUAST ,EQii 9> PfilH! IPU. . - . 

1F(HLAST .EQ; 3> PRINT 1012, HIND, HnAz*CF(i) , ALOW, AHTbH 
IFINLEV ,GT, Ol PRINT 1013 
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1 | 



43 

44 

45 

46 

47 
46 
49 

30 

31 

32 
? 3 ._ 
34 
55 

36 

37 
36 
39 
60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 


C 

c 


c 

c 


IF(IROT,EO,0> print 1014 
IF(IB0T,«E,0} PRINT 1015 
DO 600 lAslrNAbG 
ELGDU = ANGltA} 

IFdPRINT.NerOi PRtMT 10 l 6 » ELSDL*Bf<iJ 


BEGIN 5 gD calculations 
NCON » Q 
call SETIV 

IF(HLAST.,GTL0J NC0^(J|2 _ 

call RK 4 «i 4 i 8 lat:vTO’-'— — 

200 CALL PCH<l4rlOLV5Dj 
IF<T-T3D» 2lO»249r245 
210 call APSTEP8S2ao;s326) 

220 call RK 4 I 14 *IoLV 5 D} 

call APsTEPiJ 2 Jflif 2 a 5 l 

225 IFtABS<VT(3)} ,Gt, VLIM) GO TO 220 
GO TO 350 


BEGIN 3-0 CALCOLATtONS 
245 KStEP » 6 
250 call PCM(6iaPLS5P) 

call APSTEP8S23o;sS20} 

820 call RK4I6i8DUV3D) . 

call APSTEP8S33DiSS2$| 

525 IFtABSIVTOn ,gT'. VLIM) GO TO 320 
850 call APCAL 
RKSTEP 9 3 
DT ■ AMAXKDTiDtAiI 
KSTEP « 2 

call RK4(6i80U3D> 

360 call PCHI6i90LV3D) 

call iPsTEPtJ3eo»ss?ei 

390 call RK4(6«80LV3D) 

IF(ALT,LT,HLTM| SO TB 400 
call IPSTEPl5330*S3?0j 
400 call ipcal 


SAvE(lA.ai) R ANSltAInCFtl) 
SAvE(lAi02> e apt , „ 

SAvE(IA(03> a APA 
SAvE(IA|04I a AP8»eF(4> 


label BAI 


PAGE 


2 
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U/kBEL BAl 


PAGE 


65 

66 
87 
68 

89 

90 

91 

92 

93 
'S4' 

95 

96 

97 
96 
99 

lOQ 

101 

102 

103 

lOA 

105 

106 

107 

108 

109 

110 
111 
112 
113, 
llA 

115 

116 

117 

118 

119 

120 
121 
122 
123 
12 A 

125 

126 


SAvEUAjOS} B BOT,- ’ 
savE(IA«06} a boa 
S|j(E(iA,07) A BoiACFtAJ 
.J BQ jL . ' . 

SAvlMA;li9) A B0BL«CF71) 
s’AvE(lA,io) A §OA2«OS21) 
SAvE{t*iil> A tPT 
sAvEnAiJ2> f iP8*eri4» 

SAvStl*il3> A IpAZt^BFU) 


IF(JSPENTlEfi;0 .QR. S&P'.EQtO) 60 T6 609 
DO SfOQ 19 a1iMSP 
read ( 1 TAPES) e^3) 

RIADIITAPE3) HAx;(« < I.’Aji l9l»H'Axr 
riab { 1 TAPES ) R Ax; ( A t r;5 1 • I Ai;uAx.i 

HBltEm APP8) ^sB ttBlJBPMns>.Bl3i 


HPlTEtrTAPM . . 

WRlTE<tTApB»}’HA|,lA<I,5bI»liHAXi 
RRITEUTAPED 
500 CONTIfUIE 
600 CONTINUE. 

rp£ip8iuT.,NiialBa.ifej.i..uai 

PRINT 1020 . . ■ 

00 610 lAitiNANQ. 

PRINT 103Q) lSAVBm:r>iMi||3l 

610 continue; ; 

1P{IPUNCH ;bo. 0‘l 6Q.TO 706 
PUNBIl.10.90i tlA£m;xlU£EiIAal«IAfJUilAML.. 
PUNCH 10?li tSAVBcTAiUliliAliNANO} 

PUNCH 1092,’ tSAVB(iAiS9)»I8*l|NAN6 
PUNCH 1093. tSAVB(IAl09>»i«*liNAN(l.l' 

700 continue 

ONEGA ■ SAVBD 


1000 F0RHAT(iHi,8|K«aSA NAULOPS FLI6HT BeNTBRp/',2»H 1 ?ALL 0Pb 
,'GInIAa//'|50H TRaJESTBRY BUBHABY AT. BUB^i-SUY* AP_0OBI ANfi 

1001 FORBAtaHl) , . . , .. _ 

1010 FORNAmlH llllHiqtE'. i .ABQ./ylBH BAT LoXP * .Tfi 
illH launch A| * aF7i2?AH DI6*/J 

1011 format(13h Hind ■ .ahzeroi 

1012 FORMATriSH IInT “ ' | ^Tf lvirtT/aB ~ ' ' 

■;i8x»5tl) FR0M*F9;ij6M FT tO,F9.1,3H Ff) 


IlLMir 

..WPAfiT, 


m 

/n. 
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127 

126 

129 

130 

131 

132 
l53 
13-1 
l3s 

136 

137 

138 

139 
HO 
1^1 

142 

143 

144 

145 

146 

147 
146 


.somvariablE Wind sup^’liEo by iiSERj 

,18MNOK*ftOTATlN8 Hof)EL> 

,14li(RCITATlRQ MODEL) 


1013 F0RMAT<13H WtND I 

1014 F0RMAT(13h e»RtH • 

1015 F0RMaT(13h B»RTH i 

1016 FORMAT(iHi,5*H6EYaIlBD PRINT.OUT oP trajectory TOirCffuBCH EUHVAtIO 
,N »,F7,2,4H1eS;/V54>^ units ABE FiPiSiSf^REE* EXCEPT flANaE(RHr‘AND 
, Acc<ao»,,/,/l 

1020 FQRMAT«///,tSXJ21HA P 0 fi E E,23yi27HB U" R N 0 
, U TI23IU21HI M P A C T. /, 11X,28 Vjy*) ;nt;yauH»l ,5X> 
,'2J<iH#),/,2Hi2t!EL,8Xi4HTtME,8X.3HALTj7X,5HSAS(je71Xi4HnHB,0X»3HALT 
;.7x<5«RAK6BtlX73Hvl!L75X.6HrLr<^EL;3R;«KPLY/Azi7)i,«HT]RB;iX*yH9*NQBi 
,'5K,2HAZ77»iRJ5KTBE8>i6XVyflfSE^l7Xt-»HtFT5,'6X*4H<NHi,7Xi5K(SECr.7Xr 
,4H{FT>l6X,4H^NH)i4Xr«H<FT/SECi»4x,5H^DeG),4X*Sfi^DEa)S7j(,SH^SEC), 
;5X,AH(NM>,5Xi5l!l{DE8)i//) 

1030 FORi1AT(F8,2»fl2.2,fl2;o.F9,l,El2.8iFl2'.0»F9,liPir,27P972*PlO',*. 
,Fl2V2AP9;i;Pt.3> ■ 

1090 FORMATiiN. launch BlBVAtIONS tOEQ) FOR [HpACT lYARSE, S-0 FLTafiL AN 

,'D AP06EE’AiTI’i?V!3|Tf 4X*P872TIT ’ 

1091 FORfiAtdtMPlET RAnIbB !N NM » 7/ i < F7 1 278 (1)1 i F7i 2 }) T 

1092 FORMAT ( »Bb 6' PLt/BL IN DEl ' l/ifFT ,8l«<lXiP7.2J ) } ' 

1093 FORMAT(»ApoBiE ALTITUDE IN FYT,/, (Pf ,o;6UX|P970 I ) ) 

RBTURN 

END . . . - _ " 
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19 

2P 
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22 
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24 
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32 
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40 
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LABiL BOCaL PAGE 


e 

c 


e 


SUfiftOUTlNE SDCAL .. 

CtL§Ut:ATlON AT SUPhtOUT POINT, 


lNTtG6R «KSf|P 

RIAL ipfrjpn.tPB,^ 
DOuetE PREClIlON P^S 
COMHON/SfoRX/ z1r6; 

COHMON/STqROI/ QOl'iA 
COhMON/'STORJS/ X0,V|! 
COHMOH/'STOMii.liilSl 
COMMON/STOR#«^ BTlt? 
COHMON/SToROX/ TipTI 
■’ DHjfii 
COMMON/5TOR9I/ 5U4| 
CBmMON/STORIO/ Hg,wa 

il ■ 


tOPg»ARl»AB2 

Op^PlfQPI.HMlTPjK 

*0HEOArGHjRB»Ra,ESOiSsHl*ESOl 
aOtR.RXYiALhSRjHLiH': 

:>/A(8>rVBL tAee;voLiMi.vi:iH 

» .. . . ._ 

TNPlfDT{.AiPTAli.TBO,RI<STiPtEPBIS,EPTJNY, 

CMAXiKmejtNQOOo . 

if(l4i5i;pPSil4i;pPQYlJl»Emi4r/BPW.|20). 

XMLASTjH]tfDiWi^AlJALOLAJlJ55iNtEyS'AL'fV(100>, 


COHHON/STOR15/ APTJA 

i I? At I 


1bOVL»BOEL,BOAZiTPT; 


Sip5{4»+OME9A*B{8) 

SSeS(9)<>0HEGA»S(i) 

VT(4) • HTlU)i£l. 

VT(2) » RT!t4>«Sl If 9fU»!»S2 ♦ Blt(6)«si6) 
VT(3) ■ RTIl7)«St * iTI(«)»S,2 ♦ Rt 1<?)«SUJ 
BDT “ T 
BOA s ALT 

SI » VT{t)*»3 * 9TtBj|««2 

BOVL » S9RTIU ♦ Vl(31***L _ 

ARl*-VT(3) 

AR2«SQRT{S1J 

BOEL * DATANtARl/ARn 

ARlPVT(2} 

AR 2 iiVT(l) 

B0A2pPATAtJ2l4Rl,AR2i 

IFcBOAZ.UT.OS') B 0 At»BOAZ*TP| 

Si « RE/R 

call AZRAN(81«sa)lSlt8r*)#St*B(3ii:ff 

BOR » SR 

RiTORN 

END . . _ 
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label’ effect 


1 

2 

3 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 


subroutine erFBCT(tROT) 

C UNIT-MINO EFFECTS, CORIOLIS DEFln AMP RANGE DERIVATIVE; 

e 

REAL LATL,LONGL,LATiLONG»HASS,HPRoP;MPL,HDbT»HACH»LC(l,LCP,IYX,Ipt. 
: IPAI.IpRiMsEP 

character NAMEmBO 

DIhENSION SAVE(61,9I 

COMMON/STOROC7 ZERDi ONE, PI i QP I , HPUTR I ,CF( 9 ) • 

COhMON/SToROI/ G0,0f4iQMEGAfGH.RG.Ra,ESO,ESbl,ESOl 
COMHOR/StORO«/ GDLi«rL,LATLjLO&QL,THETALi^iL*GDL«T,LAT, LONG, theta, 
; PHI78AHAL.P5IL,GAMA;PSi,ELGbLazQW,EU6D,A2GD,EL, 

V AZ,SR*Z . , 

COHMON/STORP9/ MASS gHpRoP , BPl, , SUQM ^32 J ,MC0 T, ThR, HACHTCA , AD, 03 , 
SLOPe;cHp;LOG,LDP, l?YiSV,SZ 

COliMON/STORlO/ HBfHR/HLASTrHIND.MNAZULOU, AH1GH,NLEV« ALEVdOO I . 

Aj,T«[|OP)*yBLitipOi±Y6LNllOOj _ .. 

COmMON/'SToRIS/ NP«AX;NSnAXrNPL»NAN5jNVAR*NP,RPl,jaPEljT;jS,JSB,NC0N 
COMHON/SToRtS/ AEtVAPAZ, APA, A gR;BoT,BdA,Bofi,BoVL;BCEEiBOAZ, 1(?Y, 

'c0M«ON/STORl»/ i?A^ii;if APB2»lTAPE9rlTAPe4;iPRINT,lP0NCH, INDEX 
C0n«0fl/STOR?9/' ABRAYIROO) 

common NPST,NSEP;NSP:NMa4j,KSP<20i;RMASS,T3D,T'AVl!3>7ifjQI61)^B(4»i 
PLMt20>,'pHtt52I,TSEF{82i;Mslp<32>;TSPl2d)iSPHt80i,'At9r,i4J ' 
EQUIVALENCE tABRAYtir,'SAVECi;iH 
C 

PRINT 1000 

PRINT IplO, HAHE;MRLMS0,AZ8DL«CF<l|,WifJD»WIND»CF<6) 

IR<1ROT,EO,OJ PRINY IQli 

IF(1R0T,NE,P} PRlNt 1012 ’ 

saved « OMEGA 
del * -0,5«CF<2} 

CFC = CFI2>»CFM> 

C 

DO too JAljl^ANfi , .. 

ELGBLT ANgUJ ' 

SAVE(J»1) b EL0DL*CF<D 
C 

c no-vind impact and CORIOLIS deflecYion -4.-..- 

IF{IPRINT,NE{OJ PRTRT lOBOr SAVE<j;i> 

ONEGA ■ l.opclO _ 

NCON ■ 0 ■ ‘ ' 

call 8BTIV 


?AQE 1 
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call »PL*S« . 

UR * JPR , _ __ 

PN * iPAi 

J p ( iPm wT.WjlQ } iRtW T lozir SAVBlJllI 

OHEGA » SAVSS" 

CALL SBftV 
call 8PL*SH 

SAvBtJi'4) * ilPR«cSjt|PAl^UUN«B0S(?Ny>*Cfj4} 
SAVEtJi?) • HPR*s!8UPA25-UNiSlN(PNT ji'CPlij 
I F { I RQTjiiiai_4G.i£itt I : 

OHEQA « tlOB-lO . _ 

QO TO 2Q - ' 

iO UN • tPR , 

PN ■ IPA? ... 

20 SAvE(U, 2 » UNf5F(4) 

SAVE(«ijLaJ..» PN«CEXll 




HSAD»HINB RANGBi 

IF(IP»IN^NBJ05 9R!Nt ipfji . 

NCoN •• 2 
call JETlV 

H8a2 I.-E5J - . 

call mcqnst 

CALL IPLASH . 

TIRM » lPR*l*IS(IPill>t8lN<iyj)4C03<rpAZJ*CD3{PN) ) • 

SIGN ■ TBRK^ABSjTeIb^ 

TIRM I S0RTltPB«ipR*UN*UN-2,0*tPR«UNAC0SUPAZ-PN)) 

HIAD..I, aiaMtlfiBHAXJlIMlL 

SAvE(Ji7J » NEAD>*crr4J 


tail-nino rangb. 

irnpRiNT.Nejo} eptsT io>4V sAV'njiii ' . ' 
call SETtV . . . ..1 

WNa2 P ESj. t.Pl ! 

lF(UNAZ,HE,TPn WNAI*WNaZ-TPL . 

call HCOHST .. 

call 9PLISH 

TERM • lPR»Hlnt)(|PXZit8iN<gN>*068ttPAzT *C08<PN yr 

tIrIH J SQRTiip^llpl*iN*UN-8 ,0»tRfi»UMCOSTrPA7«P 
tail * (9fotl'9feRHWPfU0 

SAVE<J|8J « TAlL*cr{4) 
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p o 


85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
lOD 
101 
Ida 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C CRoSS-HIKd SAN6E, — o»-»- 

IFtlPRlNT.NBiO) Pr!nT loasy 9AV6<j;i) ■ 

CAUL 9ETIV 

WNAt * p8i * HP! ■ ' ~ 

IF(WNAZ,QE,TPn WNXZiWNAl-TPI 
CAUL WCONST 
CALL SPLASH 

SIGN 9 1,0 - 

TiRM ■ SORTtlPR*rpi*UN»UN-2,0«lPN*t)N«COSfn5'APipFJn 
cSoBs B (sisnbte&mI^ITnd ■" ‘ 

SAVE<J,9J • CRDSS*er<4> 

c .... 

C TOwER' TILT, 

IF(IPR1NT.NB',05 print 1022i 3AVB<J»ij 

eLGDL ■> ANQtJ) ♦ del. 

NSON » 0 
CALL SETIV 

call splash 

sign ■ (tPR-UN)/ABS(IPR-UN) 

TERM • SQRmPB*IpR*U)il.uN-i,0«lPl«*UN«c6SVlPAXiPNh 

TTILT • {SlQN*TE6Hi^DEL . ' 

SAVE«J,6J ■ TTILTbCFB 
100 CONTINUE 

e 

iF(ipRiNT;Ne;o} print togs 
PRINT 1030 
DO 110 gfltNANS 

PRINT 1090, esAvEja*Hii«ir9j ' 

110 CONTINUE 

IFdPUNCH .es, 01 GQ TO 120 

PUNCH 109o, NA(gG.(5A0E«j;i}»J»i;NANeI 

PUNCH 1091, {SAVE<a#7l, J*lrNAHe> 

PUNCH 1Q92, (SAVeiJilhJillNANQI 
PUNCH 1093, {SAVe(Ui»|iJ»l<NAN8> 

PUNCH 1094, tSAVe(a,S),J*lrNAl!tOl 
PUNCH 1095, tSAve(U,4»,J*i;NA(l6| 

PUNCH 1096, tSAVfi(J#9|,J*lj:NANG> 

120 PRINT 1040 

CFC " Cne}^BFt6J ■ “ ■ 

DO 130 J»1,NAN(S 
SAVE(J<2) < lATEtJi'ZUCFie} 


4 
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127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 
136 
139 
1^0 

141 

142 

143 

144 

145 

146 

147 

148 

149 
l5n 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 
167 
166 


SAve(iJi4) ■ 
SAvC(J«5) ■ 
SAvE(J,6i ■ 
SAvE(Jt7i > 
SAvE(Ji8{ ■ 
SAvE(d»9> > 
PHINT lOSo, 
130 continue 

ONElitA i saved 


0AVEtu«4)«cria) 

SAVE(j;f>ACFt8] 

XAVE(j«et«cn8) 

SA6E»j;7UCFC 

«a?eij:b>»cfc " 

lAVEfJiflftCFC 

tSAve(U,l),I*U9> 


1000 F0RKAt(lHl,2»HiilASA WALLOPS FLIGHT bBfTfERf/TzyH WSlLOPS' t SLA' nD, VIR 
,'GlNlA*7/t60b UNiT«WINd EfTECTS, CORiDlIS DEFl-ECTtON ANb RaNsG DEBI 
I'VATlVE,//) 

1010 fbrmaT<i3h vehicle I ,A60,/,18h PaV loAd * rrfi2,'4H 

• 13K launch » ;f7,2|4H DBG*/,13H hInD ■ |P7;2;9R Pt/sSC 

•,'F7.«,9H H/seO.) 

1011 F0BMAT(13h BARTH | , l6l4NOtg-ROTAtlNQ "riODEtY ' ' " "" 

1012 F0PMAT<13h earth .| ,l4MR0TAT)N(i HoDELj 

1020 FORMATIiHi.SAHDETaIlBD TRAUECTORV roR No mind case with ELEVATI 
,*ON *,F7,2,4bl DBG;/,54H UNITS ARE F3P4s.DE0rEE, EXCEPT RANGElNH) AN 
;D ACC<G0)../;29H earth 18 NDNi.R0?AT1RQ MDDGLm//) 

1021 F0RMaT(iHi,58HDETaHBD TRAaEcTOR r FqR No rind case with ELBVATI 
,0N ■,P^,2,4H DEGi/,S4H UNITS 'aR 6 F«P«S-BHQ'REf ,'BRBE pT FFyNlSfTNMt’ aR 
,D aCC(G0>,,/,29H EAaTH IS ROTATINB H0D6L.,//I 

1022 F0RMAT{iHi,36HDETaIiED TRASECTOR? PoR RANGE D£RV CASE WITH ELBvaTI 
.‘ON P,F7.e,4M DEG;/.?4 h units are r“P-S-DBQftEG, EXCEPT RANQEINHi AN 
i‘D Aec<Q0J.i/O 

1023 F0RMATtlHi,5»HDETAlLEQ TRAUECTgR7_PoR HEAD «IND CASE HlTH ELBVATI 
,’ON p,F7,2,4h DEG;/;J4h UNITS ARE'rSl^-DEQREf; BXCEpT SJHfiETTIHi AN 
,D ACC<00)..<'/» 

1024 F0RMAT(lHl,?6HDETAtLED tRAUECTORV PqR TAIL’ HIND CASE WITH ELEVATI 
,ON «,F7,2,4H DEG,/,S4h UNITS A«E FWpffS-DBofiEE, EXCEPT RANGEINmI AN 
I'D ACC(G0I.,//} 

1026 F0RHATIJHi, 56HDETAILED TRAaECTORY F oR CROSS WIND CASE WITH ELBVAtl 
,'ON p,F 7,2,4H DE6;/,»4H UNITS' A^B P*P-S-DeoffEtr, EH'CEPT WgETNR} AN' 
I'D ACC(60),i//) 

1026 F0RMAT41M1) 

1030 F0RMAT(/7,13H PiP-S SVSTEM»7/, 

,‘1X,9HLAUNCH EL» 7Jl , 15HN0-HI ND lMPABT*7jl,2oflC5RlOLlS DEFLECTION, 
i5*,tl«RANGE DERlV,16?;2eHU»lIT « HiND EFFECTS,/, 

,3X,5H<DEG),7lt,9H8Afi8BlNMT,3!if;?ffA2tBBSl.5«,9HNairrfnNHrr3nf7 

,aHEAST(NH},6X,SHtNH/DEO);lOX;i2HHB«b<NM/PFs),5X,12HfAtL<RH/FR8I, 
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label SFfECT f-AGE 


1«9 

170 

171 

172 

173 

174 

175 


176 

177 
176 
179 
1«0 
1«1 


1«2 

1«3 

114 

185 

186 
187 


,5X,13HCR0SSINM/Fes1i//I 
1040 F0RMAT(/7.13M M-R-S sysiBHj//, 

,’1*,9KLAUNCH BLl7!<j 1JHNO-.HIND IHRACTJ7 Xi 2DHC0ri0LIS . DEFLECTION, 

,-5X,ll0RANGE DEJJliA8rZ28HtJNIT = WIND EFFECTS, / i . 

,3X,5H<DEQ> »n#9HBAN9et«M)»3X;7HA2tDEG)*5X,$HNQRTH<KH)4 3X, 
.8«EAST{KM),6*I.8H<KM(^DBS);iDX*i2HHeAb<KH/Hpa)f5X,12HTAlL<KM/HI?S) , 
,‘5X,13BCR0SStKM4MPs}i/>> 


1D50 

1090, 

109l' 

1092 

1093 

1094 

1095 

1096 


FORMAT {F8.2iC14;S#Plil2iF14,2*Fli:2;Fl9,T,igx;P10;5;7x<Fi0.5i'7x» 

P 4 0 ^^5 H 

FORMATMN. LAUNCH tLEUATlONS IDaQJ FOR Ul^iT-WlND EFFECTS, RAN8E OE 
Riv AND DORIDLIS pTM-|/1T?7T9U „ TTT: , ,,,, 

FORMATCHEAD UBlT-filND EFFECT (NH/f T/SEC ) ' ./* (F7,3, BtlX|F7, 3 ( ) 5 
F0RMAT<«TA1L UNlTl\i(iND EPFECT «NK/Ff 7SEC 1 1’.7» !F7t3, BUX|F7, 3 ) H 
FORMAT<tCROg» DNJT-wiHD BFREPT <NH7FtZSEC lF7, 5, 8<lXiF7.3 j h 
FORMATORAnGB DE810 tUM/DEfi)',/i.<p7,3.8{iX.FTu3LU .. 

F BRhA T I ' 50R 151 Is. DBFLECTtOH TP NORTH tNM > i . 7 * rF7.3, 8tlX i F7 , 3i IJ 

FDRMAT< 'CORIOLIS DiFLECTlOW TO EAST INM)’,/, »F7i5| 8(lX,F7t3 )) i 
RETURN 

end 
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LA0gL FZCURV 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


13 

14 

15 

16 
17 
16 
19 
So 
21 
22 
S3 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

n 

42 


C 

e 


c 


c 

e 


SUBROUTINE FlCUROOROT) 

CALOUIATION or WlNb*WElGHTINQ FACTOR AND BALLISTIC WIND FACTORV 

RBAU LATL>UCi&Oi:,UAY,t.0NG»MAS3<NPR0g,HPL,'H'DOT?HAGH,LC6,LCP, lYt» !Pf , 

. IPAMfR 
CMaBACTER NAME«69 

dimension RGtlOOIf'FZdOOlf DfZllOt) 

COMMON/STOR 90 / ZERbyONr,Pl 7 QPl»HPl»Tl?Ii'CF( 9 ) 

COMMOI'li'STOROl/ 69#0R*0MEQArGMiREiRaVESO,SSOI »BS'q1 
C6Hl40N/'ST0Re4/ GDLATL;LATLylOBGL,TNETAL.PH!i.»QDi:AT, lAT^LONG, TNETA~» 

PH ii 8 A H A L ; P S"IL » G A H aTp S iTeI GDL i A Z Q D I TeIgD ,1 ZG D r'E L , 

'cOMMON/SToR 09/ H*S^?HPR0P,MPL,9UBHJ32J,MBof,'THR,HACH',"CA,AD,QS, 

V SL0B6iCHP;LBG,UCPVttYysv,9Z 

COMMON/StoRlO/ We,lflN',MLA9T7WdiD*WNAZlAL0H, aH1qH,NLEvULEV(100). 

'cOMMON/StORfS/ NttxlNl®N?LTNfe^VA^fcr,'J5T>Elif;^ ’ 

C0MM0N/ST0R13/ A^ti’ A^AZ, APA, APR,BoT;BOA,eOR,®OVLiBOEC, DpAZ, IPT, 

'cOMMON/SToRiV i?A^EtVltAP«2;iTAP63yjTAPB4;iPRlNT, tPONCH, iNpEX 
C0MMCN/STOR99/ NAMBfARRATISOO j 

EOUIVALENQE tAftRAYtli;ftGl lI) »tARRAY(101);rZ(t)). 

, {ApRAYUen.DFZttrr 


PRINT lOlO, NAME;M6L»Q0.ELflDL*CF{li;AZfiDL*CFtl),WlND;WlNI)*Cn6) , 

. wnaz*cf<i» 

IFnR0T,BOi!jj BJINT loll 
iFdROt.NE.Oi PRInT 1012 
SAVEO 8 omega 

iF(lROT ,E0, OJ BM88Api;OE-10 

initiauze to Uo wind; 5 

lF(lPR|NTjNfi{0} eplNT 10^ 

NCQN » 0 
CALL SETIV 

call splash 

UN i tPR 
Pfi A IPAZ 

RQtil 8 UN 

Fl(i> a a, . . 

DFZU) * 0, 



Page l 



93302 bV iii25o75 i8';i29 


UPEU FZCURV 


PAGE 


-13 e 

44 C INITIAUIXE FCR MAXtflUH W»ND ALTITuUEi 


45 


48 

49 

50 

51 
92 
S3 

94 

95 

96 
57 
98 
39 
60 
61 
.82 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
60 
81 
82 

83 

84 


e 

c 


c 


c 


loaii auoh.ahiW 


ALOW « Oi 

wgoN *2 • 

call setiv 

call MCoN$T .... 

call splash ..., . 

DftOH * SaRTH PB»!Pili«UN»U N>2.0 » IPR»0w»COS\l P AX«PN>> 

RG(NLBV) 8 IPR 

FlcNLEV) «l,e _ 


wind UPTQ other !N.tBBHEblAtE ALTiTQiiSi 
NBnB ■ NLEV • i 

DO iOOliZaJaiKD 

AMIQH « ALEV^N} 

IFcJPHINT.NE'.O) ERtNt lozir ALOW.aHiGh 

call setiv 

call NCORST ... 

call splash 

AHUM • SORT>TPtt«lPB*UN« UN~2.0* IPR>UN»COSU PAZ-PN > > _ 

R&(N> » JPR 
FI(N> » ANUM/DKIoM 
DFZINI « FZltO) - FX<W5l) 
loo CONTINUE 

DFziNLEV) » FZiNLES?) i FZ(NEND) 


IFdPRINTiNEtO} BRTffT 10H5 

PRINT 1030 

ANUH • CP(4)WCF<8) 

DO 110 N*li NLE9 

DNOM a ALEVtW)aCF(B) , . ^ ' 

PRINT 1040., N. AL6.ViH I ;pnOMtRG tN 1 aHF ( 4 1 . R6 ( N > » ANUH , N Vl PEZ t N I . 

ALEvt Nil DNOM 

110 continue 

IFdPUNCH ,SD. Ol BQ TO 12P 
punch 1090, NLevJ{ALBV<N)*W»tjNLByJ 
punch 1091, tr2«N)TNil,NLEV| 

PUNCH 1092i {DfZ<Nl,Nil,NiJiy). __ _ 

120 omega » saved 
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lOOO rORMAtdMlitVHflASA' NAPLoRS FLIDIHY SeNTBRr/:2fH UALLOI*t3 iSU^NDt VlR 
VeiNtA.»//f38ri f;5zi."GulVB_ANp iUlBi) F *eT O<fSt/yi _ 


leio FORMA'TtiiH'llH'ICLI ■ ■■^'A' 80 i 7 dBH“ 8 Al:'TbrD“"‘l Vis,/r 

viiH LAuSe'H K i SFiijW di6;??>sh laub'ch « • .prT2«M D 
■■ TiIh'hTns ;""T"iF 7 { 5 }rH rtvisBD.zH u/ejiiFH M/sseii.^JVzii^a 


,Z FiOM NbRfril. 

1011 FOrMATUBhIIIR 

1012 FSRfiAtrUH' _ 

1080 FORMA T<igi,'|XH|ETA!i 
'.V»54H UNtTSJRB'gi^t 




i .TffilKbB-KciTlTlNl'WDSECl' ' 

I ;i 4 fRiJT*nWH 0 !lir)~'^ 

-- maEcfoi* roR'iF(TnrrTtr-Ns 'flrso' jn* hheri.’' 
^Bigg tr R ANiii»pri~THD'' A Bfc { aa vr r? ?•) " 


Ifffii FORMAT dMl»»5HBETAt£i6'jRAaE&TflPtf_(!0" ^‘*> MOH*Fi9,'0, 

‘,’4111 to'iF iQ;i g/»Mri~ UBTTS' A'gE P»M«9laiiQRIB7 6BCgPT liABaE<BH) aBd'ac 
Ve<G 8 jlf^/i... .... 

1025 FORHAKlMl) _____ 

1090 FOBMA»(^/,fr,lHf^lfTiHwTNP AL^5|«r«ffn^PXB|_JI^ 

iH'BPiziiiAKBiwTsir s]!T;/;ioi^~9H<fTiT7ir7^fftin^rTiii(Bffi rBx?4MiKM i i 


7.lH 'BPizTil4g|BrwW 

i? 4 o FORMATt i 47 r|i « ■ ■ 

10$o FOR(i)Af<<R» IPITtUDS USCOoTFif ^ "iD i^L7A_EH>_0USvfi8«^/» 

. ai/ji 9 <i>ifF 9 . 8 H)“ ~ ' rr. 

169.1 FBRBATOPfZt C0R0eT i??T^Bi0iTHlB« ^.4n j , ■ .. . 

. . FB HMATtingLTII P(tl BaOVg»i2« tgBw4>T(lX*F<.4 >T> 


RITURI9 
EBD 


LABEL FZCURV RaGE . ,3 
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LASBL IrCaIL page 


1 

2 

3 

4 

5 

6 
7 

e 

9 

SO 

11 

12 

13 

14 

15 

16 
17 
le 

19 

20 
21 

„ 22 


c 

c 


e 


SU08OUT1NE IP5AL . . . 
CAlGUCAtIoN at jHPAfT POINT, 


JPlffl|R_iKSI|P 

RBAU lPT»lPA?»tPB . , 

DGUGLl.enBCt«l0/g MSID.pjl „ 

COMMON/SToRP 8/ XG7?9f?0.R*BX«iAl,T^3HiHLlfi 

COMMOS/STOR94/ GDLATt;yLATLTLOtiaL;THeTAL. P~H! L»QDLAt AAf jLONg.YWBfAi 

V PHi;SAHAL«P3lLjGAHA«PGt,ELGJL*AZQBU*EUtt£f7AzGD«eLi 

V Aitikii 

COKMON/STOR9I/ TaT.nNPl«PTl.AiEllAliTaoiRBS!CE£»^BlGiePTl[!lV, 

; DMlEj.iDMAX.K.?TiBfS.9D5D , . 

C0K«QN/'S!oB8i/ SU)<Hr<lA«5>iDlS(l4)i]iP!5iUlIiERWU4j?8Pli.20L 
COmBOS/’STORIS/ AETjsBAZ,APS,A£l!l.;afl.T;eOAt8Qfi.^fc>Q\tt.aOELiBOAZilfiT, 

. iCAEilBR 


call AZRAN(Sa}i9(2i;$(3)»Ti 
IPX » T 
IRR » SR 
IPA2 *> 5RAZ 

return 

ESO 




LABEL IPSTBP PA8E 
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1 

S 

3 

4 

5 

6 
7 

e 

9 

10 

11 

12 

13 

14 
19 
16 
17 
16 

19 

20 
ai 
22 


SUBBOUTINE lPSTEC!(«r») 

e ADJUSTMENT FBR StEfe SIZE AS IMPACT CLOSES, 

INTfGER NkSTBP __ _ 

cOwMOll/StoRfl/' Gf;SSi'6MEiWi'SR71TI7Rir»TsoTBs'5T»6s5I‘" 
COMMON/SToROS/ X(),?0;ZOV«»fiXTVAL^’SRiHLlR ’ 

C0MMON/STORO6/ VY< 3 > , VA (BlyVELf ACSrvflLIMfVU 
CDhMON/SToRG?/ T;DT;TMPi:DTLA7OTAllT0o7Rf?ST6P|EpBlG,BpYlNY< 
; DHlSi»DHAX*KSTBEi'N(SO0D 

B _ ‘ ' 

G8’ R 9072;0 

TT « a,b*DT 

DIST ? {VeL402»TT>*TT 

IFtALT ,LTi bIST) KSTEPrIOOOO 

TT » DT*2,0 

DIST • ( VEL*fl2»TTj,»lTr _ __ 

iRiDlSt ILt, ALt] fiSTUR'N '1 

RKSTES ■ 1 

TT R <LT/(VBU*0O«bT) 

DT R AMINKTT.DHIX} 

RETURN 2 

end _ 
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1 

2 

3 

A 

5 

6 
7 

e 

9 

IQ 

11 

12 

13 

lA 

15 

16 

17 

18 
19 

30 

31 
22 
33 

24 

25 

26 
27 
2B 

29 

30 

31 
82 

33 

34 

35 

36 

37 
30 

39 

40 

41 

42 
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SUoaoUTINE L.TNT(F,XiiFl,Xl.0RDE8iHAX»KB,Ke,kKB,KKE> 

G LInGaR INTEgeOLATlON 9V AITKBN ITEl^ATlON HETWOD, 

S 

INTEGER ORDEB . 

dimension FU>;XU5.F!(i>iXUli 
dimension FFOETtSltll) 

C 

lORD • 0RDE8 

IF(MAX-ORDEB) 10ilB*20 . .. 

10 N 8 ORDER •-FAA-. .... 

order » 0R08R-8-1 

20 K1 K KB 

DO too KK»KKP*KKB 
XNEW « XI (KKl 
DO 30 KsKliKE 

KSAVE » K ..... 

IF(XCKI,6T,XNEI)0 go To 40 . . . 

30 CONTINUE 

40 IZERO ■ K 5 A 9 I - (0ftDeR*lJ/2 

IFdZiRO ,LTJ 1) itgRO*! 

IF(< IEERO+OBDER) ,0T; MAX) IZER0«MAX-0RDBR 

DO 50 JJ'liORDgR . ... 

N s JJ* iZEgO 

FFDETtltJJ) F tF<l2BR0)#<XtN)-,XNEHl-F(N)«(X(lZER0>-XNBH))/ 
';<X(N)-X(1ZEB0) ) 

50 CONTINUE 

IRORDER.EQ.l) GO TO 70 

DO 65 I»2, ORDER 

U * I • 1 
L e II ♦ IZESO 
DO 60 Hft, ORDER 

J » M * IZEBO ^ ^ - 

FFDET(t»M) ■ (rFOETniiIl)*(XlJ>-XNEHJ b FFDEtUIi M )* (XIL I-XBiiB) 1/ 

;<xtj)*xa>) ... ... 

60 continue 
65 CONTINUE 

70 FKKKI ■ FFDBT^oBDERIORDSR) 

KI « KSAVE o 1 

100 continue 

order » lORD 

return 

end 
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1 SUBBDUTINE LlSrdUF.G) . . 

2 C PRINTS DAT* ITBRED 5N FILE NPq,’ 

3 e 


4 CWABACTEg.y4«E«8S_., 

5 ceHMON/STORft/ NAMBjARRATtaOO} 

6 e 

7 REWIND NFC 

8 100 READ<NFCfl 09 Qf G nDkU^I NAME 

9 PRINT lOOg, NAME 

. . ID- GO TO 405 

41 UO REWIND NFC 

12 B 

13 1000 rORMATlASO) 

14 RETURN 

15 END 



99869 Cl 0t-2i-76 


LABEL PCM 


1 

3 

i— 

4 
$ 

« 

? 

9 

.. ^.. 

ig 

u 

13 
1 ^ 

14 

1^-- 

14 

U 

18 

19 

23 
2V 
2? 

24 
24 
29 
26 
23 
28 
2 ? 

43 

41 

43 

44 
44 
35 
34 
37 
34 
39 - 
4q 
41 
43 


■ C. 

c. 


- -t90 

c 

c 


110 


SUBROUT5ME -JieM<NF,S0LV> 

MODIFIED ADACSlBASHFORTH PRE&ICTM-CflRRECTOR METHOD, 

INTEGER RKSIftP 

DOUBLE eRECiSTON DPS+DPC#DRF- . --- 

DIMENSION e<14j,Ca4),PC<l4) . , ' 

COMMON/SIOR#*/ T*DT,TNPl*DTi,A*DI.Al,TBO.RKSTEP,’EPBia,ERTlNY, 

, DMIN-,DMAX,KST-£eiN4IOOD 

C0MH^G4Uft*A7 _S { lA4-^4^U^5)^4>RS44.AJL;DECU41-;-eaH.UA4 2JLL 

C0MM0N/SI0R18/ NPMAX,NSMAX«NPU*NANQ#tII/AR,NP.NPt-< JSP£NI,JSj JSC.NCDN 
COMMON/SIORI*/ lTAPEl,lTAPE24nAPE34JlAPE4»IPRINT, IP.UNCH, INDEX. 
PARAMETER 0955,0/24,0# C3-«‘59,0/24,Df C3*37;o/2i,.0... C4»9, 0/24,0, 

. C5919. 0/24,0, C6i5.0324.-0# C7»l, 0/24,0+ Caa251.0/27o,Oj. 

. C9919.0/270;o 

PREDICTOR ARB ITS MODtF-IER, 

DO 100 t»l,RP 

pPF s Cl,*F(T,4>,C2#iF.<U344CJ»FlI»2),|i4ftFIl,i) 

P<1) = DRS<|J * 0T*DPF 

S(l) = 3U) = C8<*DPC(I> ■ ■' 

GONTINUC.. ' 

FJNCTIO'4 F(a,T+DT) , 

T ! T*OT ■ 

CALL SOLi/ISI ■ ■■ 

CORRECTOR, FINAL WLUE-, AND- JAaILIMUM. fiRROa . _ 

EMAX s 3, ’ • 

DO 110 I«1,NF 

OPF = C4aF( J»5>*CS*F{l,41.C6#FU#3>+fl7(»Fn,2) 

C(I> » 1RS(|> ♦ DT»DPF 
PC< I) » R(II e C<n 

sm s 044-4--*- 09*P«-I4 - .... ' 

6MAX = lyiAX * ERW< I)*ABS<PC< m ~ 

CONTINUE • ■ 

EMaX = cifAXANF 

STOP Sizi-.AJUUSTMENT ACCORDING TO ACCURACY REQUIREMENTS, 

IFU T*&^4lrf«#»W4- -tSE . T NPA-1- GQ-T-Q-404 — 

IFtEMAX vQTi ePBla) QO TO 2-00 

IF<NGOOO ,UI> KSTEPJ go to 500 • - 

IFIEHAX tLT* EPTINV) GO TO 300 -- 


Page 
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45 
44 

43 

46 
4Z 

44 
4? 
50 
54. 
5? 
54 

54 

55 

56 
52 

58 

59 
6'J 
61 
62 
64 

64 

65 

66 
61 
63 

69 

70 

71 

72 

73 

74 

75 

76 

77 
74 
79 
4U 

ai 

42 

84 

84 


GO TO 530 

n 

w " 

• C . DECR£4of= TH5 31EP SIZE, 

200 TT 5 Dr/3,0 

IftTI ,uT. 33IM) GO TO 500 
T 3 T-OT 
OT 8 tt 

1F< (TtTTa4,3),GE,TMPl) DTs C T'tPl-T ) /470 

HO 210 - . . 

S( I) = n»GU> 

210 CONTINUE 

call RK4(NFI90LV> 

RETURN 

2 

.. c. increase. tag Step- SIZE, 

303 TT = DT«2.0 

IF(TT ;gI. 9BAX1 GO TO 500 
IF(<T*TTJl4.3r,G6.TNPl> GO To 5fl0 
DT » TT 
DO 310 1»1,)4H 
QPSUl a Sill 
310 CONTINUE 

call rk4(nf;3olv) 

NGOQD a 4 

RETURN • , 

C 

...C.. RaAS£..CM4NG53,. 

400 TT 3 DT 

DO 410 H IJNF 
opsm 3 s(i) 

413 CONTINUE 
HKST6P a X 

DT *-TN'a4iK4*O-l,OSs07J e T 

CAuL RK41NFI3olV> 

T 3 TNPl 

IF(A4S( T-T8a> ;lT. TBO»ltOe-Q6) CALL BOCAL 
IF(IPRI'JT ,4a, 0) CALL TPRINT(5,1) 

RKSTEP 3 3 
OT a.Tr 

call refaze 
CALL RK4(NF:S3LV) 

NGOOD • 4 



LABEL PCM 


Page 
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8!> 



RETURN 

86 

C 



_ 82 

— c 

, 

4NT.£flRATjON- flOHPLETE UITMIN EBBQH. BOiWDS-, 

a? 


500 

CALU SOLVCSf 

8? 



DO 510 UliNF 

9J 



DPSfl) 




DPCU) * "PC<h 




F(lil) » r<J»2) 

.—93 



. rj-J..«24-a_E’t4j3) . ... . 

94 



FU«3> « Fn»4) 

95 



F(I|4) 3.ra,5) 

96 


510 

continue 

92 



NGOOa a Mooqo+l 

98 

c 



99 

- 


4F_{l2Bi!4l_..tja. a) CAU TPiJlAmS*-!! 




RETURN 

sol 



end • 




* 


PASS 


- U8EL- PCM 


3 - 
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Label »data Hiat JL 


1 

2 

3 

4 

5 

6 
1 
8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


e 

E 


B 


C 


subboutine bcata . . . 

reads rocket data asd writes dn a tehporaru Tape, 

IRTEGER ORDfft _ . 

real MSER 

logical iPRt , . 

CHaBACTEB NAHE*89 
CHaBACTER L1H|E»69,CQLa1(80} 

DJhEnSJON Fl99},!i(9?rim?T)»XH99) 

COHMQN/SToBil!/- ZB&OiME ,P UqPlt.HElilPUCF(»l, 

COMHOR/STORPI/ Ge,MiQHECArGM«RE|RqfE5nieSQl*E6Ql. 

COMMOR/STOR83/ PBESOi0BNSOrTBMPB4'39UNOD<f?BB5Bi?)fTEHPB{?>» 

. HALT849>eTSi;OPet9!»gQNST,PBESiDBNStTEHft*SOUNp 

COHMON/STORU/ KPRB ViHPRBVf MPBE V, aBRGV^LPRBVUaJ ^VPRBW tl4 ) , 

oaQBML<i 

COHMOR/SToRIW/ lTAREi;iTAPi^lTARB5nTAPB4.l^RlNT, fpUNCHi INDEX 
COmmON/SToR??/ NAMBiAaRA?<IPOl , „ 

common NPSTrtJSBP;NSp;ii|M{14],KSP(20l;PMASS,f3O|TAV(5)‘,'4tJG(6t)(f8(4i j 
V PLM<?OJTPl!lTt88hTOER(92UHSlP{3a)TT8P<?Q)tSPM«80lfA<9T,14) 
data Sl.CG,S-l!{T/2»*0‘,/ 

EOUIVALENCE iITAEE2»iaKlLTAPe5ilHi — 

EQUIVALENCE f ABRAYIH I ARRATUOO ) i^> »t AR«A v (199? »F1 )V 

. (ARRAytzfei.xn 


I B NPMAXt2 
DO 1 J»ifl 
parij? f 1 j8E?id 
TSEP(J? 9 l,OE*10 
MSEP(J> * 0, 

1 continue 

DO a JxlrNSMAX 
TSP(J) « l.OB-^iO 

2 CONTINUE 

IPRa * IPRINT .EB, B 

RBADdR^ioooi Name 
REaD{IR»1OO0) lire 

READdRilOOli RMAS31T3D 

READdRilOOPi LIRE . . 

READ tIRi loot) RpgTKPHTlMnnrNPSTl 
READ(IRit009) LIRE 



ORIGINAi; PAGE 
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LABEL SdATA '(•AGE Z 


A 3 

44 

45 

47 

48 

49 

90 


-91 

52 

33 

94 

95 

96 


-97 


60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
76 
79 

ao 

ei 

62 

83 

84 


HiAP{?RJ»OOii BsepJ(T88Ptlt,lni;N96P) 1. 

RlAPtlBilOOl} USE ^ 

RlA5l»iiS00|} NSBp:tHBBPUhI»liNSBp> , , 

RKAti. ■ EHAiS*CF<5] „ . . 

DO So !«1.'N|8P 

* M8KPl>I>*i:f«SJ ..... 

‘io coNfjBue ^ , .... ... 

IE(VN0T,lPR8> 60 T0 11 .... 

PRINT 9(i0o . . t 

PftiNT 9001. !>HAS g.T8P ^ 

PRINT S0B2i NPST*(P«Tm»J«t?SPST> _ 

PRINT .9003, flSBpi<T5g|<I|.t»*rNBEP) . .._ 

PRINT 9004, NSSpJcMOBPnlitiliNSBPj ... „ ... 

11 CONTINUE - ..... ... — 

B. PBGuu.R.fi«4ais.; — : : 1: , zz: — ~ - 

D? lOl NR«lfNP9T . ■ 1 

IP<lPR?l PRINT 50l0t ktP .. 

RSAOdRiJgoBj USE - 

AfiD.J 04 - : 

Dfi 110 J*l,14 . . 

RiABaRilOOfll USE 

IFU.RE.aj #0 TO 10* . . 

DEc6DB(LfNEiitO02> COL . . 

DO 100 111, 70 , .... 

IE tBllUll.,Eai.'SLU4UlitDLlI.« 11*50., iiy.*AHfliJ3aiJU.*2U£fliliJJ , = 

• ADD 8 BU)9HBES0 ... - , _ . ... 

too continue . . , 

104 RBAB(IR;1d 01J HAR,tA<IiJ>»t»llKAX) . ... ... . ■ 

IB(«AK,GT,1J GO to l09 ... 

Mix '4 2 ... - - . . — ... 

A<2xJl-a.-ltliJ.i . 

105 NM<J) 8 MAX .. - . . 

IF^J.UT.t ;oe. J.fiT,8k QD TO ISO ... 

IR<»IAX,NB,2J go TO.tlQ . * » - -■ — - — ■ ~ 

IF<Alli7>;i,TiZ8Rt( .AND, A ( 2, 7 Ji LT: ZERO > Go TO 107 , ...... 

106 IB(4(ti6.J.LtiZER0 ’.AND, A < 2,8 ) .LT'.'ZERo > SO TO 108 __ 

GO TD 110 , ... 

107 Ail,7» V SLCO ... .. 

A<2i7» • SlCfi ■ . . . . 



3 
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label'h'dat* "'raoe 


89 

86 

i08 

17 

le 

no 

2? 

90 

91 

92 

93 

94 


95 

96 

97 

98 

115 

99 

8 

100 


101 

102 

103 

104 

120 

105 

106 

107 

108 
109 

(30 

110 

111 

112 

113 

114 

140 

115 

116 

150 

117 

lie 

160 

119 

120 
121 
122 
123 

170 

124 

125 

126 

ISO 


GO TO 10* 

A«2i8» f SlJt 

contisue : _ 

si;c6 t MNM)ijr);7'i 
SIY¥ « A(NMtil>,e) 

DO il> ikiiCiSAR " ' 

A<I»01> t A«.0i> « eHt(NPj 
A<lt02) « AiT«02) * »PD 

A<Ii03) .^.AUijjjtSFiajL , 

A<I| 06 > * At!, 06 ) ? PHT{N?r 

A<I»1B) * AH.lpU'BP^U . 

A(I»18J » ADtil2)#Bhl> 7 

continue 

INTERPOLAtiOKS BBGJr; 

lF(NM<*)tEgy.;Ai Di_AB8(«<j;i6)»Ma«*M.Lr.7EFl0l GO Tfl 1*0 
DO 129 lAl.tjPAO 
X<I) « AHifit 

XKU « Kill) ■ ■ 

continue 

DO 150 J»7,9 1 

HAX ? NHLJ)-- ' ' __ J J 

DO 130 I*1»HAX 
FK) ■ AUiiji 

continue 

call LJNT(FirX.Fi;x'!iOBDEN<B>*HA)(;i:NMjfl)fii’;NH(i>» 

MAX ■= NMU) 

DO. 140 I*1l«1X _ ’ 

A<iij) • ml) "■ 

continue 

NH(J> » NMIM 

continue 

GO TO 180 

coNTiNye . . 

MAX « NMU) 

DO 170 1*3, H»x 
A<I»7) « A<117j 
A(I,B) « AltJSj 

continue 

NM(7) " NHUl " 

NN(6) ■ nh<11 

continue ■■ 
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LAB6LJD«li- . PAGE i 


127 

126 

129 


130 

131 

132 

133 

134 

135 

136 


137 

136 

139 

140 

141 

142 

143 


144 

145 

146 

147 
146 
149 
no 


191 

192 

193 


134 

l55 


196 

137 

196 

199 

160 

161 

162 

163 

164 

165 

166 
167 
166 


DO 230 K*9«t9>2 
= K 4 { 

IP(NK<K) .E0f9 >aNd. ABS(AaiKWA<2'iK>>^Lt,7ER0) GO TO 210 

DO i9_0„iiljSSAft 

F4U f Allitj) 

xnj • AUiKl 

XKI) ■ *(!»») 

190 CONTIWyE . . ■ 

CALL.jaHTi.rfX> P liXllBQOEW ( 3 > KJJUL/Mt 4 JJ. 

MAX ff NH{4) 

DO 800 niiHAX . ... 

A<IiJJ f Flin . 

200 continue 

NH(J) 9 NM{in 

GO 10.238 

210 continue 

MAX * NMt4> 

DO 820 I*3iHAX ... . . . 

A<IiJ» I AdiJI 

220 continue 

NMldl .6 NtUa J 

230 CONTINUE 

C INTERPpLATieNS eHd‘. 

NRlTEdHl (eili.lBltfll 
IFCIPR2> PRINT 5011i 
DO 890 J*lil4 

MAX » NMU) . . . 

WRITE<IW5 HAX>tA<i;j»;i»l>t1AX) . . 

IFCIPR2> print 5912, NP> JtMAXfl A< l*J > , 1*1, MAX) 

290 CONTINUE 
500 CONTINUE 
C 

C SPENT stages; .... 

NSp f 0 
400 CONTINUE 

ReAD<IR<lOQO,END945D) LINE 
NSP » NSP+1 

IF(IPR2I PRINT 5520# NSP 

READdRitiOOl LIHE . 

READ(lR;iOOt} TSPlNSPI.SPMCNSPdBfS) 

srminspi » gPMUNSpipens) .... ‘ M 
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1«9 

170 

171 
122 . 

- 173 
' 124 

175 

176 
i27 
170 
179 
160 
lai 
182 
163 
184 

165 

166 
187 
166 
169 

190 

191 

192 

193 

194 

195 

l9? 

196 

199 

200 
201 
2Q2 
203 


WRlTEitW) B»3> 

iR(tPR2) PRINT 5021i TSp«NSP»^SPMtNSP)r'BT3> 

00 416 J«4||| 

REAlll^lOgj } Mp,, , .. 

REAB<IRitOOtJ HAiOlMli Jl*t»liM*xT 
IKCMAXiQy;!) Q8 TO 4B0 
MAX 0 2 

A<2iJ) ! AC^Zli 
405 NH(J) ■ MAX 

aBimiHt_HAli. U<I; j>7i.i |7MAI<> 1 

1P{IP62> PRINT 5912; N5Pij7«*XMA»K.Jj;i*i;MAXJ 
410 coNTiNue ■ ■ '■ 

GO TO 400 
450 continue 

B • ' 

1900 FORHAT(A8o) 

1801 FORMATJVJ ' “ ~~ 

5000 F ^RMATaBlJitHiNeUt data S8T NO, APfeR' intBrpolatjons to 

; get. ABROpYtlAMlC DATA REFESR8D TO 5N6'trHE vARIASIB AND ONE MaCH N 

Vo; variable,} 

5001 E0RMAT (/#1 4M RBcReT MASS »rlPEl5.5*i2M7 ’SURN-OUT TIHB 4 i1PE1375) 

5002 F0R«AT</»19a NDmIe^ Sp PHASES **left?H, START 't i 

• aoiisEis.s}}) 

5003 FBR«AT</;25I;I NOmBeS of spent STaSeS d,i2»19H» START TIMES 
;aoaREi3,5)i> 

5004 F0RMAT</#25H NUMBER OF SPENT STASaB «,i2»l5H» THfllR MASSES*,/, 

• UOaPElS^B) ^ 

5010 F 0 RMAT{//, 18 M PH*sS NO. j" ’ 

5011 F0RMAT{/i9H ARSaT S*i4(tPES3,5) » 

5012 F0RMAT{/H9H ARRAY A,,NP;j?MAX";siS;/, iia(ipBl3,S) 1) 

5020 F0RMAT(//,16H spent STAQS NOV "ViaSlSH i*— 

5021 format </,l4(|l TSp;sR.*i;B«3l»;/;3tlPEi3,5j j 

return 

END ■* “ 
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labeu UeEaze" 


PAGE 



1 

2 

3 

4 

5 

6 
7 

e 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
IZ 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


e 

9 


s 


SUBBtOUTJNI 8EFA3E . 

UBOATING THE DATA AT PHASE CHANOB-i” 


' U^lW|yBAA4.^ 

COHHON/STOROi/ 


laiEMR BKST|EiOMa8... ^ ^ . . 

REAL MASS,HCl0P,MPLiHp0T»MACHlWiUCP> lYK.MSEp 
pouBtE pneciKTCiti DtsiEPc ... „ ^ 

C0MHON/STOR57/ T;DTtT6|Pt;Dfl.AiD?A ! » TBoiRKSTE?,EPBIG, EPTINT , 

DHl^iDBAXiK. 


. .^»Stl4?iDPClil4)»6RH<14jV|P(6>20» 

COMMON/ST'dRlSA MAsSt HBB'nP.HPLfSMBMtsat.MP of .THRtH ACHiBAiADtQS. 

; su)gEisHR;LBa;LU!; 

COMRON/SToRtt/’ KPRBV,HPRfiV.yMP8aV.,APRBV.fLPREV<14),VPREVtl4|i 
; 0BDBR114> . 

COMMOS/SToRSa/ NPHA*;y?BA.X£NeLttiANIllJyAR,NP.NeiAa.SPENtiJS,jSB»h)cON 
CBHMOR/STQRtjl/^ I T APS 1 (61 i MU 

: PUB(2D)tPHTJ521iTBEP(32);MSBP<321.rTBPt2O}iSPHi20hA(997l4) " 


NR P NP * 1 
NPl = NP * 1 . 

TNPt *, PHT.<6JR11 

R 8 AJlti.TAaEJJ_ (fill lIi.ali. 4 Li 

B«31»B<3}y2.kP 

B<41 * Bt4)*W 

DQ iO.JfltU 

PO i I*tvNyil8 

f 0i 


1 CONTINUE 

REAatITAPE3> MAX;iMr;j>»lPl»HAX). 
NN<J) « HAX 
LRrEV<J) * i 
VPREViJ) f AUtJ) 

10 CONTINUE 

MDOI « 0,._ 

SOBM(NP) 9 MASS - .S<i;3> . 

S<71 • A»liSJ*S.Y 

S<8» ♦ AUi91*SZ . . 

DBS*7MSi7J . . . 

DPSi0l"SI0) 


JSl s JS*.l . . 

IB(JSl.QT;NeMAK*£) RETURN 
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label refazI page 


43 IF(ABS(T«TSBP<asi}5,9T*ltO0-O5> QO TO 10# 

44 JS n JS ♦ 1 

45 SUBM(NP) s SUBM(ldp) w MsBPtJSl) 

46 100 DO igO 

47 iB(iBStT-TSP'{Ji)' ^Fr-k “iioe-o5)"QQ To iao 

48 JSC 0 jse*l 

49 DO ilO |*1|6 

so sp(iijsc) B S(n 

51 tlO continue 

?g 120 CONTINUE _ 

53 return 

54 END 
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LABEL *K.4 ^AflE 1 


1 

2 

3 

4 

■‘S' 

6 

7 

e 

9 

.IQ. 

U 

12 

is 

16 

17 

le 

19 

20 
2i 
22 .. 

23 

24 

25 

26 
27 

,2B 

29 


32 

33 

34 

35 

36 

37 
36 
89 

40 

41 

42 


e 

B 

e 


8 


SUbBDUTINE RK4eNP.S0Li/>... . . . 

integration by RUNiB«KUTTA-S!MPSON MBTHO® Of FOURTH ORDER, 

RALSTON eoeeFlClENVS ARE USED FOR MlN|HUM 'ERR oR BOyNDS',’ 


INTEeSR SkSTBP ‘ ■ ■■ ■ 

double precision DRgiDPC^DBF 

double precision Gl<t4)«e2tl4|[»fi3(li>»i4flij 

COmMON/SYORO?/ TlDYiTniRlXDTLAjBfAliTBOZlfRSTiP/EPBtS.’EPtlNV, 

'cOMKQN/.5TQRfi|^ Stl»>°? *i?!rlr DR s??4ffnPC^^^ 

data AlV o', 4/ ■ - *’■ .... 

data 8^02,8?/ 0HSSrj726» .Qk2Y6t776.e,...ClSl875966/ . 

p“IL‘ o.iruw,.; 


CALL.SOyUll . 

DO 150 JTI.RKSTEP .. 

TT « T . 

DO 110 1*1, Nr .. ... .. 

Qim p OTPFlIlJl „ - 

S<ll • DPSdJ * ilPOtfll . 

110 continue .. 

T B TT # AleDT 

CALL S0UVI5I 
DO 120 1*1, NR 

G8(I) a DT«Fni5J . 

son 9 DPS(i) * 82 »aim * b3*Q2(iJ . . 

120 continue.. - 

T B TT ♦ BlfPT 
CALL 80LV(5) 

DO 130 III, NR . 

G3(I1 P„5T»Fm5) , . 

S(II * DRSdl ♦ QlRStlD - C2*02M| * 05*.Q.Siitj 

130 CONTINUE.- _ . .. . 

T e TT * DT . 

CAUL SQLV(9J ... 

DO 140 1*1, NF . 

G4<n n 9T*ni»S> . . 

s<n • opsni * Di*5im - d2*q2<i 1 + 03isi(lj -+ D4*a<un 

DPsni .B ani 

. DBcn» * 0, 

140 continue _ 
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Label ^k4 ' 


PAGE 


43 !P(8K3TEP,ESU) 30 TO 150 

44 - JJ A J * 1 

45 CALL SOLVCJU} 

46 150 CONTINUE 

47 return 

46 END 
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L abIOoIket I^AGE 1 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 . 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


SUBiOUTIf^E BOCKET(KODE) 

B interpolates betwebn thrust and maCh tables, 

B 

LOGICAL KODB .... .... , . 

integer rkstep 

real «ASS,MeB0P,!4PLiHDPT»MACH,LCQ,LcPrIY¥,M5Ep 
COMMON/STORO?/ PBESOIdBnSOjTEMPO i 8OUNDO|PRBS0<9 )fTEMPB t9>» 

;• HALtBI»},TSi;OPB<9>,BDNST,PRSS*DENS,TEMpiSOUNO 

COHMON/STOR05/ VT(5>IVA ( 5) ZVELf ACCiVOLIMfVUl” 

COHMON/STOR07/ T;DT,TNPliDTLA*D.TAIij:Bl).RK.aiE?jiE£.ei5j.6RTiijy, 

• DHlf4«DHAX;KSTEPrNQO0D 

COMMON/STOR09/ M*SS,MBRoP,MPLjSUBM?32( ;MDQT,ThR,HACHj6*iAD,QS, 

*' SLOB8iCHp»LeQ;L5PilTYJSV«S2 . ...^ 

common NPST|NSEp;NSP;NM|l4J,KSPt20)»RMASS,T3B|TAy,.(3JiA>lO.<61),B(4)f 
•' PL^^<20);PHT^5^llT^EB(3^^iMSEP<3a)TTSp(2fl)t3PMU6)^A(9»^l4> 

e ' ‘ 

C INTERPOLATION BETHEgN HaCH NO TABLES, 

CA n 0, 

IFISOUND ,ESV 0.» fiO TO 12* 

HACH i VEL/SOUliiD 

CALL RTEPP(HACM,J.Si*liO«?120l . 

Jt B J * 1 „. 

FRAC • (HACy * A(JI4)l/(A(t)l^41.; . . 

CA • AIJj5» • FR*CiMA(Jli5>-AtUi5Jl . .... 

IRIKODEI GO TO l50 

CMP B Auaoj * PRA5»fA{jizio)-A(j;ig)) 

SLOPE « A(gi,i2> * rRAG»{A(|ji;i8>.Aej«i e)> 

LOP « A<J,14} * PRAC«(A(Jlrl4>«Atj;i4>T 
GO TO 190 
110 CA B A(li5> 

IF(KODE> GO TO 190 " 

CMP B A(t,10i 

SLOPE » A (1/12} ._. , L_ 

LOP B A(t,14J ... ■ _• 

GO TO 150 . ;; ; 

120 CA B A<NH(5J,55 

121 IF(KODE) GO TO 150 . ’ „ 

CMP * A{NH(10>,10> ... 

SLOPE « A(NH{12)*12) ._ 

LOP = A(NM<14),14) 

3 j;'"‘ 
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43 

44 

45 

46 

47 
46 
49 
90 


93 

94 

95 

96 

97 
. 90 

99 

SO 

61 

62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 


c INT6RPOU4T10N BEtHllBR TlHE TABliBS; 

150 IP(AeS(AUil}>A(2;i)lt 2DO*220r2dO 
200 CALL PtBRp<TiJri;$2l«.‘S220j 

Jl * *,1 . 

DELM * A\^l*5)-Af48J 
FRAB • (TsAlilil) iAbeLf 

tHRiA<J«2)4FRAC*1At«li|2)»AtJ*2H-Btl>*^RBS 
lE{tHR.Lt;o;} 'ThIiO,' 

MRrSP ■ A<sJtS> + F8AB*DELH 

intKODB) RSTUHEl 

MDOT • delh»PEL;t 

LOG B A(d,7| ♦ r0ASMA(Jl»7>*A<Ji7h . . .. 

I¥Y i A{d,0} ♦ F8AB»U{jlna>-A4d»8»J „ 

ftiTURN 

8iO TWR » 

MRR8P » 4(1(3) 

!P(K008) RETURB 
KDor « 0) 

LOG I 4(1,71 

lYY I A(1,8) ; _ 

eiTURBL- 

2C0 TNR B 0, 

MPR6P B o; 

IR(KODB) RETURB 
MOOT » 0, 

LOG • A(NH(7J|7) 

I?Y 8 

return 

END 

1470 EQUAUITY or NON-gOUALitS BOMPARTSObl NO'T B6 HE*NIN(Srt)L 
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^3302 01 11-2987S 18,150 " ' LABEL StERP ' PAGE " ' "i 

1 SUBBOUTINE BTERP(ViL;1 ,*♦*} 

2 C IRTERPOLATtOM PROGBtH CALLED RV RoCKET. 

4 ISTE08R OROBR .. ' 

5 RIAL HSEP 

6 • COmHON/SToRU/ KPRev,HPREViMPBBV,APREV,LPREV(l4)»VPREvei4l, 

7 , 0SDERU4) 

8 CBhMON NPSTt^SEP*N5P',NH(i4),KSP(2ei.lMASS,f3D,TAV(3)7A(J6(61)fB(4)« 

9 , PLM«20)jPHT(S2).TSER(32l*H3eP<32>tTSP»20>»SPHi2ai»A(99,l4i 

JO c . . 

11 Lb LPREV(l) 

12 MAX = NM(IJ 

13 IP(V,0T,A(L,1> .and, V.LB,Aa*lJ{ ) 1 8o TO ZOO 

14 IP(V,LE,A(l,t> J BETURN 1 

15 ip(v,QE.A(«AJt*i>> Return z 

18 . IF(V - VPREvnn Itl 0 , 20 o;il 0 . 

17 iOC JB » 2 

16 JE P L + 1 

19 GO TO 120 

20 110 JB p L *’ 1 

21 JE ■ MAX 

22 120 DO 130 J»JB,JE _ 

23 JSAVE = J 

2,4 IF(V , LT. AJJiU) QO To lAO 

25 130 continue 

26 140 L B JSAVE - 1 

27 fOO LPREVin = L 

28 VPrEVU) » V . . . _ 

29 return 

30 end 
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LABEL (FtImAT "PAQE . 1., . 


1 

2 

3 

A 

5 

6 
7 
6 
9 

IQ 

11 

12 

13 

14 

15 
19 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


Q 

8 


8 


StfBiOWTINE 

ROTATION hath* PrOa TopOCBNTBiH.«Bo»CEN> TO INERTIAL, 


DnuBLE PftECjHbN D gSJDPc ftr YX^yg 
c8HMOiil/Sfom/ G?,5 rX0ME5X7W«I7 


rUVBsq 7B so T » E s q 1" 


C0MMOlil/StOR(r5/ x8,’SD?20VR,S)?T|SLT^SR;HLrfl' 

C0hM0NPS70R|I|> rTU?> 

COMMOS/SToRfll/ sU#irf(i4i5>J15P8rii»)"rDPCltl4)lERWU4)fsPt6,20l 


X • 3.11) , . _ 

Y e «2) 

Z 0 S<3) 

Q B X«X ♦ Yif 
R B DSORTIO. I Ut) 

RKy I bSORTljfl) 

AUT ■ R «■ R i ; 

R{. $ liOOO/fi . 

R*Yl ? iJO.DOifRRY/ .. ' 

RTI14) f -Sl?)»R«Yt 

RT149) * SU)*BXXI 

RTU6) « 0. . 

RIU7) S_s£llilfiL 

RTUB) t -SU)nRl 
RTI49) A -SD5>«Bt, 

• RTia) » RTt {95*9Tt<5k 
RT1<2) ■« -RTtlf JARtjH) 

RTU3) » RX¥»Rl 

RBrypR 

END 
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LABEL SETIV PAGE 


1 

2 

3 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 
19 
t 9 
17 

le 

19 

20 
21 
22 

23 

24 

25 
2« 
27 
2B 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 


e 

e 

6 

6 

C 

D 

e 

Jt 

6 


SUBROUTINE SETIV 

initiali?es all the necessary data bepore starting phase no IV 
THIS routine initializes step dt st Integrating first by rk4 with 
PT=PTUA One step AUP then with DTeDTliA/l TW_Q STEPS. *ND 
cohparinS at t»t*dT InTeSRation vaRiablIs obtained DIpFERENTLYV 
IF error Lies betoNc the error bounds a Nem Step dt is itqraYbd; 
NCON $ f>l PIKES LAUNCH AND INERTIaL ORIQIns. 
k 6 TAKES LAUNCH IN GD ERAHgt 

• 1 takes launch in ac prahe.’ 

♦ 3. ONLY ASsOrBS THAT. WIND rS P rB sBnt, .. 

integer BKSTEP.OBDaR 

real LATL.UONGL-,LATtUONB»HAS3jHPRoP,MPL.HOQTiHACH,LCG,LCP,lYt»MSEP 
double precision DRS'iDPCfARllARS 
EXTERNAL SCLV5D 

DIHENSIQN. SA9EU41.. . 

COmHON/STOROP/ ZERDi:0NE,Pl?0PIrHPr»TPl.CF(9) 

COMMON/sroROl/ G«.6R;0HESAiGHVRB.Ra,ESd|Rsci| »ESQ1 
COHMON/STORC2/ PaESO.OENSOrTBMPO.SOUNDO,l*RBS0<9) .TEHPBt9», 
i H*LTB49».TSLOPET9I»OoN8T,PRES*DeNS,TEMPJSOUND 

COhMON/SToROS/ lifl.VoVZO.R.RXT'iALTtSRiniJlH 

COMMON/SToR 6»/ GDLATLVLATLirLOHOL.rHE TAL.PHlLlBDLAT.l ATrLON S.THgTA. 
♦' PHlJBiMAL;PSIL|GAMAfP9UELQOUA?QDL,ELGO.A?QD(iBL, 

I az.sraz 

CGhHON/STOROS/ Vt{5»;VA<3) jVEL.ACO.vgLIMfVl’l” 

COHMON/STOROe^ RTI{?| 

COHMON/STOR07/ TjOTiTNPtlDTLAiDTAl iJTBOiRKSTEP.BPBIG, EPTINY. 
DHlNlDHAXtK3T8|*]ieoOD 

COmMON/STORCB/ Sfl4);ni4.3)*DRl(t4);DPCTl4>lERWll4‘>TsP(6.2or" 
C0MM0N/ST0R09/ MASS.MPR0P.MPL.9 UBm{ 32J rMOof.THRiHACHTcA*AOiQ3, 
SLORSlCMP^LCGJLDP.lTYiSV.SZ 

COHHON/STORie/ WE.HN .MUASIrWlblD^WNAZlALOtt, aHIQIS.NLEV , ALEVt lOO J , 
i' ALTW(l0O>;vpLEtl0O) kVRLN(lOO) 

COMMDN/SToRII/ KeReYiUPRBV»MeBEY,APRE.ViLeBESf<ll).yPREVtl4), 
i OBDERUA) 

CBmMON/SToRIS/ NPMAX*NSMAXjNPLiNANG*NVARjNP|Npi, JSPENTi'jS.JSq.NCON 
CDhM0N^SToR 14/ lTAPfil?ITAPE2;iTAPE3,lTApi4;iPRlNT.lPUNCH,INDfiX 
common NPST|NSEP;N5p;NM(t4)|KSP<20l,RMASS,f3U|TAV(3>U(jQ<61).B<4 j. 

. PLH<20),PHTt62J iTSEP(92)iMsep(3BUTSP<20>.SPMl2{l>.A(9f ,14) 

IFINCONJ 10,20.30 “ - -- - - • 
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-13 

44 

45 

46 

47 
4B 
49 
30 
51 
32 

53 

54 

55 

36 

37 

38 

59 

60 
61 
62 

63 

64 

65 

66 
67 
66 
69 


72 

73 

74 

75 

76 

77 

78 

79 

80 
6l 
«2 

83 

84 


10 continue 

e Fix I.AUNOH AND IQERTIaU ORIGINS, 

G ENTER here IP LAT LONG ARE CHaNgEd; 

ARl f GDUATL 

AR2 « ESQUtlBSlNlAftlJ/DCOSCARlM 
ARl « DATAN{AR2| 

LATL • ARl 

TSETAL T -LATU ♦ HRI 
SK ■ DS|N(Ae$) 

CN a DCOS(AQti 

phil»lonol 

IF(LONOl;LT,0.} PHTU*TPI*I-0NGL 
ARl * LONGL 
SU > OSlN(ARX) 

CL ■ DCOS(AB)f> 

BBS 0 OP|,ATL«LATL . . , . 

ARl » Bps 
SEPS ■ DSlNlXRi) 

CEPS » DCOSIARI) 

RE P RQ«SoRTtE8Ql/UTO«E80»CM«CH)i 
XO n RE*CNtt6i: 

TO » RE«Ch*BL 

ZD p REpSh 

1F(1PR!NT;nB4'2> RETURN 
PRINT 5000. WCbN 

PRINT 5081, (SDLaTLILATL, THETAL.EPS 
PRINT 5002, LONGL, PHtU 
PRINT 5003, UejiX0,TO,zu _ 

return 

e 

20 continue 

G INPUT lAUNCy A&QLES ARE IN QSOxDGTtc FAaHe! 

G enter here ir Launch angles aRE changed, 


SEL c DS!N(A«1) 
cel ■ DCQSUEli 
ARl » AZ8DL 
SAZ » DSlNIAjlj 
CAZ ■ OCOS(*ill 

ARiPSeLf5EEflaSB?5*e8LilCAL„ _ 

AR2pDSQRT{1«EDQ6ARI»A|U> 

GAMA * DATAdiABl/AlSZI 



gAllj£‘ 
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LABEL SETIV 




"i 


85 

86 
87 
8B 

89 

90 

91 
9' 2' 

93 

94 

95 

96 

97 
9B 
99 

100 

101 

192 

103 

104 

105 

106 

107 

108 


111 

112 

113 

114 

115 

116 
117 
116 

119 

120 
121 
122 

123 

124 

125 

126 


ARiaCEL* 8 AZ 

ASaiiClPS^CELXCAZ^SBRSsSEL 
PSI = DATANZtAftlJASe) 

IP(PSJ,L'T, 0 ; j PSlcPSlATpI 

EL • SAmA 
AZ » PSI 
ELGD ■ ELGDL 
AIGD » AZGDL 
5 

30 continue 

c enter here if lat, long, and launch angles are not CHaNQEP, 

EL p GAHA . . 

AZ « PSI 

aRi = gaha 
ST P OSIN(ABl) 

CT p DCqS(ABI) . 

ABl = PSI - - 

$U P DSIN(A 81 ) 

CU P DCOSIABl) 

ARI s EL 

SEL e DSINIAPIJ 

CEl = DCQS(ABl) _ 

ARI = AZ . . „ 

SAZ p DSINIARII 

CAZ = DCOSI 8 P 1 ) 

ARI = SEPSpC/r*CU ♦ CBPSpST 

AR 2 = l.ODO 6 ARl »»2 

ELGD = HPI . 

IF(AR 2 ,LT, 1 , 0 D« 15 > QO TO 31 
AR 2 = DSORUAR 2 ) 

ELgD « DATAatAEl/AR 2 i 

31 ARl s CT*SU 
AR 2 PCEPS#CT#CU-SSPSp 8 T 

AZgDpDATAN 21 AR 1 #AR 2 } 

IFTAZQD.LT.Oi ) AZGDsAZGD*TRI .. 

B 

I . . 

0 INITIALIZE THE INTEGRATION VARIABLES ANP PARAMETERS, 

c iterate for good initial STEP dt; 

DT P DTLA _ 

RKSTEP * 1 - 

110 IFIRKSTEP .fJE, 2 » GO TO 129 



page 



label sE^iV Page 


4 


«33ba bi' ie«i94 


127 

126 

129 

13P 

131 

132 


193 

134 

185 

l36 

.137 

138 


139 

140 

141 
j42 

143 

144 

145 

146 

147 
146 
149 
iJo 


Hi 

H2 

153 

l?4 

195 


1*6 

197 

196 


199 

160 

161 

162 

163 

164 

165 


166 

167 

168 


DT R DT/tH 

8 

120 Nf) R 1 

NRi « N.P,* 1. 

T o TAV<1) 

T8Pi « PMT<S(»1} 

TBO » T5D 
AtT *,TAV(2) 

R 6 RE 0 ALT 

y el R T AV f8 1 J. IPIb gLlLT’i VOLlHt .V E L _] 

VSq » vei;r*8 

UNcO, ) WS»0- 

RBWINO tTAPBS ■ 

RIAD(ITAPE3J (B(|),'lRl*4> 

B<3JrQ< 3)/2,0 » Bt4jR0<4}**2 

PO 42’ .. 

DO 124 UliNPAR 
A<ItJl 9 0< 

124 continue 

READ(tTARE3) HAX;{AaiU}»I4i;MAX) 

B HAX 

UeREVtJ) ? i. . . . 

VPRBV<J> » AH»J» 

125 continue 

HOOT « o; ^ . 

SUenUNP) R BKASS « pPL » Atl»3> 

HBREV * 1 

HPR5V *0* 

KPREV » 1 
APRIV ■ Oi 

K8TEP«iO I NSQSDrI 

J9b8 I JSCrO 

• INDEX *0 . . . 

S<1) » BtCMSPL . 

5(2) k R*CMr8L 

5(31 • R8SM , , 

VIbVEL*CEL»CAZ 

va = VEL*CEL*8»Z _ 

Vla»VEL*S6t „ 

call BTIMAT-- - 

S(4I •« *TI(il*91 * I!!I{4I*92 ♦ PTl lT) *V3 - QHBSAOSte 
S(5J 4 PT1(8J*V1 * pTK5)*V2 ♦ ♦ oMbgA*8U> 
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169 

170 

171 

172 

173 

174 

175 

176 

177 
179 

179 

180 

181 
182 

183 

184 

185 

186 

197 

198 

189 

190 

191 

192 

193 

194 

195 

196 


199 

200 
201 
202 

203 

204 

205 

206 
207 

209 

209 

21c 


S(6l « RTI(3>*91 + 9ri<6MV2 + RTIt9)*V3 
S<7J 9 0, 

S<8) 4 0| 

S<09> » •Cl,*6U*cr*6M n CT*9l,«SU * flH«CL*ST 

S<lO) ■ iCU#Bf«SL*5H 4 CL*CT«SU * Ch«SL*§T 

S<li) 1 Om*CU*Bt + 8M«3T 

S<12) • Cl»SM»IU = CU«9L 

S(l3) ■ 9 u«SH#9U + CUkCu 

S(14) ■ »CM#»U 

DO t3P I6t,i.4 

DBsiii ■ sn I 
130 CONTINUE 

IPCRKSTEP ;go, 3) OO TO 900 
call RK4I14j»0LV9dJ 
IROKSTgP .Bo, 2J eo TO 150 
DO 140 :«lit4 
SAVE(I) • SID 

140 continue 

RKSTER 9 2 
GO TO IIP 
190 EKAX ■ 0 i 

DO 160 t*l«18 
ERR s 0, 

IF(ABS(SU> j;GT;i;ea»0«) ERR*ABS(14 0 "SAVEm )/'SI 
ENAX > AKAXlUHAM.eRR) 

160 CONTINUE 

ENAX « ENAX/15*0 

IR( lPRINT;Eai2i PRINT 5Q06* DT*2 ,EMaX,SAVeVs 
IFiBMAX b 1,06-05) lTO.l70rl80 
170 RKSTEP « 3 
GO TO 118 

190 IF(DT - DHINJ 170 , 170, *190 
190 RKSTEP 9 1 
GO TO 110 

900 IE( IPRINT,NB;2} RETURN 
PRINT 50lg, NCON 
PRINT 5011, ILRDL.OAMA.EI 
PRINT 5042, A?RDL,8SU'AZ 
PRINT 5pl3, T,DT;Nf»iRKSTeP 
PRINT 5014, 5 


label 9ETIV PAGE 


5 




♦3302 01 11-25575 


ie;i54 
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211 

C 

212 

5000 

213 

5001 

21A . 

50Q2 

215 

5003 

210 

5005 

217 


216 


219 

5010 

220 

2Z1 

soil 

222 

5012 

223 

5013 

224 

5014 

225 


226 



F0RMAT(lMl,i5HTHB 8«IQ1N AT THE LAUnCh SITE !s FIXEO. NCONiniai 
F0RMAT(/»25H GDLATL»LATUITI*ETAL»BPS »»4< 1PE13,3 > > 

FqRMA_I(/!»i3H l.qN5LiFHIL_?*2(lREtg;^n,. 

format(/^*iah Re,)iiqiYO,zg «i A(1 Pei 3»5) j 

FOR«AT(/t' JTEJIATJOR FOR good step size — DT,EHAX«'»2B?0,tge/, 
S WITH DT ONE STEP AND S W|TH DT/2 TVio SfgRS.'i/j 
.(toElS.Sr/iUBlS.S)} 

FBrmAT<///» 70H TNE VBytckE 

.isitializeoa hieoNiiiSi 

rORMAT(/rl6H EUDLTGAKAiEI 
F0r«AT{/'t 15H AZGDLjPSI.Ai i >3llPEl3 ,5 ) ) 

F0RflAT</»l7H T.'DTfliiPjftXSTER •> lPEl3,5,lPEia;9, 2I>) 
rOR»AT(/*l3H SU...|t«> »|?»<10<1PB13JS>)) 
return 

END 


ATTITUDE AND INTEORaTION HARlAEUES AiRE 
'■f3(l>Ei3i5}>‘ 
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LA86L” S0LV3D PAGE 



1 

2 

3 

A 

5 

6 

7 

8 
9 

in 

11 

12 

13 

14 

15 

16 
17 
16 
19 
2f> 
21 
22 
23 
2d 
26 
26 
37 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 


SUBBOUTINE 3DL*30(KJ 

G EQUATIGNS of Monous IN 3 dimensions, 

e 

REAL HASS,«eROP,MpLfHDt)T;HACH,LCQ,LCP. lY4l*HS6p 
DOUBLE PREClilSN dPSJDPC 

COhMON/STOPOI/ QO,ORIQHEGA;QM«ME,RQ,FSQ,SS»IiESQZ 
COHMOH/STOR02/ P9ESp , DEnSO 7 TBHP 0 , SOUNDO , PRES'-' (9 > # TEmPB 1 9 J » 

, hal^bi? >.tslope(9», CONST, Pres»d6ns, TEMP, Sound 

COHMON/STOR03/ XO, Y 0 i ZO , R, SXY , ALT , SRiHL iH 
COMHON/STOR05/ VT (5 > , VA < 5 ) rVEL , ACC , VOL IM.f Vlj” 

COMMON/STOROS/ RTitpl 

CDmhON/SToROS/ S{l4>iF(l4,5>,DPS(l4),DPCU4)<8ftW<l4)?aP{«»20> 
COHMON/STOR09/ MASS|HPROP,MPL|9U0M(32 )i‘hDoT|ThR,MACH#6A|ADiOS» 

, SL0R5»CMP*LCG,L0P,iVy,SY,3Z . 

COmMON/STORIO/ WB,BN;MLAST|WlND*MNAZtALOri,AHlQH.NLEv;ALEVUO0), 

, ALTHtlOO)*VBLBtlOOIiV^LN(lOO} 

COHMON/SToRia/ NPmAX;nSmAX,-NPU,NAN6,NVAR»NP|^PI, JSRENTIJ1S»J^«NC0N 
common NPST,NSep;NSp;NM(14},KSP<20},RMASS,T30,TAV(3);AN8(6i>»8e4>» 
, PLMt20),PHTtS2J,TS£P(32UMSEP(32)fTgP»20)»SPM<flS>»A<9f*14> 

C 

call HTIHaT 
call atmsph 

IF(NCON,EQ,2j call WTERP 
GR p QM/tR»E) 

51 a S(4) 4 . 0MB6A»S(2t 

52 a S(5) - OMBGAaSUi 

VT(1> a RTlUJaSt * RTI{2>"S2 ♦ RTI(3)*St6) • UN 

VT<8> a RTI(4)»S1 * R?I{3)6S2 + RTl<^»aSH6i ■ WE 

VT(3> a RTUT)»S1 X RTIl«)«S2 ♦ RTl(9)*Sl6i " - 

VSO p VT<i>a»2 * VT<2»**2 * VTt3)**2 

VEl a SQPT(VBq) 

OS a DENS»VSO»B(3) 

call rocket true, J 
HAsSaSUBH( NP j+MPBOP 

WTI s 1,0/MASS " . — - . .. 

ad c CAeOS 

AF a THR - ad 

SI a AF/VFL 

AFCl c Sl»VTtl5 

AFC2 a S4*VTt2} 

AFC3 a Sl*VT?3l ■ " ’ 

F(1»K) a 5(9} 
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Label soLvsd' page 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IS 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
41' 
42 


subroutine gOU¥5D(K) 

C EQUATIONS OP MOTIONS IN 9 DIMENSIONS, 

C 

real HASS,MeftOP,HPU»MDOT*MACH,LCG,LCl**lYT|HSEF/NA,NB;LSK,LEXC(}, 

I • JDCON 

bOueiE precision DP5,DPC»ARlrAR2.Vi;y2.V3 
COmMON/SToROI/ G9,GRiQHEGA«QM,RE.RQ,ESO,ESnliESQl 
C0 mM0N/ST0R92/ PRE50,DEnS0?TEMP0 ,SOUNDO,PRSS8(93«TEMP8<9> , 

. HAlT0(9>,TSLOPE49),OoNST,PRSS*D5NS,TEMP,SOUND 

COhMON/STCROS/ XO,TO;ZO,R,RXY, AUT.SHiHLjH 
COhMON/STOROS/ VT(S),VA(S>, VELiACC,VdLi«rVLlM 
CDM«ON/SToRP»/ RTK9) 

COmMON/STOROR/ S«l4>»F(l4,5),BPS(l4>iDPCU4>«ERW{14);sP(6,20J 
COhMON/STDRPP/ HASS,MPRoP,HPL,SUBHt32»^MDOT, 'HR.HACHVCA, AD,03, 

, SLOPEiCMP»LCB,LDP,I?YiSY,SZ . 

COhMON/STORIO/ WE;RN,MLAST;W!NO;WNAZiALOM,AHlGH,NLEVi ALEV(IOO), 

, AUTWaOtOiVELEtlOoUvGLNCiOq) 

COhMON/SToRU/ NPHAX;NSHAX,'NPL»NANdl,NVAR»NP,''Pli JSPENT,JS, JSC.NCON 
common NPST»NSEp;NS'*JNH(l4),KSP<20j,RMASS,T30,TAV(3)VANG<61)fB(4j> 
, PLM<20>',PHT(52) ,TSEP(32)iMSSp<32)»TSP(20)»SPHI20),A{9f,i4J 

EQUIVALENCE (S^9 1 ,pMli ) , (S (10 ) ,RM8l> i <S(il5 ,Rn 31> , (S t 12 ) ,RH12 > , 

; eSH3)IRM22)*(S<14)iRH32) 

C 

RMi3 » RH21*RM32 - RM31«RM22 

RM23 » RM31»RM12 - 4M11*RM32 

RM33 s RH11»RM22 > RM21«RH12 

CALL RTIMat 
call ATHSPH 

IF(NC0N,E0,2J call WTERP 
OR » aH/(R»6> 

WEC e HE ♦ OHESA»RXV 
VI ■ S<4) - RTI(1>«4N - RTt(4)*WEC 

V2 ■ S<5» - RTI(2)4WN - RTI(5)*wec 

V3 « S<6) - RTI(3>»4N - RT1(6)*WEC 

VA(i) B RKlliVl + R42t»v2 ♦ RH3i«v3 

VA{2) B Rm12#v1 + BM22*v2 + Rm52*v3 

VA(3) ■ Rk13»V1 * R423*v2 ♦ RH33*v3 

VT(3) = RTM7)*V1 4 RTI(8)»V2 + rtI(9)*V3 

VSQBVAd J*«2+VA(2)«#2*VA<3)**2 
VEL"SORT(VSe} 

OS=DBNS*VSO»0(3) '■ ■ 

call RQCKETt'.FALSEi) 
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M*Sfis8UBM(N9}+HPSa9 

MTlKli’O/'HASS 

ADsOAtOS 

AR ■__mi . 

ARl ? VAI3>^VAU> 

ACFA * DA.TAl!ltAftl) . „ 

ARl » VAl2)/9Atl) . 

BtT* • DATAtJtARl) 

NA=5l.ePE.#ALE:i*l5S ' ~ ' 

N0=SUie£lflmAOS 

LExeGtABS<ai8>»LBCl 

JQC6N • HdOTSURXBQ . „ 

PDcfiER « 0;?hCHPfQS|Bf4)/veU ; JP80N#kBX?IQ 

JDcftN « JOCe|N**.fl ■ "■ ^ J 

TER(i >f 1,0/lSY 1 .‘J 

s* 8 iizuiaatL ' 

SI K KBhTERH 

AS a NB«OBC0 H*Sz ' ‘ ““ r.11' 

A« n NA*JDCDN*!SY 

LSH»ABS(UCQ»LCR) „ ^"1" _11_ 

F< i,\i) * SiA) . . „ 

Ft a+Kl ft- SIS). , 

Ft 3<K) f St9) , __ 

Ft 4,KJ a QRkRTnZ) A WTt«{RHi*aA?^RHia#*Q-RFIl3*AH) 
Ft S,K) * GSYRTJtBj * Wj!*tRM2J,aAP«M^'»*0«RR23*AHi 
Ft 8,K) • Ga«RTH9} ♦ WTl*tRM31»AP'^RM3?*AQ^RW33*AHT' 

Ft 7,K) a «tiA*l.SH ft PpCoBFaSt ' ” 

Ft 8.8) 1 fcll*l;Sl! ■ gMoeFas z 

Ft 9.K) • RMlZaSt 5 BH13«SY 

FtlOiKJ a RM22*S? » «H 23 "SY 

Ftll.K) * RM52fSl • Im 33*SY ' ■Z;iT". ' ' 

Ftia.K) • -8Mii*az 

Ftl3,K> a -BH21*BZ _ ' 

Ftl4,S) a -BM3ia3z 

return 

end 
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1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

13 

13 

1* 

15 

16 

17 

18 
19 

30 

31 
33 

33 

34 

35 

36 

37 

38 

39 

30 

31 
33 

33 

34 

35 

36 

37 

38 

39 
4(1 
41 
43 


C 

0 


c 


SUBROUTINE SPENTUBOTI 

calculation of tbajectort for the spent stages, 


INTEGER RKSTEPiOBDER 

RIAL LATL. LflNGL.L At. long. MASS, HPR oP.MpL.flDnf* MACH, LCS, lCP.it?. I PT» 
; IPAI.IpR.MSeP 
double precision DPS;DPC*AB1*AR3 
CIaBACTER NAME#08 
EXTEPNAt SOLV3D,90LV?0 
dimension SAVEI28,in 

COhHON/SToROO/ zero, one, P lrQPI*HPl .TPI,CF(9) 

COHMON/SToROI/ g8, GR. omega. GH. re, RO.E sn.EsnUESQt 
COMMON/S7OR03/ PBESo.OEnSO, ‘TEMPO, SOUNDO. PRES' I9> .TEHPB 1 9) , 

; halts ( 9 ),TSLOPE<9),DONST, pres, DENS, TEMP, SOUND 

COhMON/STOROS/ XO,YOi'ZO,R,RXYi AUT.SRiHLlH 

COhHON/SToROA/ GDLATL,'LATLiLOWGL,THETAL,PHILl6&LAT,LATtL0NS^HETA, 

, PHi;GAMAL»PStL»GAMA,PSl,ELGnL<AZQDL,ELOD,AZGO,EL, 

, AZ.SRAZ 

COHMPN/STOR05/' VT(3 ) iVA (3 ) ,VEL. ACC,'VOlIM»VLIM 

cohhon/sToros/’ rti 

C0HMON/STORO7/ TjoTj TNP iIDTUA, DTA 1 , TBO, RKSTE'* , EPS I G , EPT I NY, 

I DMIM,PMAX»KSTEP,NGOOD _ .. 

COMHON/SToROB/ Sf l4)iF<l4,5J,DPS«i4),DPC<14>»ErtH(l4 J,SP(6,3of 
C0HM0N/STOR09/ HASS , HPROP , «PL , SUBm < 33 ) , MDO t , ThR, MACH.'CA , AD,' Q3, 

; slope;cmp*lcg;lcp,iyy,sv,sz 

COmMCN/STORIO/ we, l!|*i;MLAST*WlND*WNAZi alow, AHIoH.NLEV; ALEV (tool, 

, ALTW(100UVBLE(100)I:V^LN(100> 

COmHON/STqRII/ KPREV,HPREVfHPRGV,APREV,LPREV<14),VPREV<14), 

, 0RDERU4) 

C0hM0N/SToR 13/ NPMAX»NSHAX;NPL,NAMG,NVARfNP,Npi, JSPENT.JS iJSGiNCON 
C0HM0N/STOR13/ APTJ APAZ.APA, APR;BoT.BoA,9oR,BoVL,BOEL, 0OA2, IPT, 
IPAZilPR 

C0HM0N/STOR14/ I T APel.,’ I T APEZ. IT APe 3 , 1 TAPE4; iPR I NT, I PUNCH, INDEX 
C0 h« 0N/ST0R?9> NAMEiARRAY(600} 

common NPST»NSEp;NSp;NM(l4),KSP(30t ,RMASS,T3P,TAY(3) I ANG<61)»B{4r» 

i • PLM(30)-,pHTfS2>.TSeP(32UMSBP<32)f.TSP*30>,SPM(30),A(99.14» 

EQUIVALENCE t ARRA Y t ) I ,'SAVE (1*1 M 


SAvEO « omega 

iraRQT.eo.Ol OMEQAfJi;0E"10 
REWIND 1TAPE4 
DO 700 IP=i,NPL 
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43 

44 

45 

46 

47 
40 
49 

90 

91 

92 

93 
34 
95 

56 

57 
9S 
99 
60 
61 
62 

63 

64 

65 

66 
67 
60 

69 

70 

71 

72 

73 

74 

75 

76 

77 
76 

79 

80 
61 
82 

83 

84 


m.sl 

PRINT. i.QOO .. 

PRlHl iQ>Oi NAMEiHPt4G0,AZQDU*CF<i) 

■ifx f LjsTtefl'ig) P8 i[ii im. .. 

IE(l6i,ilSTliol8)_p|i&T 1012» Wlljtf*WN4 


iBCNLtVyST'.il >rInT 1013 


WlBtf;WN4Z«CF(i>7A[.QW.AH!GH 


8BTUaB._„, 

io IFtlRllNtiiaiOi-RRlST 1020 
DD la6.IltlLl)AS&_ _ 
ECgBL.p AiNfiltAJ... - 


Dfl-Jtal_KtljJllE . 

Bill no; . . 

Bt 2 J- l.Di ... _ 

e<4» 8_ii; 

P I AD i W AP.E 4 } . J 8 P I K J iSP.H { K ) f B 1 3 ) 

6131 9 B19J^3j8. ... 

QB ^ . 

M»X i 2, . 

00 ao J«ti 8 l»AR 

AiliJl 10 . 

20 CBNtlSUE 

IE(J,t;T ,4 ;a 0 i..Jffl.T.«U. 00 . TO 21 . 
READi.r.TAREi 1 . .nixStM 

21 NN(J> ■. MAX. 

UBRev<JJ » 1 . _ ... 

VBREV(J') > AdfJl 
SO ceNTIBUE . .. 

REaB(«TAPE 4 ) ( 9 (l)n>l* 6 ) 

DB 40 . 111 * 6 . 

DRSIIL I S(l) 

40 CBNTIBUE 
NR 4 . SPHKXaI 
NRl f> HP ♦ t 
J 5 i N.P . 

Jlc JtJJSHAX _ _ _ _ 

T-«.T.lP_l!U ... "T 

TBPli _l‘_t 8 T<tlllL..„ . 
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LABEL spent 


85 

DT B DTLA 


86 

SUDM(NP> = 3PM»K) 


87 

KPREV s 1 


88 

HPREV « I4aLT0(1) 


89 

mrr3v » 1 


90 

APREV = ZERO 


91 

RKSTEP * 3 


92 

KSTEP B 10 


93 

NQOQD ■ 1 


94 

CALL RTIMAT 


95 

SI ■ RE/R 


96 

call AZRAN(St»S(l>rSl»S(Z>iSl#S<3),Tl 


97 

WE B 0« 


98 

WN ■ 0 , 


99 

IR(NC0N,Eq;2} gall wTgRp 


IQO 

SI « + DMgGA*3{2> 


101 

SB • S(5) - DMEGA»Stlj 


102 

VTd) = RtUD»S 1 ♦ RfI(2)*S2 t RT1<3 

}*S16 j 

103 

VTta> B RTI\4)«S1 ■* RTI{5)«S2 + RTM* 

)»S16i 

104 

VT(J» a RTI(7)#S1 *. RTI{8)*S2 t RTK9 

)aS(6i 

105 

VSQ » VTa»»"2 ♦ VT(21»*2 * VT<3»#*2 


106 

VEl « SQHTfVSQ} 


107 

ARl a -VT{3) 


108 

AR2« SDRT(VTtlJ*«2+VTj2)«*2> 


109 

ELV = DATANlARt/AR2> 


110 

ARl B VT(2) 


111 

AR2 « VT«1) 


112 

AiH s DATAN2tA81rAg2) 


113 

1F(A2M,LT.0,> AZHbAZM+TPI 


114 

SAVE(K.1J a TSP(K) 


115 

SAvE(Ki2) a SPH<K)*Se 


116 

SAVE(Ki 3) a ALT 


117 

sAvE(Ki4) a SR«cr(4) 


118 

SAvE(Ki5> a SRAZ«CE(1| 


119 

SAvE{K« 6J a VEL 


120 

SAvE<Ki7J a 6LV»6FtH 


181 

SAVE<K(6J a AZH»CF(1I 


122 

0 


123 

index = 0 


124 

IF(lPRINT.Ne^,Pl BrTvT 1017, K»ELGDL4CF<t) 

125 

call RK4(6,S0L93D) 


126 

50 call PCH(6,gaLV3D) 
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LAB6L tPSNT " WT 


127 

126 

129 

1?9 

131 

132 

133 
184 

135 

136 

137 

138 

139 

140 

141 
1<2 

143 

144 

145 

146 

147 
146 
149 
t»0 

I9l 

152 

1J3 

l!4 

115 

156 

iSe 

159 

16a 

161 

162 

163 

164 

165 

166 
167 
166 


call IPSTEP{SP50,«6tl| 

60 call aK4«6#SOL93D) 

1R(ALT,LT.HI,1M> 00 TO 70 
call l'FST|P(S60/i66} 

7b CAUL tRCAL" 

SAvE(t(i9) « T 
SAvE(K«10)« KR«CP(») 

SAvECKill)« SraZMC2<1) 

500 CONtjUUe 

00 T0 509 ^ 

PI I NT 1001 
print 102Q 

505 DO KAiisip,.., , , : .. . 

PaiNT 1030, 3LiDL*erU>,K,fSAVE<K»t)Il*l*ll) 

510 CQNTiBUE 

000 continue 

56b continue ' ■■■ ■ ■ “ ■ " ■ __ 

ONEGA f SAvgfi ; ';r . " 1 ■ ‘„ri. „ 

'lOOO FORMAT<lHl,J»HNAaA WA'UUOPS FU16NT BE(rfB fff77 29H WinTlTOTS ITliJtm); vjl4~ 
,’GPNlAi/'/»37H tiiAJEeTBRY 0U«HAR¥ p68 SPENr¥TA&ESi7/71 " ' ‘ 

lOOl FOR«A.r.{lHlJ _ _ 

1010 FORMATUSh vehicle I ;a 80V/-,J8H “ftp LOAD LOS tPARBN 

Vt vEHJCL 6>,^;13H UUSC'H AZ » fR7;a?*lT'BES' 1,l)7/r ■ ■ ' 

1011 F0RMAT{13H aiND " - --- - 

1812 F0rMAT<13h WlNb 
;i8x»5« FeoMve9ri;6»i 

1813 F0eHAUl3H..at_NL 

1914 F0RMAT<13h earth , -- 

1815 F0RMAT(13h iXRTH i I'lAlilRClTATiNfl flffBEtiT 

1816 F0RMAT,(7/7,19H MisSjNfl I'NPOT DATA;?/J[4iH 08 a 3‘ TABCl NBT BtfPPtflD’F 

,•00 Spent stijoeI) '■ “ 

1017 FORMAT(lMi;55HOeVAl|,feD prtlNT>*OUT dft SpENfafAiS'e i ,T27HH ' fllRBCfSP 

7YW1.3ZH NHEgl Parent VEHICL E L AUNcW IL tiir'«TF7~.f.’4Hr5ira74';; 

;54H UNITS A§Sli5i-S4BEQRlEr EXCEPT RANOEtfH) AND ABcmr/.P) 

1020 FBrMAT(43x, 20H3 E 'P A R A t } o' N749xU6er“ S ft 'A C f 
t.A,80X*‘77(lHP),iexj23UH*»!r/t7H LAyNeK;g,5HNolr“6j<MHTfMr,SXf6^^^ 
;G«T»6X,3MALT*7X,9HRA8QE,6X,-2HAZ;6x;3Rvein;6X,9HrLVli:74ftr6T(n:T7*z7'' 


, i4HZER0l ^ 

I ,r7i2;0H Ffi!sfic;rPr^S*l3H -DBS' •' 

Pf TO’.SPYiilH'DTT"""' 1" — ” ■■ 

. 1 .3pWVABiABLni BtRb SffP ft LiEV BY Qaew't • ~ ‘ ’ 
I •leeNON-ROTiTlNG MOfiBC) 


1030 


• un I ' Ar 7 ni«wcfUC»o Aj<nj|«. 4 OX#' 9 nvcuv 0 Xi VnrUT/'SL««A»enr LT/AZ* 

,’l6x*4NTIME,?S,9HfiAft8i;6x?21llAZj/;6H ELTDE8y7l2'x75H<seCinr*7lHl8CM^ 
I . «xt.4|(ETii i* . iHltlM )76x;5HtDEfi):3xS8HtFf?sBet;ix;SH7Ba0n5>[.fi-rtniffl 
;>ilOX#5|:Useq],4x;4Jlil?BJ.6X*5HiDl5ll72l„. 

FORMAT<F3,2,l5#Fll,2,fl0,2*F10,9,Plti2,Fi0l2_*F12,2pjFtO,¥»Fl![J'27" 
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LABEL SPESt PAg¥ 


l69 

,F9»a»F10.2J 

170 

RETURN 

I7l 

END 



LABEL "splash PaQE 
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1 

2 

3 

4 
6 
6 
7 
6 
9 

10 

11 

12 

13 

14 

15 
18 

17 

18 

19 

20 
21 
22 

23 

24 

25 


5 

e 


5 


SUbBOUTINE splash ' ■ 

CALOULATES TRAyECTBBIES POR WIND EFlpEcTsi 


INtECER «KST8P _ _ _1_ :: 

extIcrnaL sou<> 3 d, Solved 

C0HHON/STORO3/ X8,So:i!0,R,eX¥>ALT,SRjHLlff 

cohmon/storSS/ vt(8jJvA<8)>-veL»Acc;voUMfVCiS,„ 

COmMON/SToRO?/ T;D'fjTljPl»DTLAiDTAl<,TB0rRKSirEP,6PBIS,eRTIH¥, 
V DMlNiDMAX«KSTeeiNiSOllD 


call RK4fl4fSOLV5D) 

100 call PCMI 14 »S 0 LV 5 D) 

IR(T-TBO> lOOfllO.'llO 
UD KStEP b 6 
200 call RCMC6»gOL¥3D) 

iF(Vt<2j ;lt5o,oi lift TQ.aoe 

K $ T £ P s 2 

300 call PCM«6iSBLv3D) 

call IPSTEP{$ 30 o;S 409 J 
430 call RK4(6iSOL03D) 

IF(ALT .LT. HLIHJ SO fO 900 

call JPS!EP«S4ft0;j«051 . 1 , 

900 call IPCAl 
RETURN 
END 
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UABEt-TMINT PAOE - - I ■ 


? ' C . 

4- -c. 

s 

6 

7 

■8 

—9 

iq 

u 

1 ? 

u 

Vi. 

18 

12 

18 

19 

23 

24 ■■ 

2 ? 

24 

28 

29 

26 

22 

28 C 

29 

30 

31 

3 ? 

33 - 

38 

35 

36 

37 

38 C 

39 
<10 
81 
8 ? 


SU3WOUT1NS T><HINT{KtL) TPRlNT 

detailst pumt-out or the trajectory. 


1NTE3ER KKSTGP 

real UATL,U3Hau.l,AT,l.Q*lfi,HASS,HPHnPit1PL,MD0T,'MACH,LCG,LCP,IYY,MSEP 
DQgai.E PRECISION DPS,DPC,ARi.AR2 
dimension CSP(2) 

COMHON/SIORBU/' ZER0,0NE,PI,QPI,HP| ,T!JJ,CP09) 

Mwaoil/Slosai/ SO-»^6H-.OMEOA^GMtSEjaQ+eSQ,ESOi;ES01 

COHHON/STORA3/ XQ , Y 0 • ZO , H , RXY , A^. T , SR iHL I H 

COMMON/STOR98/ RDUATL, LATLi LONGL* THEIAI, PH IL;gdUI*. 1.AT ^LQUQ^ THETA, , . 
. PHI.QAMAL.PSIU,GAMA*P3l,EI,GDL;AZflnL,ELQD,A2QDtEl, 

. AZ.SRAZ 

COMHON/STnR0!J/ VT ( 3 1 , V A ( 3 ) , VEl. , ACC* V8L I H. VLl H 

...CUHMOkl/SIOBIB&X STi-(ji) . . 

COMHON/STORO-A/ T , DT , ThiPi , dTua , DT AI , TSO. RKSTEP , EPBIQ* EPT 1 NY, 

, . DMIN.QHAX,«StEP,NQOQD 

COHMON/SIORB^/ SU8).F(14.5),DPSU'l)<;PPC«18),ERHllA)*SP{6,2Q^ 
COMMON/STOR09/ MASS , IIPROP , HPU , SUUM ( 32 ) , MDOT , THR.HACH* CA* AO, US, 

. SLIPE.CMP.UCG.LCP, lYYiSY.SZ 

- common/ STor 1,1/ WE,UN.*MLASl,WlNU,lilNAZ;AtOW,AHiaii,ULEV-aLEia44laJ, 

. ALTw(i0O),VEuE<10q),VBLNUO01 

COHMON/SToRta/ NPhaX,NSMAX,N®L,NANQ*NVAR,MP,NP1,JSPENT*JS, JSC, NOON 
C0MMQN/SI0R48/ I TaPEI, ITAPE g, I T aPS3, ITAPEA, IPRXNT, JPUNCH, INDEX 
common ''BST,-;iS£P.NSP,NM{14> ,8SP{3BJ .3MASS;T3n,TAV{3I* ANQ<61J,B(8J, , . 

. PUM<20) ,PHT(32I ,TSEP(32) irtSEPt32) ,TSP(20) *SPM<20UAIS9,U)„.. 

-]W.TA-Ill«e«,c3P/ Q, lH0«-4.i(l/ . 

CALL FXoeT<71, 1*1.0) 

iniNDEX ,EQ. 0) PRINT lUOQ* CSPJIPR4NT) 

I o T 

IF(1 .EO, 01 no To 1 

-IF(L- .ec,,.!) GO To.l . , . 

IF<I .E3-, IBRev) return 

1 IPREU s T 

INDEX 3 INDEX ♦ 1 

IF(M0DnN!)EX*.3a) .EO.D) PRINT IDOU* CSP(2) 

-ElA-a 0, . _ 

IF<T ;GT, TROl GO TO 15 

ARl 3 S0RT(UA{2 )*o2 + VA(3)»«2> 

AR2 a ABStVACl)) 
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LABEl. TPRImT 


.44. .. 
_ 

M 

46 

— ja. - 
46 

M. .. 

.. . S4- - 

54 — 

• 52 
52- 
.* 54 
-.„55- 
-56 -- 

U 

SB 

49 

.. .. 68. . 

61- 

. 62- . 

64 

. - 64 - 
.. .65 . . 
66 
6Z 
6S 

63 

7Q 

- 74- . 
-72- - 

74. . 

74 

— n — 

76 

77 

78 

79 

- aa 

82 

S3 


ETA = DaTAN-(4«?J/AR2) 

Sl3S(4)*0MES8aS(2) 

SSsSiSJ-QMEQdttSili _ - - 

VTUl = «TH1)«S1 ♦ MT1(2)*SZ ♦ HTK2)*S(6) ■ WN 
. VT(2J 5 RT1M)*S1 * RTl(5itS2 * BTI{fi)#S-{6J c HE 

. V7«31 s BTl<7JttSl * HTI(81*S2 S HTHSJ«S<61 

- 45 VT£1) « VT14> « Hi! • .. . . 

. . «T<21 s VT(21 * WC 

SI. 3-VT{Wj»*2..*. VLT4-2-7**2 .. 

vei.o = <;6 ht</t( 3)«»«2 * si» 

ARl s -V;T<3) 

AR2 = S,-78T(3i) 

EUR tt DAIAN/ AH1/AR2) 

ARl » VT<2J . . - • — 

AR2 » Vt(4> 

AZl « UATA.^2tARl* AR2> 
lF(A2l,tT,o;i AZUaZIaTPI 
SI » RE79 

call AZR4N(S4«S(1)iS4*S(ZJ*S1»S(3>,H , . . 

U a 0.5fcQ*5/8i3) ' ... 

OAUEA-b 1) . — ^ 

ACC * 0, . .... 

3 a -2 

Do 20 Ix4>3 . . .. - 

JajtJ 

Si B RTI(J)o('’( 4,K1 + RTI {J*U*E{b,Kl a .RTl(J*27tF.(6«K} - 

, ... 10l!EGA«92)o(RILC.JJ«tSilj4RXUJ±llilit2j) ..i. - 

, atOBOHESAoCRTl t J*l)*SU>-RTl£J!«St5n 
ACC. R ACC St»o2 ..... 

20 continue , . . . _ 

ACC R S3BT(AeCI 

PRINT 1340, T,SRaCF(4>,SflAZ«CFlli,AL?;ELE*CFll),AZI*CFlU4.VeLD*. 

... , ACC:7aO,.!tACH»ruR,AD40^£tAAC|;UJ*4AU:4,HASSJACO,WtUVlE 

laoo FORMAT! Al.! TI'^E RANGE BEARING AUIlTUDE FLT/EL FL.T/A2 VEL 

,0C1TY ACC OACH THRUST DRAG. D/P3ES- - EJA QAuFA HEIGHT MIN. , 

.. ,DN WlNDfil.A/) . . . 

4310 FOHMAT(l*,FI*2,2Fa,2.F10,Oi2F8«.2»FiiIia4lXt+E642iLEfl«fl,J’7.a»Fa+04. 

„i;7.^24.F7^0,£a*l.,.2F5il4- - 

RETURN . . . „ 

END 


Page 2 
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LABEL (JCONST PAGE 


1 

2 

3 

4 

5 

6 
7 
B 
9 

m 

11 

12 

13 

14 

15 

16 
17 
16 
19 


G 

e 


e 


subroutine WCOfiST 
constant wind table, 

COMMON/STOROO/ ZBR0,0NE,PIfQPl*HPl»TPliCF(9) 

COhMON/STOPIO/ WE,HN»MUAST,'WIND*WNAZi4L0N, AHlGR.RLEV7ALEV{idin, 
ALTWdOO >»VBLBClllO ) i' VELNUOD) 


WAZ r AMOD(WNAZ+ei, TPl ) 


ALTH(t) 

ALTW{2) 

ALtU(3) 

VElE(I) 

VElE(2) 

VELE{3) 

velmu 

VElN{2» 

VElN(3) 

return 

END 


O.l 


ALOW 
ALOW * 

AHIGH 
0 

WINDoSiNtWAZ) 

VELEt2) 

0 

WIND«COSCWAZ) 

VELNf 2 ) 



LABEL rtTERp' (‘AQE 


1 


'93302 01 U*25b75 


iail72 


1 

2 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
26 
26 

27 

28 
29 
3n 

31 

32 

33 

34 

35 

36 


e 

? 


8 


SUBtOUTINE HTERP 

INTERBOLATEB north ftND EAST WIND cOHPS FROM ilTABLEj 


INTEOfiR ORDB# 

C0m«OS/SToR 95/ Xf,?0,?0,R.RXV,.AI,T^SR;HLlH , 

COHMON^StORtO/ WB,rtiq;HLA9TfWINDiWNAUALOll,AHlO&|NLEViA.LEVtlOfi5i . 

: AUSil<*()Oj»VeL,EtlflO)rVgLN(lo3j _. . 

COMMDN/SToRU/ KeftBV»HPREV;MPBEV,4RREvafREV(l4).tVPREVtl4h 
; q8D8M14> 

EOuIVALENCE.jtAlJLiAL 


H = MPREV . , . 

IRCA.QT.ALTHtM} ..AtiCt A.t.Ei'ALTH(M+lH Gq TO ?00 

IF<A ;gt, AU.TW^MLAST) J QO to 220 .. 

'.LT, ALTWHJ) SB TO 2?0 . 

1F(A * APREVl. iafliS.14;ilO_. 

lOD J8 > 2 

JB » M*1 .... 

GO TO 120 .. . . 

110 JB » N*1 „ .. 

Jfi ■ MLAST _ 

120 DO t3Q . . .... 

JSAVE * J . . 

IFCA .'LT, ALTWfJI) SO TO ISO _ 

130 continue 

140 M : JSAVE -1 . . 

ZOO HBrEV » M 

ABREV 4 a . .._ • . 

FRaB “ tA-ALTWtMhT(ALTWtH«U-ALTMfMh 
WE • VELE(H) * FBAe*JVELE{H*l)-VELS(»1H 
WN • VELN(M) + FaAC»IVELN<N*l}-VELN(Nn 
ZIO RETURN 

220 WE * 0. . 

WH n 0. .. 

RETURN 

END _ . 




